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[Text]  Chapter  1.  Definitions  and  Initial  Formulas 

1.1.  Some  Definitions  and  Formulas  for  Determinate  Functions 


The  determinate,  i.e.,  random,  functions  or  temporal  processes  3(t)  examined 
in  the  present  section  serve  as  the  analytical  representation  of  radar 
signals. 

We  shall  employ  the  symbol  to  designate  a  pair  of  Fourier  transforms, 

e.g. 


(1.1) 


which  in  expanded  form  means 


^(0=  -J 


(1.2) 


The  letters  t  and  f  designate  time  and  frequency.  Here  and  below  we  shall 
not  stipulate  the  existence  of  improper  integrals,  since  it  is  assumed 
that  all  of  the  integrals  of  the  determinate  functions  examined  are 
absolutely  convergent. 
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The  real  temporal- process  3(t),  representing  high  frequency  oscillation 
with  slowly  varying  amplitude  S(t)  and  phase  2'n'0(t) 

3  (f)  =  S  (0  cos  (2*/,/  -j-  2it0  (0  4-  <p}  (1.3) 

can  be  written  in  the  form 

or  in  sjmibolic  form 

where  fQ  —  carrier  frequency,  (p  —  initial  phase,  S(t)  complex 
amplitude (envelope) 

S(0  =  S(0e'''''('\  '  (1-6) 

which  determines  the  form  of  the  signal  or  the  nature  of  the  modulation. 
Accordingly,  the  function  S(t)  can  also  be  called  the  modulation  function. 

If  the  spectral  density  of  the  function  S(t),  called  the  modulation 
spectrum  below,  is  G(f) 

G(f)^S(0.  ..  •  (1.7) 

the  frequency  spectrum  of  the  process  3(t)  and  its  S3rmbolic  representation 
3(t)  will  accordingly  be 

j  ('W  s  a) = T  o  (f  -  M  ®  (-/-  f.)  =■'’ 

i(0->i{/)  =  O(/-/,)e'’.  (1.9) 
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The  overscribed  bar  designates  complex  conjugate  quantities.  When  the 
carrier  frequency  is  significantly  greater  than  the  modulation  spectrum 
width,  the  frequency  spectrum  of  the  symbolic  representation  imitates  the 
frequency  spectrum  of  the  original  on  the  positive  half-axis  f  (to  within 
a  constant  of  2). 

Frequent  use  is  made  below  of  the  scalar  product 


J  '  (1.10) 

—00 

of  the  two  oscillating  processes  3^(t)  and  32 (t)  with  close  carrier 
frequencies  f^j^  and 

I  (1.11) 


Assuming  that  the  result  of  integrating  an  oscillating  process  with 
slowly  varying  amplitude  and  phase  over  a  sufficiently  large  interval  is 
practically  zero,  the  scalar  product  (1.10)  can  be  represented  as 


J  3,  {i)  3,  (0  dt  =  4-  Re  {e'  5  5,  (0  S,  (() 

t  ‘  t  *  .  —00 


(1.12) 


Following  the  definitions  introduced  in  [18,  8],  we  shall  call  the 
quantity 


(1.13) 


the  total  energy,  and  the  function 

t  (.) = ^  J  J  (0  *  (<— )  dl  =  Re  (e"-'-'  x 

—00 

^  —00 


(1.14) 
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the  correlation  function  of  the  process  3(t).  We  note  that  the  term 
correlation  function  employed  in  the  theory  of  random  processes  is 
analogously  applied  here  to  a  determinate  process  —  the  signal  3(t). 
We  can  write 


•j)  (x)  =  IV  (x)l  cos  -f  X  (x)], 


(1.15) 


where  |'F(t)|  and  K(t)  —  are  the  absolute  value  and  argument  of  the 
complex  correlation  modulation  function  '?(t) 


^  —00 


(1.16) 


The  power  spectrum  of  the  signal  a)(f)  is  associated  with  the  modulation 
power  spectrum  ^(f) 


(1.17) 


by  this  formula,  which  is  analogous  to  (1.8): 

«>(/)=lgW=xO(/-/.)+-r«(-/-/o).  (1-18) 


The  carrier  frequency  is  defined  by  the  relationship  [8] 

■J  J  S'{t)W)dt  =  0. 

•  —00  — 00 


(1.19) 


It  can  be  assumed  in  general  that  the  power  spectrum  of  a  radar  signal  is 
synmetrical  about  the  carrier  frequency  fg,  i.e.,  the  modulation  power 

spectrum  J)(f)  is  an  even  function.  Furthermore,  the  modulation  correlation 
function  'i'(T)  is  also  an  even  real  function.  Accordingly, 

x(x)  =  0for  It, 
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while  (1.15)  is  transformed  to 


(1.20) 


Consequently,  the  correlation  function  of  a  radar  signal  looks  like 
oscillation  having  only  amplitude  modulation. 


In  some  calculations  it  is  helpful  to  employ  an  approximate  representation 
of  the  modulation  correlation  function.  This  study  employs  the  rectangular 
approximation 


(1.21) 


where  T  —  correlation  time 
K 


J  I'l'WI**. 

—00 


(1.22) 


and  a  representation  based  on  expanding  the  function  'i'(T)  by  powers  of 
the  parameter  T.  If  we  introduce  the  notation 


I 


00 

W"  (0)  =  - =  4*'  =1 - . 

f  lS(0l‘d/  J 

— W 


(1.23) 


then  for  small  T 


(1.24) 


The  parameter  3^  defines  the  dispersion  of  the  frequency  components  of 

the  modulation  spectrum  ^(f)  with  respect  to  the  "center  of  gravity" 
(the  point  f=0) ,  and  can  serve  as  a  measure  of  the  spectrum  width  of 
the  signal  3(t). 
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The  signals  emitted  by  a  radar  usually  have  a  periodically  repeating 
modulation  function  which  can  be  written  in  the  form 

00  oo  k 

2  S,JI-kT,)=±  J]  (1.25) 

— 00  — 00 

It  is  assumed  that  the  function  (t)  is  identically  equal  to  zero 

outside  the  interval  (0,  T  ).  The  latter  is  stipulated  in  order  for  the 

M 

function  in  question  to  approximate  radar  signals,  but  is  not  a  bound  which 

flows  from  the  exposition  above.  The  other  quantities,  such  as  G,  ¥, 

etc.,  referring  to  the  same  repetition  period,  are  also  given  the  subscript 

T  ,  so  that 
M 


Or  «(/)->  Sr  «(0. 


(1.26) 


The  function  (1.25)  has  a  discrete  frequency  spectrum.  However,  by  using 
the  concept  of  delta-function  6 (x)  we  can  arrive  at  a  generalized 
representation  with  which  the  spectral  density  function  also  includes 
separate  spectral  lines,  including  the  constant  component  [8] 


GO 


k=—oo 


(1.27) 


Operations  with  a  delta-function  based  on  these  formal  relationships 


8  (/-/,), 

^  f{x)h{y-x)dx=J{y), 

•-<0 

^  19,  —  I 

s  S  .■=  '  ■ 

*=—00  *=—00 


(1.28) 

(1.29) 


(1.30) 
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(1.31) 


,  |]  e«"=  I;  >(i-±y 

^=“~00  J^=~0O 


can  be  carried  out  sufficiently  rigorously.  In  order  to  do  this  it  is 
necessary  only  to  carry  out  the  entire  operation  with  the  function  with 

the  function  exp  (-  ■^“) »  with  another  impulse  function  having  unit 

area,  then  passing  to  the  limit  T  -»■  0. 

The  complex  envelope  S(t)  of  a  reflected  radar  signal  usually  represents 
the  product  of  a  periodic  function  Sj.(t)  and  a  cutoff  function  l(t;T),  where 

T  is  the  effective  duration  of  the  cutoff  function 


7'=  j 
“00 


(1.32) 


We  shall  be  using  below  the  rectangular  cutoff  function 


for  |/|<0,5r  \  .  sin  nT/ 
for  |/|>o,5r)  "^ 


(1.33) 


and  the  bell-shaped  cutoff  function 

I. (/,  T)  -  |/2 r  exp (-  2,.rV*). 


(1.34) 


The  frequency  spectrum  of  the  product  is  equal  to  the  convolution  of 
the  frequency  spectra  of  the  factors,  so  that 


OO  ^ 


*=—00 


sin  nT 


X 


(^“  *m) 


(1.35) 


kT 
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and 


M 


a 


S(')=l.('.n  .J]  SrJi-kTJ^^  5]  Or,{T{) 

*=3—00  * 

:  xe.p[-2=r(/-Ay] 


X 


(1.36) 


The  following  general  rule  can  be  formulated.  When  multiplying  a 
periodic  function  and  a  cutoff  function,  the  Fourier  transform  becomes 
solid  rather  than  remaining  a  line  spectrum.  Furthermore,  each  spectral 
line  "spreads"  into  a  spectral  band  which  has  a  form  similar  to  the 
Fourier  transform  of  the  cutoff  function. 

Adjacent  spectral  bands  of  the  frequency  spectra  (1.35)  and  (1.36) 
practically  do  not  overlap.  For  example,  the  adjacent  frequency  bands 

exp|— 2i:7’*  — j  and  -  exp  |^— ,  even  with  a 

relatively  small  ratio  between  the  signal  duration  T  and  the  modulation 
period  totaling  2,  intersect  at  the  0.002  level.  We  shall  be  using 

this  fact  later  in  various  calculations.  As  an  illustration,  we  present 
here  the  calculation  of  the  correlation  function  of  the  modulation  of 
the  signal  (1.36).  Assuming  that  adjacent  spectral  bands  do  not  overlap, 
the  power  spectrum  of  the  modulation  can  be  assumed  to  be 


S  a.37) 


The  latter  expression  represents  the  convolution  of  the  functions 


and 


27’exp[— 


(1.38) 
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the  Fourier  transforms  of  which  will  accordingly  be 


and 


where 


jkss—OO 


exp 


(-??)• 


|Sr  ««)!■'''■  - 


and 


Therefore 


J  S,.(0 Sr. ('-')<"• 


—00 

I 


vM=6xp(-^) 


ik=— 00 


(1.39) 


(1.40) 


(1.41) 


(1.42) 


Analogously,  for  the  signal  (1.35)  the  correlation  function  of  the 
modulation  is 


0  '  for  \x\'^T. 


(1.43) 


1.2.  Some  Definitions  and  Formulas  for  Stationary  Random  Functions 

Of  major  importance  for  the  theory  of  stationary  random  processes  is 
Khinchin's  formula  [51],  according  to  which  the  spectral  Intensity 
function  v(f)  of  a  random  process  n(t)  and  its  correlation  function 
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(1.44) 


are  connected  by  the  pair  of  Fourier  transforms 

(1.45) 

I 

where 

a*  =  (/!*(/)),  (1.46) 

with  the  angle  brackets  <  >  here  and  below  designating  statistical 

averaging.  It  is  assumed  for  simplicity  that  the  mathematical 
expectation  of  the  process  n(t)  is  zero. 

The  reception  interference  examined  in  this  study  is  represented  by  a 
stationary  normal  process,  usually  called  a  random  process.  We  shall 
distinguish  between  two  types  of  random  processes:  high  frequency  n(t) 
and  lower  frequency  N(t).  The  spectral  Intensity  of  a  low  frequency 
random  process  N(t)  is  the  low  frequency  function  (f)  which  is  almost 
wholly  concentrated  in  the  vicinity  of  the  point  f=0.  Furthermore,  the 
origin  point  on  the  f  axis  is  defined  by  a  formula  analogous  to  (1.19). 

The  spectral  intensity  v(f)  of  a  high  frequency  random  process,  like  the 
power  spectrum  of  determinate  signals,  can  be  represented  as 

(-/-/,),  (1.47) 


where  (f)  is  a  low  frequency  function  which  is  practically  non-zero 
on  an  interval  whose  length  is  significantly  smaller  than  the  center 
frequency  f^. 


Substituting  (1.47)  in  Khinchin's  formula  (1.45),  we  find 

00  ,  -1 

(/i(/)rt(/— t)>=f  J  c^(/)cos(2«/'*)d/jcos2Tt/,x  — 
—00  ■ , 


(1.48) 
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A  random  process  can  be  represented  analytically  in  different  ways. 

We  recall  that  the  spectral  intensity  v(f)  or  its  Fourier  transform 
^n(t)n(t-T))  characterizes  a  random  process  exhaustively.  Therefore, 
the  criterion  for  the  correctness  of  a  selected  representation  can  be 
the  coincidence  of  the  spectral  intensity  V(f)  or  average  products 
(^n(t)n(t-T)>  of  the  analytical  representation  and  the  real  process. 

A  high  frequency  random  process  is  often  represented  [2,  18]  as  the  sum 

of  sinusoidal  oscillations  which  are  amplitude-modulated  by  low  frequency 

random  processes  N  (t)  and  N  (t) 
c  s 


n  (/)  =  cos  {2r.f\l  -t-  ?)  -f  A/,  (/)  sin  (2*/,/  ?) 


(1.49) 


or 


(1.50) 


where  ({)  —  arbitrary  initial  phase,  and 


(1.51) 


If  the  high  frequency  process  n(t)  has  spectral  intensity  v(f)  assigned 
by  means  of  (1.47),  the  processes  N^(t)  and  Ng(t)  must  be  such  that 


(1.52) 


<iVc  (0  iV.  (/—.)>==- (iV,  (/- x)  AT.  (0>  = 

=  J  '•V*  (/)  sin  (2*/t)  df. 

—00  • 

The  average  product  ^n(t)n(t-T)^  is  the  same  as  (1.48). 

In  accordance  with  the  real  conditions  encountered  in  practical 
applications,  we  can  then  assume  that  the  function efV®  (f)  is  even 
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Then 


and 


(1.54) 


(K  -  ■'))  =  (/  -  ^)>  =  J  (/) 


•»00 


(W,  (/)  //.  (/  - 1)>  =  (A/c  (<  -  ■')  (<))  = 

{N(t)N{t—^))  =  2{NJt)N^{t~-i)) 

{n  (0  n(t  —  t))  '=  (A/,  (/)  A/,  —  t))  cos  2«/„t. 


(1.55) 

(1.56) 

(1.57) 


(1.58) 


In  this  case,  according  to  (1.56),  the  random  processes  N  (t)  and  N  (t) 
are  statistically  independent.  ^  ® 

It  is  also  possible  to  represent  a  high  frequency  random  process  n(t) 
as 

n  (()  =  (t)  cos  -f-  *  (/)  -f-  9]  -{-  (Ox 

Xsin{2u/./-f-x{/)-f-9].  (1.59) 

The  determinate  function  K(t)  defines  the  angle  modulation  principle. 

If  the  low  frequency  random  processes  N  (t)  and  N  (t)  are  such  that 

CIC  SK 

(H..  (')  K,  ('  - ')) = <«..  (««,.('-  ’)>  = 

=cos(x(/)  — x{/— .)]  J  (1.60) 

—00 

(A(„  (f - (A/,.  (/) /y„  (/ -  x))  = 

oo  (1.61) 

=:sin[x(/)— x(f  — t)]  ^ 

—30 

the  average  product  is 

n{t  —  t))  =  cos  2«/,x  J  (f)  (1.62) 


i.e.,  the  same  as  (1.58) 
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to  a 


The  derivative  N' (t)  of  an  arbitrary  differentiable  random  process  is  also 
random  process  [22].  The  statistical  characteristics  of  the  function 
'  (t)  can  be  found  from  the  statistical  characteristics  of  the  oripnal 
function  N(t).  By  measuring  the  order  of  differentiation  and  statistical 
averaging  (integration),  we  obtain 


{N'  (0  N'  (t  - 1))  =  -  (N  (t)  N(t-  x)>, 


(1.63) 

(1.64) 


The  spectral  intensity  functional  (f)  of  the  random  process  N(t)  is  often 
wider  than  the  spectral  lines  of  other  temporal  proces^s  examined  in 
conjunction  with  random  ones.  The  spectral  intensity  gfi  (f)  is 
approximately  constant  within  the  frequency  domain  employed  in  the 
present  problem.  The  rate  at  which  the  function  j^(f)  drops  off  outside 
the  working  frequency  region  is  of  no  significance  from  the  practical 
point  of  view.  Therefore,  a  wideband  random  process  N(t)  can  be^ 
approximated  by  a  hypothetical  process  —  uncorrelated  noise  (white  noise) 
characterized  by  the  following  relationships: 


«V’(/)  =  /^,==  const, 


(1.65) 

(1.66) 


The  spectral  intensity  v(f)  of  a  corresponding- high  frequency  process 
n(t)  consisting  of  two  sidebands  —  and 

will  be  ^ 


and 


(1.67) 


(n(On(^))  =  -a 


(1.68) 
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It  should  be  emphasized  that  white  noise  cannot  be  realized,  and  that 
it  represents  only  a  convenient  mathematical  idealization.  For  example, 
a  real  normal  process  N(t)  with  correlation  time  T  interacting  on  a 

system  with  a  time  constant  significantly  larger  than  T  behaves  like 

K 

white  noise  with  spectral  intensity 


>00 

^0=  J  {N{i)N(t-x))dx. 
—00 


We  note  that  we  are  employing  representation  of  frequency  spectra  along 
the  entire  frequency  axis  from  f  =  -«>  to  f=<»,  while  only  positive 
frequencies  have  any  physical  meaning.  Therefore,  the  spectral  intensity 
which  figures  in  our  formulas,  such  as  v(f)  for  the  random  process  n(t), 
represents  a  quantity  half  as  large  as  the  real  spectral  intensity,  which 
means  the  average  noise  power  per  Hertz  of  bandwidth. 

The  following  relationships  hold  for  random  processes  which  can  be 
approximated  by  white  noise: 

0“  .  CO 

^ — 00  —00 

00 

2-  (JJjj A^c(^)^c(0^. 

—00  ^  .  ’  * 

—00 

-  — .00 

=  ^o  [fj +f{  9{t,.tM,)dt,dt,  H 

—  00*  —00 

—00 


(1.69) 


(1.70) 


(1.71) 
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where  (jiCt^.t^)  and  are  arbitrary  numerical  functions  of  the 

variable  t^,  and  t^,  t^.t^.t^.  Expressions  (1.69)  and  (1.70)  are 

obtained  by  direct  application  of  (1.66)  and  (1.56);  the  derivation  of 
formula  (1.71)  is  provided  in  [40]. 

A  random  process  can  often  be  represented  as  the  convolution  of  white 
noise  with  some  determinate  function.  Therefore,  formulas  (1.69)-(1. 71) , 
in  part,  can  be  used  to  compute  different  moments  of  random  processes. 

Let  us  assume,  for  example,  that  the  random  process  nj^(t,a)  is  the  result 

of  convolution  of  white  noise  n(t)  with  a  determinate  function  of  time  t 
and  the  parameter  a-s(t,a) 


00 

«)=  J  n(z)s{t--  2,  a)dz: 


(1.72) 


Then,  on  the  basis  of  (1.69), 


CO 


(1.73) 


If 


a)  =  Re{  jAr(2)S(/ 


“00 


then 


<A^i(/|.  “»)>== 

:=Ar.Re{  J  5  (/+/,.  a,)  a,)d^}. 


(1.74) 


(1.75) 


As  a  second  example,  let  us  find  the  mathematical  expectation  and 
dispersion  of  the  function 
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(1.76) 


'  \V{x)=±  ^n,{t)l{t,T)n,{t-z)l(t-x,T)dt. 

which  is  a  generalization  of  a  short-time  correlation  function  [57], 
Process  n^(t)  is  defined  via  (1.72),  assuming  that  s(t,a)=s(t).  Using 

(1.69)  and  (1.71),  we  obtain 

aw''w-<r  {<)»■)= 

=  x){f  (x)-i-t(x4-T).j»(x— 'c)]rfx, 


where 

00  '  ; 

=  J  !(/.  T)dt. 

—00 

x)  = 

00 

!(/.  T}\{t~x,  T)\{t-x,T){{\^x-x,  T)dt,' 

— 00 

2 

and  iIj(t)  and  Q  are  the  correlation  function  and  total  energy  of  the 
process  s(t). 

In  the  case  of  a  bell-shaped  cutoff  function 

'!'.('*)  =  exp(— 

x)=p=exp(-|^)exp(-^). 

We  shall  be  representing  various  electrical  oscillations  which  are 
certain  functions  of  time  as  points  or  vectors  in  multidimensional 
space.  The  possibility  of  this  representation  is  substantiated  in 


(1.77) 

(1.78) 
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[18,  8,  63],  and  refined  in  a  number  of  subsequent  studies,  which  are 
reviewed  in  [43].  Accordingly,  we  shall  represent  the  interference  n(t) 
as  the  point  n=(n^,n2. . . ,n^) ,  and  we  shall  write  the  density  of  the 

probability  p(n)  that  random  interference  on  the  interval  t^<_t£t2  will 

take  on  a  value  of  n  (interference  probability  density  functional)  as 
'  [18,  8,  63]. 


p(n)  =  &expj— n*(0  d(j. 


(1.79) 


Expression  (1.79)  is  derived  on  the  assumption  that  the  spectral  intensity 
of  random  interference 


1%-  for 

[O  for  i;i>/„: 


Furthermore,  the  coefficient  k  in  (1.79)  is 


k 


The  boundary  frequency  f^  can  be  as  high  as  desired.  Therefore,  expression 

(1.79)  is  used,  in  part,  for  a  white  noise  approximation  of  interference. 

It  must  only  be  kept  in  mind  that  as  f “  the  coefficient  k  ^  «>,  and 

there  exists  no  finite  limit  for  the  function  defined  by  formula  (1.79). 
However,  since  the  coefficient  k  is  independent  of  the  realization  of  the 
interference  n(t),  the  ratio  of  the  probabilit^^density  functionals  for 
two  interference  realizations,  say  n*(t)  and  n  (t),  has  a  finite  limit  and 
physical  meaning.  The  ratio  of  the  probability  density  functionals 

ic 

P(^  1  indicates  by  how  much  the  realization  n*(t)  is  more  or  less  probable 
**  ■' 
p(n  ) 

than  the  realization  n  (t) .  Thus,  the  probability  density  functional  of 
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white  noise  is  defined  to  within  the  coefficient  k.  This  does  not 
result  in  misunderstandings  when  using  the  concept  of  "white  noise", 
since  radar  observation  theory  does  not  actually  examine  the  interference 
probability  density  functionals,  but  rather  their  ratios. 

Chapter  9  investigates  the  case  in  which  the  intensity  of  the  interference 
Nq  is  an  unknown  random  quantity  distributed  as  p^(Nq).  In  this  case 

expression  (1.79)  for  the  interference  probability  density  functional 
is  transformed  to  the  following: 


m  00 


'  t 


(1.80) 


1.3.  Combined  Transmission  of  Signal  and  Interference  Through 
Receiver 

We  shall  assume  that  the  receiver  consists  of  a  linear  radio  frequency 
amplifier  (incorporating  in  the  general  case  an  HF  amplifier,  a  converter 
and  an  i.f.  amplifier),  an  inertialess  detector  and  a  linear  video 
amplifier,  tuned  to  the  frequency  of  the  received  signal.  Input  to 
the  receiver  is  the  oscillation 


.r(0  =  n(/)-f 


(1.81) 


consisting  of  random  interference,  which  can  be  approximated  by  white 
noise  with  spectral  intensity  Nq/2,  and  the  valid  signal 


—  =  —  t.)  cos  {t  —  x.)  4-  ?]. 


(1.82) 


We  recall  that  the  real  spectral  intensity  of  the  interference  at  the 
input  of  a  receiver  is  N^.  We  are  examining  the  transmission  of 

amplitude-modulated  (or  pulsed)  signals.  However,  as  will  be  shown  later. 
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the  findings  can  also  be  used  for  other  types  of  radar  signals. 


We  shall  designate  the  frequency  responses  of  the  radio  frequency  amplifier 
and  video  amplifier  a(f)  and  A(f),  and  their  Fourier  transforms,  i.e.,  the 
impulse  responses  of  amplifiers,  h(t)  and  H(t).  The  frequency  response  of 
the  radio  frequency  amplifier  can  be  represented  as 

=  — /o)-f  (1.83) 


where  is  a  function  which  is  practically  non-zero  in  a  comparatively 

(with  fp)  narrow  frequency  band  near  f=0,  such  that 


Furthermore,  introducing  the  notation 

Aa  (/)->«„(/),  (1.84) 


we  find 


A(0=  J  a(/)e'''^'rf/  =  2/i„(0cos2Ky:.  (1.85) 

;  —00 


The  function  2hj^(t),  determined  by  means  of  (1.84),  is  the  envelope  of 
the  impulse  response  of  the  radio  frequency  amplifier. 

We  shall  give  the  respective  subscripts  1,  2  and  3  to  the  different 
voltages  at  the  output  of  the  radio  frequency  amplifier,  the  detector 
and  the  video  amplifier. 

The  voltage  Xj^(t)  at  the  output  of  the  radio  frequency  amplifier  is  equal 
to  the  sum  of  the  signal  voltage 


19 


00 

— ''o)=  J  j(/  — tfl— = 

—00  r 

=  cos  [2j:/|,  (/  —  tj)-}-?]  *5  (2) /in  (/ ^  Xj  —  z)dz 

—  OP 


and  the  envelope 


-^1  (if-'T- to)  =  J  ^{2)ha{t  —  \  —  i)dz 


(1.86) 


(1.87) 


[PAGES  24  AND  25  OF  ORIGINAL  TEXT  OMITTED.] 


2.  Square-law  detection  of  a  strong  signal 


X,  (t)  x,)f2S,  {t  -  X,)  (0. 


(1.99) 


3.  A  weak  signal 


I 


(1.100) 


The  latter  expression  reflects  accurately  enough  the  transmission  of  a 
weak  signal  through  a  square-law  detector;  it  is  extremely  approximate 
for  a  linear  detector.  The  accuracy  of  formula  (1.100)  for  linear 
detection,  however,  is  not  important  for  us,  since  we  are  employing  a 
square-law  detector  for  weak  signals  (Chapters  6  and  8) . 

The  case  of  linear  detection  of  a  strong  signal  is  fundamental  in  the 
investigations  below.  We  note  that  the  output  of  a  linear  strong-signal 
detector  is  the  same  as  the  output  of  a  synchronous  detector.  In  both 
cases  the  phased  component  of  the  interference  N^j^(t)  is  retained  at 

the  output  along  with  the  signal  envelope.  The  quadrature  component 

N  ,  (t)  is  almost  fully  suppressed.  With  the  matched  frequency  response 
si 

(1.90)  the  valid  component  S^(t-TQ)  and  the  correlation  function  of  the 

interference  N  , (t)  at  the  output  of  a  linear  strong-signal  detector 
cl 

have  the  same  form  as  the  modulation  correlation  function 


(1.101) 

(1.102) 


Finally,  the  voltage  X^(t)  at  the  video  amplifier  output  will  be 


X,(0  =  J  X,(z)II{t-'z)dz.  ■  (1.103) 

—00 


21 


For  linear  detection  of  strong  signals,  the  latter  expression  is 
transformed  to 


;r,  (/)  =  5 15  (J  - +  N,  Wl 


(1.104) 


where  H^(t)  is  the  generalized  impulse  response  of  the  receiver 

Zj 


(1.105) 


which  is  the  Fourier  transform  of  the  generalized  frequency  response  A^(f) 


(1.106) 


IS 


Thus  in  the  case  of  linear  detection  of  strong  signals  the  video  amplifier 
output  voltage  X3(t)  is  fully  determined  by  the  generalized  response 

the  receiver  H^(t)  for  A^(f),  and  it  is  irrelevant  how  this  response  - 

distributed  among  individual  components  of  the  receiver,  especially 
between  the  i.f.  amplifier  and  video  amplifier. 

Analogously,  for  square  law  detection  of  strong  signals 

=  f  5  [S(2l■T-x,)S(^,-'«o)  +  2S(^^  — t,)N^(2.)l  X 

^  XH^{t  —  z^,t  —  z,)dz^d^t 

and  for  detection  of  weak  signals 

x.(0  =  fj  +  (1.108) 

.  ■  +  W.T*.)  N.  .)l 


(1.107) 
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where  —  generalized  two-dimensional  impulse  response  of 

receiver. 


00  .  y 

=  \  K  (/•  —  K  (^  —  (2)  dz,  (1 . 109) 


which  is  the  Fourier  transform  of  the  generalized  two  dimensional 
frequency  response  A^Cf^^,  f2) 


\(h  /.)  /I (/,+/.). 


(1.110) 


i.e. , 


(1.111) 


Thus,  in  the  case  of  square  law  detection  of  strong  and  weak  signals 
the  video  amplifier  output  voltage  X^Ct)  is  determined  by  the  generalized 

two-dimensional  receiver  response  H2(tj^,t2)  or  A^(f^,  12)  and  it  is 

irrelevant  how  this  response  is  distributed  among  the  individual 
components  of  the  receiver. 

The  case  is  also  encountered,  in  Chapters  5  and  7,  in  which  the 
envelope  detector  has  as  input  the  oscillation 


.¥,  (/)  =  1 S,  (/  —  x^)  1  co,s  [2k/„  (/  —  tj  -f  X  (^  —  xJ  -j-  y]  4-  rz,  (/), 


(1.112) 


consist.ing  of  random  interference  and  a  signal  with  simultaneous 
amplitude  and  phase  modulation.  If  we  use  the  representation  (1.59) 
for  the  random  interference  n^(t),  the  output  of  the  strong-signal 

detector  ca'n  be  represented  as 

A\(/)=^15,(/-x.)l  +  iV„(0.  ■  .  (1. 

The  low  frequency  random  process  N  (t)  is  defined  by  formula  (1.60). 

CK. 
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CHAPTER  2.  Statistical  Treatment  of  Radar  Observation  Process 
2.1.  Statistical  Model  of  Radar  Observation  Process 

The  observation  process  performed  by  a  radar  operator  consists  of  detecting 
signals  reflected  from  targets,  and  estimating  certain  parameters  of  these 
signals,  such  as  the  parameters  which  determine  the  target  coordinates. 

The  interference  which  unavoidably  occurs  in  any  radio  link  masks  the 
valid  signals  and  makes  it  impossible  to  use  the  received  oscillation 
to  establish  with  absolute  certaintly  that  a  valid  signal  is  present  and 
to  estimate  its  parameters  precisely.  The  occurrence  of  interference 
causes  the  observation  results  to  be  random,  and  makes  it  necessary  to 
study  these  results  by  statistical  methods. 

Under  actual  operating  conditions  the  decisions  made  by  an  operator 
when  detecting  signals  and  in  estimating  signal  parameters  are  based  on 
certain  subjective  criteria  which  depend  upon  the  level  of  training  and 
other  properties  of  the  operator.  Subjective  operator  properties  have 
some  influence  on  radar  observation  performance  indicators.  However,  the 
limiting  capabilities  and  indicators  of  radar  systems  are  determined 
primarily  by  the  statistical  nature  of  the  problem.  A  statistical 
statement  of  the  problem  requires  only  that  we  free  the  radar  observation 
process  of  subjective  "baggage",  assuming  that  the  operator  acts  on  the 
basis  of  clearly  formulated  decision  making  rules,  and  replacing  the 
operator  with  a  decision  device.  Replacing  the  operator  with  a  machine  — 
a  computing  and  decision  device  —  also  corresponds  to  the  current  trend 
in  the  development  of  radar  systems. 

We  shall  represent  the  statistical  model  of  the  radar  observation  process 
to  be  examined  below  in  the  form  of  the  diagram  shown  in  Fig.  2.1.  The 
main  components  of  this  diagram  are  the  signal  space  E,  the  interference 
space  N,  the  received  oscillation  space  X,  the' receiver ,  the  output  space  Y, 
the  decision  device  and  the  decision  space  S*. 

The  signal  space  E  consists  of  the  set  of  all  possible  valid  reflected  _ 
signals  arriving  at  the  input  of  the  receiver.  The  set  of  valid  sig/nals 
can  be  defined  in  principle  such  that  it  includes  the  case  of  receiving 
several  reflected  signals  simultaneously.  However,  we  shall  assxime  that 
it  is  possible  for  only  one  single  target  to  be  present  in  the  space 
which  the  radar  is  examining.  A  brief  discussion  of  the  possibility  of 
the  simultaneous  presence  of  a  large  number  of  targets  (resolution  problem) 
will  be  presented  later. 
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Fig.  2.1.  Statistical  model  of  radar  observation  process 


The  reflected  signal  is  formed  as  the  result  of  interaction  between  the 
radiated  signal,  the  reflecting  object  and  the  environment.  The 
reflected  signal  can  be  represented  analytically  as  some  known  function 
(determined  by  the  nature  of  the  radiated  signal)  of  time  and  a  series  of 
random  parameters.  The  random  parameters  in  the  reflected  signal  to  be 
estimated  will  be  called  useful  parameters.  The  random  parameters  of  the 
reflected  signal  which  are  not  to  be  estimated  and  which  are  statistically 
independent  with  respect  to  the  useful  parameters,  i.e.,  not  containing 
information  about  the  useful  parameters,  will  be  called  parasitic 
parameters. 

The  reflected  signal  can  thus  be  written  as  the  function  3(t,s,u),  where 
s  and  u  are  random  vector  quantities  which  represent  respectively  the  set 
of  useful  and  parasitic  parameters.  We  note  that  the  reflected  signal  can 
consist  of  the  sum  of  the  elementary  signals  3=3^+. . .+3^) ,  differing  from 

one  another  in  carrier  frequency,  modulation  pfinciple  or  by  the  fact  that 
they  represent  the  state  of  the  electromagnetic  field  at  different  points 
in  real  space.  It  is  assumed  that  the  reflected  signal  3(t,s,u)  has  a 
finite  fixed  duration  of  T  sec.  The  signal  space,  besides  the  set  of 
functions  3(t,s,u),  also  includes  the  null  signal  which  denotes  the 
absence  of  a  target.  We  shall  designate  the  event  "target  absent"  by 
the  symbol  Sq  and  consider  that  Sq  is  one  of  the  possible  values  of  the 

random  vector  quantity  s. 

Each  point  of  the  signal  space  E  is  thus  uniquely  determined  by  the  values 
taken  on  by  the  parameters  s  and  u.  It  will  be  helpful  below  to 
introduce  as  well  the  concept  of  useful  parameter  space  S,  or  the  set  of 
all  possible  values  of  the  random  parameter  s. 
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As  a  rule,  random  parameters  can  take  on  arbitrary  values  within 
certain  intervals.  However,  in  most  theoretical  and  experimental 
studies  it  is  useful  to  idealize  the  signal  space  somewhat  in  order 
to  determine  the  basic  dependencies  and  regularities.  It  can  be 
assumed,  for  example,  that  some  of  the  parasitic  parameters,  or 
even  all  of  the  parasitic  parameters,  in  a  reflected  signal  are  non- 
random  given  quantities.  This  idealization  makes  it  possible  to 
determine  separately  the  influence  of  each  of  the  parasitic 
parameters  on  the  radar  observation  process. 

Analogously,  it  can  be  assumed  that  only  a  single  useful  parameter 
is  to  be  estimated,  such  as  the  delay  time,  and  that  the  remaining 
useful  parameters  are  fixed,  i.e.,  given.  In  this  case  the  scanning 
is  done  with  respect  to  only  one  parameter  or  one  coordinate,  such  as 
the  range  coordinate  in  some  fixed  direction.  The  latter  assumption 
makes  it  possible  to  resolve  a  number  of  important  problems  in  radar 
theory. 

Finally,  an  idealization  is  employed  in  which  the  set  of  useful 

parameters  s  can  take  on  only  a  finite  number  of  fixed  values  s^, 

s„...,  s  ,  as  well  as  the  value  s..  Radar  lines  with  a  discrete 
2  m  U 

set  of  useful  parameters  will  be  called  simplified  radar  lines. 

The  present  chapter  imposes  no  limitations  on  the  class  of 
interference  contained  in  the  space  N.  It  is  assumed  only  that 
the  distribution  p(n)  is  defined  for  the  random  vector  quantity  n 
representing  the  realization  of  the  interference  from  the  class  in 
question.  Added  together,  the  interference  and  the  signal  form  the 
received  oscillation  x(t),  which  is  represented  by  the  point  x  in 
the  space  X.  The  system  must  provide  a  unique  transformation  of  the 
received  oscillation  space  X  to  the  decision  space  S*,  which  consists 
of  the  set  of  all  possible  decisions  s*  made  during  the  radar  observation 
process.  The  decision  space  has  the  same  structure  as  the  useful 
parameter  space  3. 

Using  as  guidance  the  consideration  of  convenience  in  practical 
construction  of  the  system,  the  latter  is  assumed  to  consist  of  two 
parts:  a  receiving  device  and  a  decision  device.  Accordingly,  the 
operation  of  transforming  the  received  oscillation  into  a  decision  is 
split  into  two  stages.  The  first  stage  involves  transforming  the  received 
oscillation  into  the  output  of  the  receiver;  the  second  stage  consists  of 
transforming  the  output  into  a  decision. 
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2.2.  Statement  of  Problem  of  Optimal  System  Definition 

An  optimal  system  is  one  which  provides  the  decisions  or  estimates 
which  are  best  in  terms  of  a  selected  optimality  criterion.  The 
optimality  criterion  is  selected  so  that  it  corresponds  sufficiently 
well  to  the  nature  of  the  problem  being  solved.  This  rule  does  not 
define  the  criterion  in  only  a  single  way,  and  leaves  open  the 
possibility  of  substantial  arbitrariness  in  selecting  it.  The  extent 
to  which  the  selected  criterion  can  serve  as  the  basis  for  analytical 
solution  of  the  problem  should  also  be  kept  in  mind. 

An  estimate  is  a  random  quantity,  and  its  quality  can  be  determined 
only  over  a  large  number  of  trials.  For  this  reason,  practically 
all  optimality  criteria  can  be  reduced  to  ensuring  the  minimum  or 
maximum  mathematical  expectation  of  some  function  of  the  signal 
and  decision  (estimate).  The  general  optimality  criterion  for  the 
radar  observation  process  can  thus  be  written  as  [28] 


p  J  r  (s*,'  s)  p  (s)  p  (s'/s)  ds*t/s  =  extremum, 


$•  s 


(2.1) 


where  s*  J.s  the  estimate  of  the  random  quantity  s,  i.e.,  the  decision 
made  as  t‘he  result  of  making  the  radar  observation.  The  set  of  all 
possible  decisions  S*,  which  includes  integration  with  respect  to  s*, 
consists  of  the  interval  S*  corresponding  to  the  working  range  of 
variation  of  the  useful  signal  parameters,  and  of  the  point  s*  ("target 

absent"  decision)  outside  the  interval  S*.  For  convenience  in  writing 
criterion  (2.1)  the  decision  set  S*  is  augmented  by  the  vicinity  of 
the  point  s*.  The  useful  parameter  set  S  has  the  same  structure  as  the 

set  §*,  i.e.,  it  includes  the  working  interval  S  of  variation  of  all  of 
the  useful  parameters  s  and  the  point  s^  together  with  its  vicinity. 

r(s*,s)  is  the  value  function,  which  determines  the  relative  value  or 
significance  of  the  combined  realization  of  a  signal  with  parameters  s 
and  the  decision  s*.  In  all  cases  here  and  below  in  which  no  misunder¬ 
standing  can  arise,  the  probability  densities  and  probabilities  of 
different  quantities  are  denoted  by  the  same  letters:  p  —  probability 
density,  P  —  probability.  Conditional  probabilities  are  designated  by 
the  same  symbols  as  unconditional  ones,  and  differ  only  by  the  presence 
of  a  slash  in  the  argument,  to  the  right  of  which  is  indicated  the 
conditional  quantity  or  event.  For  example,  p(s*/s)  —  conditional 
probability  density  for  the  quantity  s*  given  realization  of  the 
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quantity  s  (signal  with  useful  parameters  s) .  The  character  ~  over  the 
symbol  p  means  that  the  probability  density  function  includes  the  discrete 
point  s*  or  Sq.  Thus 


1p(5.)5(s-s.), 


ses. 

sGS. 

s*6S*. 


^  [P(sl/s)H^^sl),  s-GS*. 


(2.2) 


The  symbol  G  indicates  "membership",  while  the  symbol  G  indicates  non¬ 
membership.  The  function  p(s)  determines  the  a  priori  probability  of 
different  outcomes. 

Simplified  lines  can  be  viewed  as  a  special  case  of  real  lines,  for 
which  the  set  of  useful  parameters  s  can  take  on  only  one  of  the 

m+1  fixed  values  s-,  s,  ...,  s  .  When  using  simplified  lines  it  should 

U  1  m 

be  assumed  that 


^(s)  =  5]  P(s,)8(s  — s,), 

1=1 

/'(s7s)=£p(s;/s)5(s*-s;). 

1=1 


(2.3) 


Accordingly,  the  general  optimality  criterion  (2.1)  for  simplified 
lines  becomes 


mm 

p  =  ^  r  (sj ,  s,)  P  (s,)  P  (s*  /Sf)  =  extremunt 
<_oy=o 


The  problem  of  defining  an  optimal  system  consists  of  determining  the 
correspondence  (called  the  system  operator)  between  the  received  oscillation 
X  and  the  decision  s* 
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s'  =  A  (x). 


(2.5) 


which  would  satisfy  the  optimality  condition  —  condition  (2.1)  in  the 
general  case.  Determination  of  the  optimal  operator  actually  means 
determining  optimal  processing  of  the  received  oscillations  in  the 
receiver  and  an  optimal  decision  rule  in  the  decision  device.  A  practically 
important  extension  of  the  problem  of  defining  an  optimal  system  is  the 
calculation  of  the  performance  indicators  of  the  system.  The  performance 
indicators  of  an  optimal  system  are  called  the  limiting  or  potential 
Indicators,  since  they  characterize  the  theoretically  limiting  capabilities 
of  the  system. 

The  solution  of  the  problem  can  be  found  in  general  form  for  real  lines. 
Then  the  decision  for  simplified  lines  is  obtained  by  applying  formulas 
(2.3).  However,  it  is  more  convenient  to  begin  our  examination  with 
simplified  lines,  since  this  makes  it  possible  to  retain  the  connection 
with  the  customary  representations  which  have  been  developed  primarily 
for  the  two -alternative  signal-no  signal  situation  [6,  15,  61,  63,  etc.]. 

Criteria  (2.1)  and  (2.4)  are  Bayesian  if  the  function  r(s*,s)  is  assigned 
in  advance  and  its  form  is  independent  of  the  decision  making  rule.  Of 
all  known  practical  criteria,  only  the  criterion  of  minimum  information 
loss,  for  which 


r(s*,  s)  =  — log/7(s/s*), 


does  not  belong  to  the  class  of  Bayesian  criteria  [28].  In  investigating 
the  radar  observation  process,  the  minimum  information  loss  criterion 
provides  essentially  no  advantage,  but  it  makes  problem  solving  much 
harder.  Therefore,  it  is  assumed  below  that  the  function  r(s*,s),is 
independent  of  the  manner  in  which  the  estimate  s*  is  obtained. 


2.3.  Simplified  Lines 

The  oscillation  x  is  input  to  the  receiver.  After  appropriate  processing 
of  this  oscillation  it  is  necessary  to  make  some  decision  s*  (j=0,l, . . . ,m)  . 

In  the  final  analysis,  without  separating  the  functions  of  the  receiving 
and  decision  devices,  the  decision  is  based  on  the  fact  that  the  received 
oscillation  space  X  is  divided  into  nri-1  regions;  Xq,  X^^...,  X^. 
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Depending  upon  the  region  into  which  falls  the  point  x,  representing  the 
received  oscillation  x(t),  a  corresponding  decision  s*  is  made.  This  means 

that  correspondence  (2.5)  takes  on  the  appearance 


/  s’  =  s*  for  (1  =  0, 


(2.6) 


and  it  is  necessary  only  to  define  the  algorithm  for  dividing  the  received 
oscillations  space  X  into  different  decision  regions. 

The  probability  P(s*/s^)  that  for  the  signal  s^  arriving  at  the  input  of 

the  receiving  device  together  with  interference  the  decision  s*  will  be 

made,  which  is  the  same  as  the  probability  that  the  oscillation  x  will 
fall  in  the  region  X.  when  a  signal  with  parameters  s^  is  present, 


^(Sy  /s,)={  p(x/s,)dx. 


(2.7) 


The  conditional  distribution  function  p(x/s)  of  the  quantity  x,  viewed 
as  a  function  of  s  with  a  fixed  value  of  x,  is  called  the  likelihood 
function.  For  additive  interference  the  analytical  expression  of  the 
function  p(x/s)  is  determined  directly  by  the  given  noise  distribution 
p(n) . 

The  quantities  in  (2.7)  with  the  same  subscripts  for  the  signals  and 
decisions  define  the  conditional  correct  decision  probabilities.  The 
quantities  in  (2.7)  with  different  signal  and  decision  subscripts  define 
the  conditional  probabilities  of  various  types  of  decision  making  errors. 

An  error  in  which  the  decision  s^  (i=l,2,...,  m)  is  made  that  one  of  the 

non-null  signals  is  present  when  there  is  no  signal  at  the  input,  is 
called  a  false  alarm.  The  false  alarm  probability  F  is 

Fl=  (s;  Is.)  =  1  -  P  (s;  |s;).  (2.8) 

»=l  - 
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An  error  in  which  the  decision  s^*  —  "no  target"  is  made  when  one  of 
the  non-null  signals  is  present  at  the  input,  is  called  a  missed  signal. 

We  shall  assvime  that  the  value  function  r  is  defined  such  that  the  equation 
p  =  nay  is  the  optimality  criterion.  Before  defining  the  algorithm  in  which 
we  are  interested  for  dividing  the  received  oscillation  space,  we  shall 
formulate  the  following  lemma. 

If  (j)(x,0),  4)(x,l),...,  (|)(x,m)  —  is  a  set  of  functions  defined  in  the  space 
X  such  that  the  equality  (})(x,i)=(J)(x, j)  for  i?fj  can  occur  only  on  the  set  of 
points  of  measure  zero,  and  Xq,  X^^ . X^  are  non- intersecting  regions 

filling  the  space  *  then 

'i=o 


^  J  tp (X,  <) rfx  =  max  (2.9) 

Xl 


only  when  the  regions  Xj^(i=0,l, . . . ,  m)  are  defined  through  the  system  of 
inequalities 


?  (X,  0  >  f  (x,  /),  (/  =  0,  1 , . . m,  /  /).  (2 . 10) 


Let  us  compose  the  function 


<D(x)  =  sup[y(x,  0),  <p(x,  l),...-,9(x,  m)], 


where  the  symbol  sup  (supremum)  denotes  the  largest  of  the  values 
contained  in  the  square  brackets.  Then 

(®(x)rfx>^  f  y(x,  i)rfx.  (2.12) 
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The  equals  sign  occurs  in  the  latter  expression  only  if  for  all  1 


J  ®(x)rfx==  f  f  (x,  /)  rfx,  (2.13) 


which  leads  directly  to  (2.10). 

The  stipulated  condition  that  the  measure  of  the  set  of  points  at  which 
(})(x,i)=4)(x,  j)  is  zero  can  be  omitted.  Then  the  and  >  signs  in  system 
of  Inequalities  (2.10)  must  be  combined.  For  example,  for  all  j>i  we 
write  the  sign  >■,  and  for  all  j  <  1  we  write  the  sign  >.  Keeping  this 
remark  in  mind,  we  can  later  not  concern  ourselves  with  weather  or  not 
the  aforementioned  condition  is  satisfied. 

We  now  transform  expression  (2.4)  as  follows: 


P  ==  S'S  '■  ^  ^ (2 . 14) 

/=0  Xj  <=0 


Based  on  the  proof  above,  we  conclude  that  condition  (2.14)  of  the  maximum 
of  the  quantity  p  is  satisfied  if  the  system  of  inequalities 

S  r  is* ,  s,)  P  (s)  p  (x  1  s^)  >  5]  r  (s’ ,  s,)  P  (s,)  p  (x  j  s,) 

•.  ^  '  (2.15) 

(A  0,  I  ,  .  .  ttlf  k  /).  ' 

is  used  as  the  definition  of  the  region  (j=0,l,...,  m)  . 

The  problem  is  solved  completely  analogously  when  the  value  function  r  is 
defined  such  that  the  optimality  criterion  consists  of  minimizing  the 
quantity  p.  In  this  case  it  is  necessary  only  to  introduce  the  function 
inf  (inf inum) ,  which  denotes  the  smallest  value,  instead  of  the  function 
sup.  This  makes  it  necessary  to  replace  the  ^  sign  in  the  final  result 
of  (2.15)  with  ^ . 

This  separation  of  the  regions  makes  it  possible  to  define,  in  general  form 
for  simplified  lines  (l.e.,  based  on  a  general  optimality  criterion). 
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optimal  processing  of  the  received  oscillations  in  the  receiver  and  an 
optimal  decision  making  rule  for  the  decision  device. 

The  boundary  between  the  receiving  and  decision  devices  is  somewhat 
arbitrary,  and  is  determined  by  considerations  of  convenience  in 
practical  circuit  implementation.  It  is  possible  to  have  an  arbitrary 
unique  variation  in  the  output  of  the  receiver  with  a  corresponding 
variation  in  the  decision  making  rule.  It  is  also  possible  for  various 
operations  on  the  received  oscillations  to  be  related  to  either  the 
receiving  or  the  decision  device.  It  is  only  necessary  that  the  joint 
action  of  the  receiving  and  decision  devices  ensure  that  system  of 
inequalities  (2.15)  is  satisfied  when  decisions  are  made.  For  example, 
when  the  oscillation  x  is  received,  the  receiver  might  output  the  sequence 
of  the  quantities  P(s^)  p(x|s^)  for  all  values  of  i  (the  Woodward-Davis 
system) .  In  this  case  the  decision  device  multiplies  the  output  by  the 
corresponding  coefficients  of  r(s*,s),  performs  the  summation 


t-o 


(2.16) 


and  compares  the  values  of  the  sums  with  one  another.  The  decision 

s*  is  made  that  that  signal  s,  is  present  for  which  the  quantity  E  is 

greatest.  The  final  result  remains  the  same  if  the  operations  of 
multiplication  by  s^^)  and  summation  are  done  by  the  receiving 

device  or,  conversely,  if  the  receiver  is  relieved  of  multiplying  by 
P(s.)  by  transferring  that  operation  to  the  decision  device. 

In  theoretical  investigations  it  is  convenient  to  use  the  concept  of 
"likelihood  coefficient"  [63],  which  we  shall  designate  by  the  symbol  A 
and  which  is  defined  as 


J(s)  =  -4iM..  ‘2 -17) 


Like  the  likelihood  function  in  the  numerator  of  (2.17),  the  likelihood 
coefficient  is  a  function  of  the  signal  parameters  (or,  in  the  general 
case,  a  function  of  the  point  in  the  signal  space)  for  a  fixed  value  of 
the  received  oscillation  x. 
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Keeping  (2.17)  in  mind,  the  following  appearance  can  be  given  the 
system  of  inequalities  which  defines  the  division  of  the  received 
oscillation  space. 


Region  of  null  signal 


m 


2  [/•  (S* ,  s)  —  r  (s' .  s^)l  P  (s)  A  (s,)  + 

/=i 

+  P  (s.)  r  (s*  ,  Sa)<P  (s,)  r  (s* ,  s,)  (/=  1 ,  2, 


(2.18) 


Region  of  arbitrary  non-null  signal  X 


2  ['■  (sj ,  s,)  —  r  (sj ,  s,)]  P  (s^)  A  (s,)  + 

i=i 

-j- P (So) /•  (Sy  ,  %)>P (s,) r (s* ,  s,), 


[r  (s* ,  s,)  —  r  (s* ,  s)]P  (s,)  A  (s,)  +  P  (s,)  r  (s) ,  s^)  ^ , 

’/=r 

S  ['■  (si  •  Sy)  —  /•  (si ,  s,)l  P  (s^)  A  (s^)  +  P  (s,)/  (si ,  sj 

<=i  /• 

(&=1,  2 . .  k^l).  • 


(2.19) 


(2.20) 


The  optimal  decision  algorithm  obtained,  based  on  the  general  optimality 
criterion,  is  fully  applicable  if  for  all  real  operating  conditions  the 
value  ftmctions  r(s*,  s^)  and  the  a  priori  distribution  P  (s*)  are  defined. 

However,  in  most  cases  this  definition  is  lacking.  First  of  all,  the 
so-called  "a  priori  difficulty"  [8]  occurs,  which  consists  of  the  absence 
of  enough  a  priori  data  to  define  the  function  p(s).  Second,  statistical 
theory  assumes  the  value  function  r(s*,s)  to  be  given,  and  does  not  deal 
with  selecting  it.  The  value  function  r(s*,s)  can  in  principle  be 
defined  (albeit  not  uniquely)  based  on  the  specific  tactical  conditions 
under  which  the  system  is  used.  In  this  case  the  function  r  will  be 
different  for  different  systems  and  for  different  operating  conditions 
of  the  same  system.  On  the  other  hand,  as  the  function  r  changes,  the 
meaning  and  results  of  the  decision  may  change  significantly.  If  the 
objective  is  to  reduce  theory  to  analytical  formulas  for  system  performance 
indicators,  it  is  necessary  to  investigate  systems  with  a  specifically  assigned 
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value  function,  instead  of  operations  with  the  function  r  in  general 
form. 

Different  optimality  criteria  which  are  encountered  in  solving 
practical  problems  are  usually  special  cases  of  criterion  (2.4)  and 
(2.1)  with  various  assumptions  regarding  the  functions  p(s)  and  r(s*,s). 

We  note  that  the  minimax  method  employed  in  game  theory  [1,55,62]  does 
not  allow  the  "a  priori  difficulty"  in  radar  to  be  overcome  satisfactorily. 
According  to  this  method,  the  observation  process  is  viewed  as  a  game 
between  an  observer  and  nattire.  The  value  function  provides  a  quantitative 
measure  of  the  result  of  each  realization  of  the  game.  It  is  assumed  that 
nature  selects  a  strategy  which  is  least  favorable  for  the  observer,  i.e., 
the  function  p(s).  With  this  assumption,  based  on  Bayesian  criterion 
(2.1)  or  (2.4),  the  observer  chooses  his  strategy  (decision  making  rule). 
Under  actual  operating  conditions,  the  a  priori  distribution  may  differ 
significantly  from  the  least  favorable  distribution.  Accordingly,  the 
decision  made  by  the  observer  will  not  be  the  best  ones.  In  statistical 
radar  theory,  the  mlnimax  method  should  be  viewed  as  one  possible,  rather 
arbitrary  method  for  assigning  the  lacking  a  priori  data  required  in 
order  to  use  the  Bayesian  criterion. 

In  accordance  with  existing  engineering  practice,  we  limit  ourselves  in 
practical  applications  to  the  Neumann-Pearson  criteria  and  an  ideal 
observer. 

The  ideal  observer  criterion  [18,  26,  62],  which  provides  the  maximum 
average  probability  of  making  a  correct  decision 


m 


1-0 


max. 


(2.21) 


is  obtained  from  general  optimality  criterion  (2.4),  if  we  assume* 


r(s],  s,)  =  V 


(2.22) 


where  6 . . 


Kronecker's  symbol 


for  /  =  /, 
,0  for  /v^/. 


(2.23) 
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After  assumption  (2.22)  is  made,  there  still  remains  a  fundamental 
difficulty  which  is  associated  with  the  ambiguity  of  the  a  priori 
distribution.  The  simplest  way  out  of  this  difficulty  is  to  assign 
an  arbitrary,  but  sufficiently  likely  a  priori  distribution,  such 
as  a  uniform  distribution 


(2.24) 


where  ]i  —  probability  of  presence  of  non-zero  signal,  which  must  also 
be  assigned.  Then  the  systems  of  inequalities  (2.18)  and  (2.19),  which 
determine  the  separation  of  the  space  X  into  different  decision  regions, 
take  on  these  forms: 

—  region  of  null  signal 

A(s,)<^^m(/=1,  2,...,  m),  .  (2.25) 


—  region  of  arbitrary  non-null  signal  X^(i=l,2, . . . ,m) 

A  (s,)  ^  A  (s)  (/=1 ,  2, . . . ,  m;  /  ^  /):, 


(2.26) 


The  decision  rule  based  on  the  Neumann-Pearson  criterion  leads  to  a 
similar  result.  The  Neumann-Pearson  criterion -is  usually  defined 
for  the  two-altemative  situation.  We  shall  introduce  an  analogous 
definition  for  the  case  of  an  arbitrary  number  of  alternatives. 


We  shall  assume  that  according  to  the  Neumann-Pearson  criterion  it  is 
necessary  to  maximize  the  average  conditional  probability  of  making  a 
correct  decision  with  respect  to  non-null  signals  .  “  , 


<=i 


assuming  that  the  false  alarm  probability  is  equal  to  some  given  quantity 
F.  The  formulated  criterion  can  be  written  analytically  as  follows: 
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(2.27) 


P=A.p(s;isj-i-jP(s;is,)=  • 

<=i 

m 

=1  A,/7(;c/5jc/x-f  JjJ/7(xls,)rfx=  nax 
X,  l=i  Xi 


assuming  that 


f  p(xl5jrfx=l  — f,  (2.28) 

X, 


where  Ag  —  an  arbitrary  constant  coefficient. 

According  to  the  above,  of  all  of  the  possible  algorithms  for  dividing 
the  received  oscillation  space  X  into  regions  Xq,  X^^,...,  X^,  the 

greatest  value  of  the  sum  in  (2.27)  occurs  when  the  definition  of  the 
arbitrary  region  X,  of  a  non-null  decision  s*  consists  of  the  system  of 
inequalities 

A(s,)>A„ 

A  (Si)> A  (s^)  (/=!,  2,...',  m,  i^i),  (2.29) 


and  the  region  X^  of  the  null  decision  s*  is  determined  through  the 
system 

'A(Sj)^Aj  (/=l,  2,  /7l).  (2.30) 


As  Aq  varies  from  0  to  “  the  region  X^  varies  from  0  to  X,  the  false 

alarm  probability  accordingly  varies  from  1  to  0.  We  select  the  quantity 
Aq  such  that  (2.28)  occurs,  i.e.,  such  that  the  probability  that  system 

(2.30)  will  be  satisfied  when  the  null  signal  s_  is  input  is  exactly  equal  to 
the  assigned  quantity  1-F.  The  coefficient  A^  thus  defined  will  be  called 
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the  critical  (or  threshold)  level  of  the  likelihood  coefficient. 


Systems  of  inequalities  (2.29)  and  (2.30),  assuming  that  the  coefficient 
Aq  figuring  in  them  is  the  critical  level  of  the  likelihood  coefficient, 

facilitate  simultaneous  satisfaction  of  equations  (2.27)  and  (2.28),  and 
provide  the  sought  separation  of  received  oscillation  space  X  into  regions 
Xq,  S^,...,  X^.  Dividing  the  space  X  into  any  other  regions  X^,...,  X^ 

makes  the  sum  of  (2.27)  smaller,  and  consequently  reduces  the  quantity 


m 


since  the  first  term  AqP(s*|sq) 


must  remain  unchanged. 


The  Neumann-Pearson  criterion  cited  above  is  obtained  from  the  general 
optimality  criterion  (2.4)  if  we  introduce  the  following  initial  assumptions 
condition  (2.22)  and  a  uniform  distribution  of  a  priori  probabilities 
(2.24)  with  a  probability  of  occurrence  of  a  non-null  signal  of 


1. 

m 

/n+A, 


■ 


(2.31) 


Thus,  from  the  practical  viewpoint  it  is  completely  irrelevant  whether 
assumptions  (2.24),  which  are  arbitrary  to  some  extent,  are  assigned  with 
respect  to  the  a  priori  probabilities  when  the  ideal  observer  criterion  is 
used,  or  whether  the  Neumann-Pearson  criterion  is  used,  which  is  also 
based  on  an  arbitrary  false  alarm  probability.  In  both  cases  the  division 
of  the  received  oscillations  space  is  determined  by  systems  of  inequalities 
of  the  form  (2.20)  and  (2.30).  We  note  that  when  there  are  specific  data 
concerning  the  a  priori  distribution  and  substantiated  considerations 
regarding  the  choice  of  value  function,  allowance  must  be  made  for  them  in 
defining  the  decision  making  rule. 


2.4.  Real  Lines 

We  shall  modify  general  optimality  criterion  (2.1)  somewhat.  We  first 
establish  the  following  correspondence: 


j/(s*.  s)p(s7s)rfs*  =  s)/7(x/s)</x, 


(2.32) 


where  f(s*,s)  is  an  arbitrary  function  of  s*  and  s, 


X^  is  the  region  in 
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space  X  which  is  uniquely  mapped  onto  the  region  O*  of  space  S*  by  using 
the  operator  of  system  (2.5).  In  order  to  prove  (2.32)  we  partition  the 
region  O*  into  elementary  regions  with  volume  As*.  We  assume  that  sg  G  a*. 

We  give  each  elementary  region  a  corresponding  ordinal  number  1,  2 . m. 

In  any  elementary  space,  say  the  jth,  the  function  f (s*,s)p(s*/s)  can  be 

assumed  constant  and  equal  to  f (s*,s)p(s*/s) ,  where  s*  is  an  arbitrary 

J  J  J 

fixed  value  of  the  variable  s*  in  the  jth  elementary  region.  Then 


'  (2.33) 

m-*ot  V  /=:i  ;  ) 


Designating  X,  the  region  of  space  X  corresponding  to  decision  s*,  and 
keeping  in  mind  that 


P  (SjJs)  As* 


■■jp(x/s)dx, 


(2.34) 


we  obtain 


, I  j  ^  ^  } 

jn->oa  V  J—l  Xf  / 


(2.35) 


The  latter  is  transformed  directly  into  the  expression  in  the  right  part 
of  correspondence  (2.32). 

Based  on  (2.32)  and  (2.2), 


f /(S-,  s);(sVs)</s-=/(s;.  s)P(s;/s)+f  f(»'.  s);7(s78)rfs-= 
^  $• 
s» 

s)p{xls}dx-{-  J  /(s*,  s)/?(x/s)rfx= 

X.  x-x, 

—  \I{sr,s)p{xls)dx, 


(2.36) 
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Accordingly,  criterion  (2.1)  can  be  represented  in  the  form 


P  —  ^dx^dsr  (s*,  s)  p  (s)  p  (x/s)  =  extremum. 


X  ~ 


(2.37) 


We  shall  assume,  as  above,  that  the  optimality  criterion  consists  of 
maximizing  the  quantity  p.  Then  as  the  oscillations  x  arrive  the 
decision  which  must  be  made  is  s*(s*  6  2*),  for  which 


X 


Jr  (s^  ,  s) p (s)  A  (s) ds  =  sup  J r  (s*,  s)  p  (s)  A  (s) rfs  J,  (2 . 38) 

s  ,  '  f 


where  sup[(j)(x,s*)  ]  designates  the  greatest  value  of  the  function  (J)(x,s*) 
taken  for  a  given  x  with  respect  to  all  possible  s*. 

The  solution  obtained  for  the  problem  of  defining  an  optimal  system 
in  terms  of  a  criterion  of  the  general  form  (2.1)  can  be  written  with 
the  help  of  system  operator  (2.5)  in  a  fashion  analogous  to  that  used 
for  simplified  lines: 


s*  =  Sjj  for  X  0  Xjj, 


(2.39) 


where  X  ‘ —  the  set  x  (defined  for  every  value  of  s*  from  decision  space 

X  X 

2*  on  which  equation  (2.38)  is  satisfied. 

The  attempt  to  use  the  general  solution  of  problem  (2.38)  encounters  the 
same  difficulties  which  were  mentioned  in  the  preceding  section. 
Accordingly,  it  is  again  necessary  to  employ  special  types  of  criteria 
based  on  certain  assumptions.  A  large  number  of  more  or  less  sensible 
optimality  criteria  can  be  proposed  for  t^he  detection  process  in  a  real 
line.  We  shall  limit  ourselves  to  two  criteria,  which  are  continual 
analogs  of  the  Neumann-Pearson  and  ideal  observer  criteria.  The  physical 
meaning  of  the  criteria  is  retained,  and  by  substituting  (2.3)  they  are 
transformed  to  the  corresponding  optimality  criteria  for  a  simplified  line. 
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The  Neumann-Pearson  criterion,  which  equipment  developers  consider  to  be 
most  applicable,  can  be  defined  in  generalized  form  by  the  equation 


p  =  A  (s^/s  J-f  j  j  6  (s  —  s*)  p  (s7s)  dsds‘  = 


max 


S5‘ 


(2.40) 


assuming  that 

I  p  (x/sj  rfj:  =  I — F.  ■  (2.41) 


The  nature,  of  the  optimality  of  criterion  (2.40)  is  that  for  a  given 
false  alarm  probability  F  the  mean  conditional  probability  density 
of  making  a  correct  decision  is  maximized:  s*=s 


p  f  I/’ (sVs)ls._,</s=  max . 


Based  on  the  above  proof,  and  keeping  in  mind  that  s*-s,  (2.40)  can  be 
transformed  into 

P=  J  J‘p(x/s)dx  =  max  (2.42) 

^9 


and  the  decision  making  rule  formulated  as  follows.  The  decision  s*  is 
made  if  for  all  s  G  S 


A  (s)  <  A,. 


(2.43) 
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The  critical  level  of  the  likelihood  ratio  is  selected  such  that  the 

probability  that  inequality  (2.43)  is  satisfied  assuming  realization  of 
the  event  s^  is  exactly  equal  to  the  assigned  value  1-F. 

The  decision  s*  that  a  signal  with  parameters  s  is  present  is  made  if 


A(sJ  =  sup[A(s)J. 


(2.44) 


Analogously,  we  define  the  ideal  observer  criterion  for  a  real  line 


(s«)  P  J  J  —  s) ;;  (s)  p  (s'/s)  dsds‘  = 

SS‘  , 

=  j*  ^  (s.)  P  (x/s.)  (Ik  -f-  f  p(%)p  (x/s)  (Ik  = 

**  x-x. 


;  max  • 


(2.45) 


The  nature  of  the  optimality  of  this  criterion  is  that  the  mean  (with 
respect  to  the  a  priori  distribution)  probability  of  making  a  correct 
decision  is  maximized.  With  a  uniform  a  priori  probability  distribution 
criterion  (2.45)  again  leads  to  expressions  of  the  type  (2.43)  and  (2.44). 

During  the  radar  observation  ptocess  the  situation  is  often  encountered 
in  which  the  signal  (target)  is  detected  reliably,  and  it  is  required  only 
to  provide  the  most  accurate  possible  estimate  of  the  parameters  of  the 
detected  signal.  Obviously,  in  this  case  it  is  best  to  use  the  estimate 
for  the  minimum  mean  square  error,  which  is  designated  6  below,  setting 


r(s',  s)  =  (s'-s)*  H  />(si)  =  0.  (2.^6) 

in  (2.1). 

In  this  case  optimality  criterion  (2.1)  takes  on  the  form 

P  =  *^*  =  .f /?(s)cfs  f  (s'  — s)V(s*/s)rfs*=  min  (2.47) 

or 

p  =  8*  =  f  /7  (x)  rfx  f  (s  —  s')*  p  (s/x)  ds  =  min  , 

X  s  '  (2.48) 
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where  p(s/x)  is  the  conditional  probability  density  for  the  set  of 
useful  parameters  s,  which  is  called  the  a  posteriori  probability 
function 


.  I  X  p(s)p(x/s)  p(s)p(x/s,)A(s)  (2.49) 

p(xr 


According  to  (2.48),  when  oscillation  x  arrives  at  the  system  input 
a  decision  s*  must  be  made  which  satisfies  the  equation 

X 


j"  (s  —  s*  )*  /;  (s/x)  £/s= inf  [|  (s— s*)*  p  (s/x)  rfs] .  (2.50) 


The  number  inf [4)(x, s*) ]  indicates  the  smallest  value  of  the  function 
s* 

(l)(x,s*)  taken  for  a  given  x  with  respect  to  all  possible  s*. 


2.5.  Discussion  of  Results  and  Definition  of  Object  of  Investigation 

In  all  cases,  regardless  of  whether  a  real  or  simplified  line  is 
employed,  and  regardless  of  the  optimality  criterion  selected, 
optimum  receiver  output  can  consist  of  the  likelihood  coefficient  A(s) 
or  any  other  function  of  the  set  of  useful  parameters  Y(s)  which  are 
mutually  unique  with  respect  to  A(s).  Furthermore,  the  circuit  of 
the  optimum  receiver  is  always  defined  by  identifying  the  operations 
required  to  form  from  the  received  oscillation  x(t)  an  arbitrary 
mutually  unique  likelihood  coefficient  function.  In  particular,  we 
can  emphasize  that  the  problems  of  designing  a_  receiver  for  maximum 
accuracy  [criterion  (2.47)]  and  for  maximum  detection  reliability 
[criteria  (2.21),  (2.27),  (2.40)  and  (2.45)]  are  identical.  We  also 
note  that  for  simplified  lines  the  output  A(s)  is  sufficient,  while 
the  necessary  optimum  output  of  the  receiver  is  a  discrete  sequence 
of  values  of  A(s)  for  fixed  values  of  the  parameter  s:  s^,S2,..->  s^* 

The  decision  making  rule  in  the  decision  device  is  defined  for  the 
selected  receiver  output,  and  depends  upon  a  number  of  factors:  the 
type  of  line  (simplified  or  real),  the  a  priori  probability  distribution 
p(s)  and  the  optimality  criterion  selected.  It  can  be  stated  that  by 
freeing  the  receiver  to  a  substantial  extent  of  the  need  for  making 
allowance  for  initial  assumptions,  we  transfer  all  of  the  responsibility 
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for  their  selection  and  consideration  to  the  decision  making  rule  in 
the  decision  device. 

The  decision  making  rule  for  a  real  line  consists  of  the  following  for 
the  cases  of  greatest  practical  interest  [criteria  (2.40)  and  (2.45)]. 
The  critical  level  of  the  likelihood  coefficient  (or  critical  level 

of  the  output  Yq)  is  established.  If  the  likelihood  coefficient  A(s) 

in  the  working  interval  of  S  does  not  exceed  A^  at  any  point,  it  is 

decided  that  no  signal  is  present.  If  A(s)  exceeds  Aq  on  some  set  of 

points,  it  is  decided  that  a  signal  is  present  with  parameters  s^, 

corresponding  to  the  point  with  the  highest  value  of  the  function  A(s), 
i.e.,  the  maximum  likelihood  point.  This  decision  making  scheme 
during  the  detection  process  can  be  called  a  threshold-type  scheme, 
since  some  function  of  the  set  of  non-null  decisions  is  compared  with 
a  threshold  (critical  number) . 

When  the  general  optimality  criterion  (2.1)  is  employed,  the  decision 
making  rule  can  also  be  reduced  to  a  threshold  scheme.  In  order  to  do 
this,  we  must  form  the  function  E(s*)  of  the  set  of  non-null  decisions 


•  X  (s')  =  J  [r  (s*.  s)  —  r  (s',  s)]  p  (s)  A  (s)  d%  -f 

s  (2.51) 

4-P(s,)/'(s*,  sj 

-and  define  the  critical  number 

=  P{\)  r  (s*.  s,).  '  (2.52) 

A  decision  is  made  that  no  signal  is  present  when  E(s*)  ^  E^  on  the 

set  of  all  non-null  decisions.  Otherwise  it  is  decided  that  a  signal  is 
present  with  parameters  corresponding  to  the  point  S*  at  which  the 

function  E(s*)  reaches  its  greatest  value. 

If  the  cost  of  a  miss  for  all  signals  is  the  same 


T 


0»* 
^00  » 


Sg 

S  =  Sg. 
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and  the  risk  of  all  false  alarms  is  the  same 


^.)= j 


it  is  best  to  employ,  instead  of  (2.51)  and  (2.52), 


and 


S  (s')  =  j  [r  (s*.  s)  -  p  (s)  A  (s) 


Simplified  lines  are  a  special  case  [with  condition  (2.3)]  of  real  lines. 
Therefore,  the  decision  algorithm  in  simplified  lines  can  be  obtained 
from  the  corresponding  decisions  for  real  lines  by  using  (2.3). 

Changes  in  specific  operating  conditions  lead  to  changes  in  the  a  priori 
probabilities  P(Sq)  and  p(s),  and  may  also  be  accompanied  by  change  in  the 

risk  function  r(s*,s).  Therefore,  when  the  operating  conditions  change 
the  threshold  level  Sq  must  be  changed,  and  perhaps  certain  other  parameters 

of  the  decision  devices  employed  in  automatic  radars  as  well.  In  non¬ 
automatic  radars  the  information  output  by  the  receiver  is  passed 
directly  to  a  human  observer,  rather  than  to  a  decision  computing  device. 

The  operator  makes  independent  allowance  for  specific  operating  conditions 
(risk  function  and  a  priori  probabilities)  in  accordance  with  his 
experience  and  training.  The  study  [16]  cites  a  typical  example  of  how 
a  ship's  radar  operator  will  approach  the  evaluation  of  the  information  he 
receives  differently  depending  upon  whether  the  radar  is  operating  on  a 
stormy  night  near  the  shore  or  on  a  clear  day  in  the  open  ocean.  It  is 
important  for  the  radar  designer  that  the  definition  of  the  optimum 
receiver  remain  the  same  in  this  case  as  well.  In  other  words,  the 
circuit  of  an  optimum  receiver  does  not  depend  only  upon  specific 
operating  conditions,  but  also  upon  whether  the  radar  is  automatic  and 
whether  decisions  are  made  by  a  human  observer. 

It  is  assumed  in  general  below  that  the  system  incorporates  a  decision 
device  which  operates  on  the  basis  of  the  generalized  Neumann-Pearson^ 
criterion,  although  this  criterion  is  not  the  best  one  for  all  operating 
conditions. 

When  the  parameters  of  the  detected  signal  are  estimated  on  the  basis  of 
the  minimum  mean  square  error  criterion,  the  decision  rule  is  to  use  as 
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the  estimate,  as  follows  from  (2.50),  the  "center  of  gravity"  point  of  the 
a  posteriori  probability  function  p(s/x)  or  the  function  p(s)A(s).  For 
this  reason,  the  signal  parameter  estimates  obtained  for  detection  based 
on  the  Neumann-Pearson  and  ideal  observer  criteria  have  the  minimum  possible 
mean  square  error  only  when  the  maximum  likelihood  point  coincides  with  the 
"center  of  gravity"  of  the  function  p(s)A(s). 

We  have  represented  the  radar  observation  process  as  an  experiment  whose 
outcome  is  random,  which  makes  it  necessary  to  choose  between  competing 
hypotheses  and  to  estimate  certain  parameters  of  the  observed  set  of 
random  quantities  x.  The  subject  of  our  study  will  be  optimal  circuits 
and  those  performance  indicators  of  radar  systems  whose  theoretically 
limiting  values  are  bounded  by  the  statistical  nature  of  the  problem 
resulting  from  the  occurrence  of  interference.  These  performance 
indicators  include  the  threshold  signal  or  threshold  signal/noise  ratio, 
the  accuracy  with  which  the  useful  signal  parameters  are  estimated  and 
the  resolution. 

The  threshold  signal/noise  ratio  is  the  minimum  signal/noise  ratio  required 
for  the  system  to  operate  with  the  required  reliability.  The  quantitative 
measure  of  reliability  is,  in  general,  the  mathematical  expectation  p  of 
the  risk  function  r.  I-Jhen  the  Neumann-Pearson  criterion  is  employed,  this 
measure  is  the  maximized  mean  correct  detection  probability  and  a  fixed 
false  alarm  probability. 

The  accuracy  with  which  useful  signal  parameters  are  estimated  is  measured 
in  terms  of  the  mean  square  error  or  dispersion  of  the  estimates. 

Resolution  is  estimated  by  the  dimensions  of  the  region  in  the  useful 
parameter  space  S  such  that  it  is  practically  impossible  to  resolve  two 
signals  with  useful  parameters  belonging  to  that  region. 

Statistical  Interpretation  of  the  radar  observation  process  makes  it 
possible  in  principle  to  discuss  and  solve  a  broad  group  of  problems 
involved  in  radar  theory.  These  problems  include  the  following: 

1.  Optimum  methods  for  processing  received  oscillations  in  receiver 
(optimum  receiver  circuit) . 

2.  Optimum  decision  making  rule  (optimum  decision  device  circuit). 

3.  Quantitative  estimate  of  potential  capabilities  of  radar  systems  in 
terms  of  basic  performance  indicators. 
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4.  Influence  of  parasitic  parameters  on  potential  capabilities  of 
system  and  on  optimum  processing  circuit. 

5.  Rational  choice  of  radar  signal  form. 

6.  Influence  on  system  performance  indicators  of  deviations  from 
optimum  of  processing  method  and  decision  making  rule. 

The  investigation  below  is  limited  to  examining  the  most  practically 
important  case  —  that  of  reception  against  the  background  of 
additive  random  interference.  The  analytical  expressions  for  the 
functions  p(x/s)  and  A(s)  which  figure  in  the  calculations  are  defined 
as  follows.  The  probability  density  p(x/s,u)  of  the  quantity  x,  assuming 
that  the  signal  3(t,s,u)  is  input,  is  obviously  the  same  as  the 
probability  density  of  the  interference  p(n)  if  the  value  of  the 
interference  n(t)  is  assumed  to  be 


n{t)  =  x{t)-3{t,  8.  a). 


(2.53) 


Thus,  for  the  class  of  random  interference  whose  probability  distribution 
is  defined  through  (1.79), 

.54) 

where  (tj^,!^)  is  an  arbitrary  interval  which  includes  the  time  interval 
on  which  the  function  3(t;s;u)  is  non-zero. 

Analogously,  assuming  the  realization  of  the  event  Sq, 

I, 

p(x/s,)  =  Aexp|^Jjc*(Od(|'.  (2.55) 


jp(x/s,  «)  =  feexp|— {x(0  — ■5((.  s,  tt)]’rf/|, 


(2 


and  the  likelihood  coefficient  A(s,u)  for  the  set  of  useful  s  and 
parasitic  u  parameters  of  the  signal  will  be 
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/• 


(2.56) 


The  likelihood  function  p(x/s)  and  the  likelihood  coefficient  A(s) 
of  the  set  of  useful  parameters  are  obtained  by  statistical  averaging 
of  (2.54)  and  (2.56)  with  respect  to  all  possible  values  of  u 


P{^ls)  =  ^  Pin)  p{x/s,  u)du, 


(2.57) 

(2.58) 


We  note  that  according  to  (2.58)  it  is  completely  impossible  to  relieve 
the  receiver  of  the  need  for  allowing  for  a  priori  data,  since  formation 
of  the  likelihood  coefficient  A(s)  requires  that  the  a  priori  distribution 
of  the  parasitic  parameters  p(u)  be  available  or  given. 

In  order  to  simplify  the  investigation  and  establish  separately  the 
influence  of  each  of  the  factors  on  the  radar  observation  process,  we 
shall  begin  our  study  with  the  simplest  cases,  and  gradually  make  the 
circuit  more  complex  as  we  draw  nearer  to  real  operating  conditions. 


Chapter  4.  Estimation  of  Signal  Parameters 

4.1.  Discussion  of  Problem  and  Methods  for  Solution 

The  problem  of  estimating  signal  parameters,  which  arises  in  radar  as 
well  as  related  areas  of  electrical  engineering  —  radio  navigation, 
radio  remote  control,  radio  telemetry,  among  others  —  is  formulated 
as  follows.  It  is  well  known  that  the  oscillation  x(t)  which  arrives  at 
the  input  of  a  receiver  consists  of  an  additive  mixture  of  interference 
n(t)  and  signal  s(t,s,u),  the  set  of  useful  parameters  of  which  s  belongs 
to  the  interval  S.  Given,  or  defined  by  introducing  certain  assumptions, 
are  the  a  priori  distribution  functions  of  the  useful  and  parasitic 
parameters  p(s)  and  p(u),  as  well  as  the  probability  density  functional 
p(n)  of  the  interference  n(t).  The  problem  calls  for  defining  the  system 
operator  (processing  circuit  and  decision  making  rule)  which  guarantees 
the  best  (in  terms  of  mean  square)  estimate,  and  calculating  the  numerical 
characteristics  of  the  parameter  sampling  accuracy  —  the  mean  square  errors. 
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One  more  initial  condition  must  be  added  to  this  general  statement  of 
the  problem  for  the  case  of  sampling  the  parameters  of  radar  signals. 

The  problem  of  the  most  accurate  possible  sample  of  the  signal  parameter 
arises  after  the  fact  that  a  target  is  present  has  been  registered 
reliably.  We  can  therefore  assume  that  when  the  parameters  are 
sampled  the  signal  is  significantly  stronger  than  the  interference 
(q2»l) ,  at  least  by  enough  for  reliable  target  detection.  We  note 
that  there  is  no  limit  placed  on  the  input  signal/noise  power  ratio. 

The  latter  condition  is  significant.  The  output  Y(s)  on  the  set  of 
possible  decisions  contains  a  large  number  of  isolated  maxima  (spikes) , 
one  of  which  is  caused  by  the  signal-noise  mixture,  with  the  remaining 
being  caused  by  interference  alone.  If  the  signal  energy  is  too  low, 
there  is  ambiguity  in  the  detection  and  sampling  —  any  large  Interference 
spikes  on  the  interval  S  may  be  Interpreted  as  a  spike  caused  by  a  signal, 
since  the  presence  of  a  signal  has  been  stipulated  in  advance.  Under 
these  conditions,  it  is  meaningless  to  solve  the  problem  of  radar  signal 
parameter  estimation  accuracy. 

The  problem  of  estimating  signal  parameters  and  the  problem  of  signal 
detection  have  a  common  statistical  model,  which  was  described  in  §2.1. 

The  difference  between  the  estimation  problem  and  the  detection  problem, 
which  also  accompanies  parameter  estimation,  consists  of  differences  in 
the  risk  function  and  in  the  initial  assumptions  with  respect  to  the 
a  priori  probabilities,  which  essentially  reflects  a  difference  in  the 
tactical  content  of  both  problems.  The  problem  of  parameter  estimation 
is  usually  associated  with  the  target  tracking  process  or  with  the  process 
of  observing  the  trajectory  of  a  detected  target. 

The  second  chapter  examined  the  problem  of  parameter  estimation  in  very 
general  features.  The  optimal  system  operator  was  defined,  and  analytical 
formulas  (2.47)  and  (2.48)  for  mean  square  errors  were  presented.  If  the 
estimate  of  s*  is  obtained  in  expressions  (2.47)  and  (2.48)  by  applying 
the  optimum  operator,  these  expressions  determine  the  theoretically 
limiting  or  potential  mean  square  errors. 

The  objective  of  the  continued  investigation  in  the  present  chapter  is 
a)  to  establish  certain  general  regularities  which  occur  in  estimating 
radar  signal  parameters;  b)  to  analyze  potential  mean  square  errors; 
c)  to  define  a  method  for  constructing  practical  circuits  for  obtaining 
best  estimates. 

For  definition  we  shall  assume  that  the  set  of  useful  parameters  s  consists 
of  the  two  scalar  parameters  a  and  3 
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(4.1) 


's  =  (a,  p)  and  rfs  =  dadp.  . 


Instead  of  the  mathematical  expectation  of  the  square  of  the  absolute 
value  of  the  error  vector  defined  through  (2.47)  and  (2.48)  we  shall, 
in  accordance  with  practical  requirements,  be  examining  the  mean  square 
errors  <6^>  and  of  the  estimate  of  each  of  the  scalar  parameters 

a  and  3 •  The  optimum  system  must  simultaneously  minimize  the  mean  square 
values  of  both  of  the  error  vector  components  a-a*  and  3-3*  (  and 

3*  —  estimates  of  the  parameters  a  and  3),  i.e.,  it  must  simultaneously 
minimize  ■(3^^  and  ^3^).  ^'Hien  necessary,  the  discussions  below  can  be 

extended  to  sets  of  large  numbers  of  useful  parameters. 

There  are  at  least  two  possible  different  methods  for  calculating  mean 
square  errors.  The  first  method  is  based  on  representation  (2.47),  in 
accordance  with  which  the  mean  square  errors  and  are  defined 

by  the  formulas 


(Sl  )  =  Jp(s)dsfp (s'/s)  («•  -  a)* ds\ 
s  s* 

0l)=jp (s) dsj  p (s'/s)  (P*  —  p)* ds\  •  2 ) 


The  structure  of  the  set  S*  is  the  same  as  that  of  set  S.  The  analysis 
is  done  as  follows .  An  ensemble  of  cases  is  examined  in  which  the 
received  oscillation  x(t)  represents  the  sum  of  interference  n(t)  and 
signal  with  fixed  values  of  the  useful  parameters  a  and  3  or  s.  The 
estimate  s*=(a*,3*)  obtained  by  applying  some  operator,  not  necessarily 
optimum,  is  random.  The  conditional  mean  square  errors  are  first  defined 
for  a  fixed  value  of  s=(a,3) 


(®l]3  =  J  (a*  — a)*P(s7s)^s', 

s* 


[»h.=S  ^ 

s*  •  /  • 


(4.3) 


after  which  averaging  is  done  with  respect  to  s 
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(«“.)=]/’ (s)  l>;  ).*.  '(s; ) = J  />  (s)  (!'  I.  rfs. 


(4.4) 


If  for  all  s  G  S  the  errors 


[6^1  and  [6^]  are  independent  of  s,  then 
as  p  s 


(8p)  =  I8;]s 


(4.5) 


and  there  is  no  need  to  average  in  (4.4). 

The  result,  of  course,  remains  unchanged  if,  fixing  both  the  useful 
parameters  s  and  the  parasitic  parameters  u,  we  find 

f  (a--ay  p(s%  u)ds\ 


after  which  we  average  with  respect  to  s  and  u. 


The  second  method  for  analyzing  mean  square  errors  is  based  on 
representation  (2.48),  according  to  which 


<«;  )=|  j (P-PT  ; 


p  (x)  (/x  \  (a  -7  a')*  P  (s/x)  d^, 


(4.6) 


An  ensemble  of  cases  is  examined  in  which,  regardless  of  the  useful 
parameters  s  =  (a, 3)  of  the  "transmitted"  signal,  the  input  of  the 
receiving  device  receives  the  same  oscillation  s(t),  so  that  the  same 
estimate  s*=(a*, 3 *)  is  obtained.  The  conditional  mean  square  errors 
are  first  determined  for  fixed  values  of  x  and  s* 
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(4.7) 


S 

5 


after  which  we  average  with  respect  to  x 


(«')=]  (sj)=j/;(x)[6j];rfx. 


(4.8) 


If  the  errors  [6^]  and  [6^]  are  independent  of  x, 

OC  X  {ij  X 


(4.9) 


and  it  becomes  unnecessary  to  average  in  (4.8).  It  was  essentially  this 
statement  for  which  Kotel'nikov  [18]  solved  the  problem  for  signals 
containing  one  useful  parameter  and  no  parasitic  parameters.  This 
analytical  method  is  also  employed  in  [8,  15],  but  without  sufficient 
foundation.  The  investigation  below  shall  employ  both  methods  for 
calculating  mean  square  errors.  In  practical  calculations  preference 
will  be  given  to  the  first  method,  which  produces  comparatively  simple 
analytical  expressions. 

In  some  studies  [31]  the  calculation  of  limiting  accuracy  reduces  to 
calculating  the  dispersion  of  the  effective  estimates.  Kramer  [19] 
presents  a  proof  that  when  there  is  one  scalar  random  parameter  a  the 
mean  square  error  in  its  estimate,  corresponding  to  the  quantity  [<5^1 

which  we  defined  by  means  of  (4.3),  is  bounded  from  below  by  the  number 


(4.10) 
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It  is  assumed  that  the  derivative  figuring  in  formula  (4.10)  exists,  and 
that  the  estimate  is  unbiased,  i.e.. 


(4.11) 


and  that  the  error  dispersion  can  be  examined  instead  of  the  mean 
square  errors. 

The  estimate  a*  for  which  the  dispersion  reaches  the  smallest  possible 
value  defined  by  the  formula  cited  above  is  called  the  effective  estimate. 
An  effective  estimate  exists  [i.e.  (4.10)  as  an  equals  sign]  when  and  only 
when  these  two  conditions  are  satisfied: 


a)  p(x/a)  =  <j>(x)/j(a7a), 


(4.12) 


where  p(a*/a)  is  the  distribution  function  of  the  estimates  for  a  given 
fixed  value  of  a,  4)(x)  is  an  arbitrary  function  of  x; 


6)  -^ln(/;(a7a)  =  A(a*  — a), 


(4.13) 


where  k  is  a  coefficient  independent  of  a*,  but  which  may  depend  upon  a. 

When  there  is  a  large  number  of  random  parameters  a,  there  are 

analogous,  but  more  complicated,  formulas  for  the  dispersion  of  the 
effective  estimate  and  for  the  existence  condition  of  effective  estimates. 
It  follows  from  the  very  nature  of  conditions  a)  and  b)  that  effective 
estimates  can  be  expected  to  exist  only  in  certain  special  cases. 
Therefore,  when  investigating  an  arbitrary  system  intended  for  sampling 
random  parameters,  analysis  of  the  dispersion  of  the  effective 
estimates  by  formulas  analogous  to  (4.10)  still  do  not  provide  the  values 
of  the  dispersion  of  the  actual  estimates,  or  even  the  theoretically 
limiting  dispersion  of  the  estimates  for  the  system  in  question. 

The  estimates  (a*,  3*)  can  be  obtained  through  various  methods  which, 
in  the  general  case,  yield  different  deviations  from  the  theoretically 
limiting  estimates.  The  most  important  of  the  possible  estimation  methods 
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is  the  maximum  likelihood  method,  according  to  which  the  estimate 
(a*,  B *)  is  the  point  at  which  the  a  posteriori  probability  function 

p(a,  3/x)  or,  if  the  a  posteriori  distribution  of  p(a,3)  can  be  assumed 
uniform,  the  likelihood  coefficient  A (a, 3),  reaches  its  maximum  value. 

The  estimation  method  formulated  above  coincides  with  the  decision  making 
rule  presented  in  §2.4.  This  method  is  important  because  in  all  cases  in 
which  an  effective  estimate  exists  it  can  be  obtained  by  the  maximum 
likelihood  method  [19].  Another  advantage  of  the  maximum  likelihood 
method  over  many  others  is  the  fact  that  the  maximum  likelihood  point  is 
invariant  with  respect  to  arbitrary  mutually  unique  transformation  of 
the  output.  This  makes  it  easier  to  construct  the  processing  circuit  for 
obtaining  the  maximum  likelihood  estimate. 


4.2.  General  Principles  of  Parameter  Estimation 

In  order  to  study  some  of  the  principles  occurring  in  the  estimation  of 
radar  signal  parameters,  we  shall  employ  a  method  for  calculating  mean 
square  errors  in  which  an  ensemble  of  cases  with  a  fixed  value  of  the 
received  oscillation  x(t)  is  examined. 

The  a  posteriori  probability  density  for  the  useful  parameters  p(a,3/x) 
is 


p{A,  P/x)  =  *,/7(a,  p)A(a,  p),  (4.14) 


where  p(a,3)  and  A(a,3)  represent  the  a  priori  probability  density  and 
likelihood  coefficient  for  the  useful  parameters,  and  k^  is  the 

proportionality  coefficient,  which  in  general  depends  upon  the  received 
oscillation  x(t).  The  numerical  value  of  the  coefficient  k^  can  be 
found  from  the  normalization  conditions 


p/x)</eu/p=  1. 

s  . 


(4.15) 


When  estimating  signal  parameters  under  reliable  detection  conditions 
we  are  interested  in  the  behavior  of  the  function  p(a,3/x),  and  consequently 
the  function  p(a,3),  only  in  the  vicinity  of  a  powerful  spike  caused  by  the 
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presence  of  a  signal.  We  shall  call  this  region  the  vicinity  of  the 
estimate  and  designate  it  a.  Obviously,  it  can  almost  always  be 
assumed  that  the  distribution  of  p(a,3)  is  near-uniform  in  the  vicinity 
of  the  estimate.  Accordingly,  we  shall  assume  below 


/>(a,  p/x)=^;A(a.  p). 


(4.16) 


When  a  priori  information  is  available  which  can  modify  a  posteriori 
distribution  (4.16)  significantly,  it  must  be  allowed  for  by 
substituting  the  basic  expression  (4.14)  for  (4.16). 

A  decision  device  based  on  studying  the  output  Y(a,3)»  which  is  a 
mutually  unique  function  of  the  a  posteriori  distribution  (4.16), 
outputs  in  some  fashion  the  estimate  a*,  3*  obtained.  The  conditional 
mean  square  errors  for  given  values  of  x(t)  and  (a*,  3*),  according  to 
(4.7),  are 


(4.17) 


In  examining  expressions  (4.17),  we  can  reach  some  practical  conclusions. 

If  the  signal  contains  several  random  parameters,  such  as  the  parameter  a 
and  others,  the  mean  square  error  in  the  estimate  of  the  parameter  a  is 
independent  of  whether  this  parameter  is  estimated  in  conjunction  with  the 
others,  or  if  all  of  the  other  parameters  are  perceived  by  the  system  as 
parasitic  and  information  regarding  them  is  destroyed  during  the  processing 
(by  integration) . 


For  every  given  value  of  x  the  minimum  possible  value  of  the  errors 
and  [63]  ,  and  consequently  of  the  mean  square  errors,  (4.8),  is  obtained 

P  X 

when  the  estimates  of  the  parameters  are  the  coordinates  of  the  "center  of 
gravity"  of  the  a  posteriori  distribution  p(a,3/x).  The  latter. 
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of  course,  is  the  same  as  the  result  obtained  in  the  second  chapter,  if 
we  use  as  our  basis  the  minimum  mean  square  of  the  absolute  value  of  the 
error  vector. 


The  best  estimate  of  the  useful  parameters,  which  has  potential  mean 
square  error,  is  thus  the  estimate  with  respect  to  the  "center  of  gravity" 
of  the  a  posteriori  distribution. 


The  estimate  a*,  obtained  by  the  maximum  likelihood  method,  in  the 


general  case,  in  contrast  to  the  estimate 


3  ,  is  biased,  and  has 

eg 

mean  square  errors  which  exceed  the  potential  values.  The  maximum 
likelihood  estimate  is  an  unbiased  best  estimate  only  when  the  coordinates 
of  the  center  of  gravity  and  the  maximum  likelihood  points  of  the 
a  posteriori  distribution  of  p(a,3/x)  coincide 


(4.18) 


The  sufficient  condition  for  ensuring  equality  (4.18)  is  the  condition  of 
symmetry  on  the  plane  (a, 3)  of  the  lines  of  the  level  of  function  p(a,3/x) 
with  respect  to  the  maximum  likelihood  point  (oi^>  This  condition 

can  be  written  analytically  as 


P  (%  +  PI  -f  P7x)*=  P  (a*  —  a,  p’  —  p/x). 


(4.19) 


In  fact,  when  (4.19)  is  satisfied 


a  b  ^  ^ 

*0  "I"  J  * J  “I" ®  •  Po  +  P/*)*^?  ~ ®0  • 

—a  .  —6 
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since  the  function 


V  ■  V 

j  P  («0  +  ®-  Po  +  P/x)  J  P  (^0  -  “•  Po  -  *P/x)  4  = 

-b  i-fc  ,  . 

*  ■  . 

=Jp(«o-“*Po+p7x)4  ‘  - 


is  an  even  function  of  a.  The  symbols  a  and  b  define  the  intervals  of 

variation  of  the  increments  a  and  3  of  the  parameters  a  and  3>  respectively. 

The  equality  of  35  and  3  is  obtained  analogously  under  condition  (4.19) 

0  eg 

If  the  a  posteriori  probability  function  p(a,3/x)  in  the  vicinity  of  the 
maximum  likelihood  estimate  (a*,  3q)  is  analytical,  condition  (4.19)  means 

that  it  can  be  represented  as  the  series 


P  (a.  PA)  =  S  «!/  A  -  “o  )'  (P  —  Po  )'• 

U  ■' 

(l-^  j  —  2ni  i,  /,  rt=0,  1,  2,,,.), 


(4.20) 


All  of  the  coefficients  a  of  series  (4.20)  for  which  the  sum  of  the 

i » J 

powers  i+j  is  an  odd  number  are  zero.  Representation  (4.20)  follows 
directly  from  the  fact  that  according  to  (4.19) 


p{^.  pa)=4-(^K+a-<*o).  Po'+(p-p;)A]+ 

+  P  («5’—  A  —  “o  )>  Po  -  (P  —  Po  )/x]}-  ■  - 


(4.21) 


We  note  that  in  conditions  (4.19)  and  (4.20)  above,  which  are  sufficient 
for  the  maximtam  likelihood  point  to  coincide  with  the  center  of  gravity 
of  the  distribution  of  p(a,3/x),  the  a  posteriori  probability  function  can 
be  replaced  with  any  other  function  which  is  unique  with  respect  to 
p(a,3/x).  These  conditions,  when  necessary,  can  be  extended  to  the  case  of 
a  large  number  of  useful  parameters.  Further  examination  of  the  question  of 
the  potential  accuracy  of  sampling  of  signal  parameters  requires  more 
detailed  study  of  the  a  posteriori  probability  function. 
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Of  greatest  interest  for  the  practical  applications  encountered  below 
is  the  case  in  which  the  signal  contains,  besides  the  useful  parameters 
a  and  3,  two  parasitic  random  parameters:  the  initial  phase  (j)  and 
intensity  e,  i.e.,  when  the  received  signals  are  of  the  form 


ej(/.  a.  p.  9)  =  .Re{S(/.  a.  P) (4.22) 


We  also  assume  that  the  interference  has  the  nature  of  white  noise,  or 
of  a  normal  process  with  a  very  wide  spectrum,  so  that  the  probability 
density  functional  of  the  interference  p(n)  can  be  represented  by  means 
of  (1.79). 

The  a  posteriori  probability  function  for  all  random  parameters  is 

'  ’  f  CO 

>(«.  p.  •.  =  ^  J  a,  P,  ¥)rf/  + 

C  L  -OO 

00 

■ j  P, 

— oo 


■  A(«,  P.  •»  T/x)  =  ^,^^exp  |— •V(a,  p)-f 

~|-j^  J  l'X(0 — ®0  *  Po  *  P»?)^^"l" 

—CO 

I  3{t,  a* ,  P* ,  <Po)J(^  a,  p^  ^)dt 


where  (ag,  3g,  <j)g,  £§)  is  the  maximum  likelihood  point,  found  from  the 
system  of  equations: 


(4.24) 


(4.23) 

I 
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Jl_ 

df 


%MJ 

In/? (a,  p,  ®,  ¥//)“  ^  [-^(O  *0  •  Po  •  “iPo  ^ 

—00  »  . 

®o  •  Po  •  9o)^^  =  0/ 

«  , 

lnp(a,  P,  e,  ?/x)=  j  Po*> 


•  • 


df 


■  .  X-^-S(^  «0.  Po.  =  ■ 

•  .  ■  • 
In  p  (a.  p,  *,  f/x)  =  j  (-« (0  —  *0  !» “o  •  Po  *  fo  )1  X 

f-OO 

Po*  9o)d(  =  0> 


TO 

-|j:ln/?(a.  P.  e,  v/x)=  J  [A-(0-*0'»(^  »o*  PI  »  )1  X 

t-00 

'  .  i  X  J(/.‘  Po»  = 


(4.25) 


The  second  term  in  the  braces  in  (4.24)  can  be  omitted  in  the  vicinity  of 
(a*,  3*).  This  possibility  is  determined  by  two  factors.  First  of  all,  the 

function  [x(t)  -  e*  3(t,  a*,  3§,  4>§)]  represents  interference  with  the 

minimum  possible  energy  for  given  x(t)  superimposed  on  the  signal. 
Accordingly,  the  mean  square  value  of  the  second  term  for  any  a  and  3 
does  not  exceed  /2q,  and  under  reliable  detection  conditions  is 
significantly  smaller  than  the  amplitude  value  of  the  third  term  near 

(a§,3§). 
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Second,  the  second  term  in  (4.24)  and  all  of  its  first-order  partial 
derivatives  at  the  point  (a*,  3*,  (j)*)  are  zero  according  to  (4.25). 

Furthermore,  all  of  the  higher-order  partial  derivatives  with  respect 
to  (p  are  zero  as  well,  since 


-  r)* 

a.  P.  =  a,  p,  (p); 


(n,  A  =  0,  I,  2,...). 


(4.26) 


Therefore,  when  the  second  term  in  (4.24)  is  expanded  into  a  series 
arranged  by  powers  of  (a-a*) ,  (3-3q)»  only  those  terms  are 

retained  which  contain  the  higher  powers  of  the  small  parameters  (a-a*) 

and  (3-3*).  Thus,  the  second  term  in  (4.24)  at  the  point  (a*,  3^)  is 

equal  to  zero,  and  is  very  small  in  the  vici^ty  of  that  point  (the  mean 
square  value  is  significantly  smaller  than  /2q).  Furthermore,  it  is 
approximately  the  case  that 


pi  *.  tf/x)=:k  ^ 


f^exp|-BV(a,  P)  + 

2««S  r  . 

+  J  Po» 


(4.27) 


Elimination  of  the  parasitic  parameters  yields 


'  p  (a,  p,  e/x)  =  V  (*) 
•  1+?* 

1  ■ 

X/ 


r  • 

f  f 

•  at.  ,J 

.  — o 


S(tt  ®i  P)^(^»  ®o  *  Po  ) 
/?(«,  P/X)  =  *,®  P)* 


J  S^t,  a,  p)'S(/f  *0  »  Po  ) 


dt 


(4.28) 


(4.29) 
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where  $  -  a  certain  function  determined  by  the  intensity  distribution 
p(e)*  For  a  Rayleigh  distribution  of  the  parameter  c 


w 

j  S{t,  a,  p)  X 


(4.30) 


Outside  the  region  a,  which  by  definition  is  the  vicinity  of  the  estimate, 
expressions  (4.29)  and  (4.30),  of  course,  cannot  serve  as  an  approximation 
of  the  a  posteriori  distribution.  However,  under  reliable  detection 
conditions  the  probability  concentrated  outside  the  region  O  under  the 
functions  (4.28)-(4.30) ,  and  under  the  corresponding  exact  expressions 
for  the  a  posteriori  probability  function,  is  negligibly  small.  Therefore, 
in  order  to  calculate  the  mean  square  errors  by  formulas  (4.17)  with  good 
enough  accuracy  for  practical  applications,  the  coefficient  k  can  be  defined 

by  substituting  the  expressions  obtained  for  the  function  p(ct,3/x)  in  the 
vicinity  of  the  estimate  in  normalization  condition  (4.15).  The  coefficient 
k  is  uniquely  defined  as  a  function  of  the  estimate  (a*.  6*.  £*) . 

Our  analysis  establishes  that  the  a  posteriori  probability  function  p(c(,6/x), 
or  any  other  optimal  output  Y(a,6)  in  the  vicinity  of  the  estimate  a 
always  has  the  same  appearance,  i.e.,  it  is  always  some  known  function  (}) 


P  (a.  P/x)  =  9  (a,  p,  flo  •  Po  »  *0  )  (4.31) 


of  the  variables  a,  3  and  the  estimates  a*,  3*,  £*,  and  does  not  depend 

directly  upon  the  received  oscillation  x(t).  For  signals  with  fixed 
intensity  the  parameter  e*  in  the  latter  expression,  of  course,  must  be 
omitted.  ,  ^ 


We  conclude  on  the  basis  of  (4.31)  that  in  the  general  case  the  coordinates 
of  the  "center  of  gravity"  a  3^„  of  the  a  posteriori  distribution,  which 

are  the  estimates  with  the  theoretically  limiting  mean  square  deviation. 
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can  be  represented  as  certain  known  functions  (J)^  and  tj)2  of  the  maximum 
likelihood  point  (a*,  eg)’ 

“eg”  II  *■  Po  »  *0  )  =  Vi  (“o  *  Po  «  *0  )• 

■  (4.32) 

■  JJ  P*  ®o  »  Po  •  ®o  )  =  9i  (^0  -  Po  »  «o  )• 

s 

Therefore,  in  the  general  case  obtaining  an  estimate  which  realizes  the 
potential  capabilities  of  the  system  with  respect  to  accuracy  can  be 
reduced  to  obtaining  the  maximum  likelihood  estimate  of  a*,  3g,  eg* 

Furthermore,  both  the  useful  parameters  and  the  intensity  parameters  are 
subject  to  estimation. 

We  note  also  that  the  dispersion  of  the  estimates  ag,  3g  is  always  the 

same  as  the  dispersion  of  the  best  estimates  a  ,  3  »  since  for  every 

eg  ug 

given  value  of  x(t)  it  is  determined  by  the  dispersion. of  the  a  posteriori 
distribution  p(a,3/x). 

In  the  overwhelming  majority  of  cases  with  which  we  must  deal  when 
sampling  the  parameters  of  radar  signals,  the  absolute  value  of  the 

integral  J  S(f.  a,  P)S((.a-}-  a,  does  not  depend  upon 

the  parameters  a  and  3>  but  is  a  function  only  of  the  increments  a  and  3 

J  S(t,  a,  fi)S(f.  a-fo,  P-|-P)d(  =  2Q’'F  (i,  p).  (4.33) 

—00 

Assuming  in  the  left  part  of  the  latter  equality  that  one  time  a=3=0,  and 
the  second  time  a=  -d,  3=  -3,  and  keeping  in  mind  that  under  condition 
(4.33)  both  substitutions  should  produce  the  same  result,  we  obtain 

W  (a,  p)  =  V  (-  -  P).  •  34) 
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Equality  (4.19)  which,  in  turn,  leads  to  (4.20),  will  hold  for  a  posteriori 

probability  functions  (4.29)  and  (4.30).  Thus,  if  the  modulation  function 

S(t,a,3)  of  the  received  signal  for  (a,B)6a  obeys  condition  (4.33),  the 

maximum  likelihood  estimate  (oi*,  3*)  is  unbiased  (coincides  with  the  "center 

of  gravity"  a  ,  3  )  and  directly  yields  the  theoretically  limiting 

eg  eg 

mean  square  deviation  from  the  true  values  of  the  measured  parameters. 


For  cases  in  which  the  signal  exceeds  the  interference  by  a  significantly 
greater  amount  (q^»l)  the  errors  (ag-a^)  and  (3^-3)  in  the  sample  of  the 

signal  parameters  are  small,  and  regardless  of  whether  or  not  condition 
(4.53)  is  satisfied  the  optimum  output  Y(a,3)  can  be  represented 
approximately  as  the  series 


K  (a.  P)  =  In  /J  (a,  p/x)  =  a„  (a  —  aj  )*  -f- 

+^ii(*“®o)(P“  Po)»  (®»  p)€®* 


which  retains  only  the  terms  with  low  powers  of  the  small  parameters. 

The  coefficients  of  the  series  are  independent  of  a  and  3  and,  under 
condition  (4.33),  of  a^*  and  3q*  as  well.  The  possibility  of  representing 

the  function  In  p(a,3/x)  in  the  form  (4.35),  which  is  equivalent  to 
approximating  the  a  posteriori  probability  function  by  a  normal 
distribution,  guarantees  the  existence  of  an  effective  unbiased  estimate 
of  the  parameters  of  the  received  signal  [19].  In  other  words,  when  the 
errors  are  essentially  small  the  maximum  likelihood  estimate  is  nearly 
always  an  effective  estimate. 

Consequently,  if  condition  (4.33)  occurs  or  representation  (4.35)  can  be 
used,  the  estimates  obtained  during  detection  based  on  the  generalized 
Neumann-Pearson  criterion  or  ideal  observer  criterion  are  the  best 
estimates. 

The  conditional  mean  square  errors  snd  calculated  for  a  fixed 

value  of  x(t),  according  to  (4.31),  depend  upon  a*,  3^»  in  the  general 
case.  Consequently,  calculation  of  the  mean  square  errors  ,  ^<5^^ 

in  the  general  case  can  be  reduced  to  averaging  [6^]  and  [6|]  with  respect 

Ct  P  X 

to  a*,  3q,  c*.  If  condition  (4.33)  is  satisfied,  expression  (4.31)  for  the 
a  posteriori  probability  becomes  the  following: 


p{a.  p/x)  =  ,9(e;,  a-aj,  P-pJ).  (4.36) 
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The  conditional  mean  square  errors  [6^]  and  [5q]  of  the  best  estimate 

Ot  X  p  X 

depend  only  upon  the  parameter  £*.  Accordingly, 

✓ 


and  only  for  condition  (4.33)  and  a  fixed  intensity  of  the  received 
signals  does  equality  (4.9)  occur 


Finally,  an  important  conclusion  which  follows  directly  from 
representation  (4.31)  is  that  the  obtaining  of  an  estimate  is  not 
necessarily  accompanied  by  reproduction  of  the  output  Y(ot,3)  for  all 
values  of  ot  and  3  in  region  O,  Fundamentally,  if  the  output  is  defined 
in  the  region  a  at  three  points,  such  as  (a^,  3j^)>  (ot2>32) 

the  solution  of  the  system  of  equations 


?  (®1»  Pl»  ®0  •  Po'*  *0  )  —  ^(®l»  Pt)» 

?  (®1>  P»*  *0  *  Po  •  ®0  )  “  ^  (®**  - P*^’ 

t  («a.  ?»•  «o  *  Po  >  ®o  )  =  ^  (?J?  Pi) 


(4.38) 


makes  it  possible  to  find  the  estimate  a^*,  3q*  and,  when  necessary,  the 

parameter  Eq*.  When  obtaining  estimates  by  solving  a  system  of 

equations  here  and  below  we  assume,  in  accordance  with  real  operating 
conditions,  that  the  intensity  of  the  received  signal  is  a  parasitic 
random  parameter.  Other  quantities,  such  as  the  derivatives  of  the 
output,  can  be  employed  instead  of  the  output  in  system  (4.38).  It  is 
necessary  only  that  for  random  intensity,  which  is  perceived  as  a 
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parasitic  parameter,  the  number  of  equations  exceed  by  one  the  number  of 
measured  parameters  and  that  the  system  have  a  unique  solution  in  the 
vicinity  of  the  estimate. 

The  investigation  could  have  been  done  in  more  general  form  without 
specifying  the  physical  nature  of  the  parasitic  parameters.  However, 
the  practical  applications  of  the  theory  which  are  encountered  below  do 
not  require  this.  In  addition,  this  would  have  complicated  the 
investigation  significantly. 

4.3.  Analysis  of  Dispersion  of  Estimates  Obtained  by  Maximum 
Likelihood  Method 

The  preceding  section  indicated  that  the  maximum  likelihood  estimate  has 
theoretically  limiting  (potential)  dispersion.  There  is  no  bias  in 
the  maximum  likelihood  estimate,  since  condition  (4.33)  is  usually 
satisfied,  or  representation  (4.35)  can  be  used.  However,  even  when  the 
estimate  is  biased,  allowance  can  be  made  for  this  during  sampling. 
Therefore,  the  potential  capabilities  of  the  system  in  terms  of 
accuracy  can  be  determined  by  the  scattering  or  dispersion  of  the 
estimates  (cXq*,  3q*)  obtained  by  the  maximum  likelihood  method. 

Our  next  problem,  accordingly,  is  to  obtain  more  or  less  simple 
analytical  expressions  for  the  dispersion  of  the  maximum  likelihood 
estimates  (a*,  3*).  It  is  best  to  use  the  analytical  method  which 

examines  an  ensemble  of  received  oscillations  x(t)  with  fixed  signal 
parameter  value. 

It  should  be  noted  that  dispersion  is  not  the  only  possible  scattering 
characteristic  of  a  random  quantity.  It  is  often  convenient  to  use  a 
different  characteristic.  In  the  case  of  a  unidimensional  random 
quantity  a*  the  concentration  or  scattering  can  be  characterized  by  the 

scattering  interval  (a^+v^^  ,  .  If  we  replace  thfe  true 

distribution  a*  by  a  uniform  distribution  over  the  scattering  interval, 

the  value  of  the  first  and  second  moment  of  the  distribution  remains  the 
same.  Analogously,  in  case  of  a  two  dimensional  random  quantity  (ag,  3g), 

its  concentration  about  the  "center  of  gravity"  of  a^,  3q  is 

characterized  by  the  scattering  ellipse 


I  r(«s-».)’ 


(4.39) 
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where  p  is  the  correlation  coefficient 


«.=(«;).  p. =<?;>. 


and  6  ^ , 

a’ 


are  the  second  central  moments, 

a3 


(4.40) 


(4.41) 

(4.42) 


If  we  replace  the  distribution  of  the  random  quantity  (a*,  3^)  by  a 

uniform  distribution  in  the  area  bounded  by  the  scattering  ellipse,  the 
value  of  the  first  and  second  moments  of  the  distribution  remain 
unchanged.  Obviously,  assignment  of  the  matrix  of  second  moments 
defines  the  scattering  ellipse,  and  conversely.  It  is  therefore 
desirable  to  have  analytical  formulas  for  the  dispersions  6^,  as 

well  as  for  the  second  mixed  central  moment  '5^^. 

We  first  assume  that  there  is  only  one  useful  parameter  a.  Then, 
according  to  the  above  definition  of  the  ensemble  x(t),  it  should  be 
assumed  that 


(4.43) 


where  a^,  are  fixed  values  of  the  signal  parameters.  The 

a  posteriori  probability  function,  or  its  equivalent  in  the  sense  of 
information,  the  output  Y(a),  is  random  in  this  case  and  can  be 
represented  as  the  sum  of  the  mathematical  expectation  and  a  random 
function  with  null  mean 

y(a)  =  (K((i))-j-r(a).  (4.44) 
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Here  and  below  the  superscript  zero  indicates  deviation  of  the  function 
from  the  mathematical  expectation 


y(a)=r(a)-(K(a)>. 


(4.45) 


In  accordance  with  the  maximum  likelihood  method  the  estimate  a*  is  found 
from  the  eqxiation 


(4.46) 


In  (4.46)  the  conventional  notation  for  the  derivatives  has  the  usual 
meaning 

^yK)=r'K)=[^r 


and  replacement  of  the  argument  a*  with  a*  in  the  random  quantity 

Y^(a)  is  justified  by  the  practical  equality  of  the  statistical 

characteristics  of  the  random  process  Y^(a)  at  points  a*  and  which 

are  close  together.  In  addition,  the  process  Y^(a)  can  be  assumed 
stationary  in  all  of  the  practical  applications  examined  below. 

Equation  (4.46)  makes  it  possible  in  principle  to  represent  the 

estimate  a$,  which  is  random  in  nature,  as  a  function  of  fixed  values 

^  d  0 

of  the  parameters  a^,  Eq  and  the  random  quantity-^Y  (a^)  with  known 

statistical  characteristics,  after  which  the  dispersion  of  the  estimate 
a*  can  be  found. 

In  the  case  of  very  strong  signals  the  error  of  the  optimal  system  (a^-a^) 

can  be  assumed  small,  and  its  mathematical  expectation  can  be  assumed  to 
be  zero.  In  this  case  likelihood  equation  (4.46)  takes  on  the  form 
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after  which  we  obtain  this  expression  for  the  potential  dispersion  of 
the  estimates* 


(4.48) 


When  there  are  two  useful  parameters  a  and  3  it  is  assumed  analogously 
that  the  received  oscillation  contains  a  signal  with  fixed  parameters 
“o’  ^0’  ^0’  maximally  reliable  estimate  (a*,  3g)  is  found  from 

the  optimum  output  Y(a,3)  by  solving  the  system  of  likelihood  equations 


da  (“o  '  Po  ))  +  •^  Po)  —  0. 

(4.49) 

^  (®0  •  Po  ))  +  4^  (“».  po)  =  0. 


The  solution  of  system  (4.49)  makes  it  possible  in  principle  to  find 
expressions  for  the  estimates  a*  and  3S  in  the  form  of  a  function  of  the 

3  0  ^3  0^ 

random  quantities -g-^Y  (a^,  3q)  , -g-g- Y'^(aQ,  3q)  and  the  true  values  of 
the  parameters  a^,  3q,  Eq,  after  which  the  scattering  characteristics 
of  the  two  dimensional  random  quantity  can  be  determined. 


Here  and  below  the  second  central  moments  of  the  estimate  for  fixed 
values  of  the  signal  parameters  are  designated  by  the  symbols  6^,  6^,  5^- 
without  the  subscript  su.  a  p  ap 


In  the  case  of  very  strong  signals,  system  of  equations  (4.49)  can  be 
written  approximately  as  follows: 


t“o— “o)-;^  («o.Po))4-|irj(^  («o.Po)>+-^^“K.Po)=o- 


(4.50) 


We  then  obtain  in  explicit  form  the  following  expressions  for  the  deviation 
of  the  estimates: 


■3^  (y  (“ti  ft))  "3^  ( y  {“••  p«))  ~  dad}  ] 


(y  f  .))-^y*  («..  h)  -iS? 

■3^  (y  (“li  ft))  ■3|i  (y  («t>  Pi)>  —  j^'3^3j’  (y  (“«>  f»))  j 


(4.51) 


In  nearly  all  calculations  of  the  potential  accuracy  of  the  signal  parameter 
estimate  presented  below,  the  optimum  output  Y(a,B)  for  a  signal  with  fixed 
parameters  a^,  3q,  Eq,  4)0,  in  the  vicinity  of  the  estimate  a*,  3*  can  be 

represented  in  the  form 

Y  (a,  p)  =  e.a‘'r  (a  -  a.,  p  -  P,)  +  c»  (a,  P).  (4.52) 


where  the  function  'i'(a,3)  is  determined  by  means  of  (4.33),  with  a  some 
number,  and  e(a,3)  a  normal  and  normalized  random  process  such  that 


(0  (»i.  Pi)®  («..  P»)>  =  («i  —  «i-  Pi  —  Pi)- 


(4.53) 
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The  following  relationships  followed  directly  from  (4.53): 


(4.54) 

(4.55) 

'  •  '  p=o 

(4.56) 

By  substituting  (4.52)  and  (4.54)-(4.56)  in  (4.51),  after  averaging 
we  obtain  these  comparatively  simple  analytical  formulas 

'  .  8' =((%-«.] 

l‘)=  *  •  * 

^  ;2aM<'(0.0)l  1-^*’ 

(4.57) 

((Po  P«n“‘.®a»|W''(0,0)|  1  — * 

(4.58) 

»P  “«)  (Po  Po))  —  0)  *  —  ’ 

(4.59) 

where 

Kp{0.0)j* 

(4.60) 

<(0.0)9'"(0.0)‘  • 

Analogously,  when  there  is  a  single  useful  parameter  a  the  optimum  output 
can  be  represented  as 


^  (®)  =  eBa*'!'”  (a  —  a,)  +  A®  (“) 


(4.61) 


and  analytical  expression  (4.48)  can  be  transformed  as  follows; 


(0)1 


(4.62) 
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When  receiving  signals  with  known  intensity  the  parameter  in  all  of  the 
formulas  should  be  set  equal  to  1. 


The  mean  values ^6^), and  of  the  obtained  scattering  characteristics 

6^,  6^,  and  are  determined  by  statistical  averaging  of  the  latter  with 

3 


■  2 
■’a* 


respect  to  all  possible  values  of  the  signal  parameters  a^,  bg.  Eq,  <Pq* 

When  the  optimum  output  Y(a,3)  can  be  represented  in  form  (4.52),  the 
characteristics  6^,  6^,  like  the  characteristics 

only  upon  the  intensity  parameter  when  condition  (4.33)  is  satisfied. 
Furthermore, 


<».’  >=  I  <«j  )=  S  s; 

0  0 


(4.63) 


with  an  output  with  the  form  (4.52)  and  known  intensity  of  the  received 
signals,  when  p(e)=6(£-l),  equality  (4.5)  holds 

(S»)=:82.  (8’)  =  8;.  (8.V  =  S»,. 


The  method  presented  in  this  section  for  analyzing  the  error  dispersion 
can  also  be  used  when  the  output  Y(cx,3)  is  non-optimal.  This  capability 
will  be  utilized  in  Chapter  8. 

4.4.  Method  for  Constructing  Optimal  Circuits  for  Obtaining  Estimate 
and  Error  Signal 

A  direct  method  for  obtaining  the  best  estimate  is  to  form  the  optimum 
output  and  find  its  highest  point.  In  many  cases,  especially  when  the 
output  is  a  function  of  several  useful  parameters,  this  method  for 
obtaining  the  estimate  encounters  practical  difficulties.  It  is  desirable, 
using  some  of  the  useful  parameters,  such  as  the  parameter  a,  to  obtain  an 
estimate  which  is  as  close  as  possible  to  the  best  estimate  without  forming 
and  reproducing  the  output  for  all  values  of  the  parameter  a  in  the  interval 
(a  -a  a  +a  )  of  the  possible  values  of  that  parameter.  It  is  assumed  that 
$  m’  $  m  .1, 

the  fact  of  the  occurrence  of  a  signal  with  parameter  a  belonging  to  the 

aforementioned  interval  is  established  reliably  and  is  not  subject  to 
doubt,  and  that  is  the  position  of  some  fixed  point  of  the  system  (the 

temporal  location  of  the  gating  pulses,  the  angular  position  of  the  axis  of 
the  antenna  system,  etc.). 
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An  analogous  situation  occurs  for  the  error  signal,  whose  derivation  is 
a  rougher  result  of  radar  observation  than  the  obtaining  of  a  parameter 
estimate.  The  error  signal  d  is  usually  defined  as  a  number  which  must 
be  proportional  to  the  deviation  of  the  fixing  point  from  —  the 

true  value  of  the  parameter  a.  When  the  reception  is  done  against  the 
background  of  interference,  the  true  value  of>the  parameter  cannot  be 
found  nor,  accordingly,  can  the  deviation  from  that  value  be  determined. 
The  potential  capabilities  of  the  system  are  exhausted  by  the  capability 
of  forming  an  error  signal  d  which  is  proportional  to  the  deviation  of  the 
fixing  point  from  the  best  estimate,  in  the  capacity  of  which  can 

serve  the  maximum  likelihood  estimate  a* 


(4.64) 


We  shall  assume  that  the  proportionality  coefficient  k  in  (4.64)  is  a 
positive  random  quantity.  If  the  number  k  is  known  a  priori,  obtaining 
the  error  signal  is  obviously  a  problem  which  is  identical  to  the  problem 
of  obtaining  the  estimate. 


We  note  that  in  the  general  case  the  error  signal  on  the  interval 

a^+a^)  can  be  expressed  through  an  arbitrary  monotonic  function  (p  of  the 

deviation  ot.-a* 

<p  0 


(4.65) 


Assumption  (4.64)  instead  of  (4.65)  causes  no  loss  in  generality,  since 
the  error  signal  d  can  be  obtained  from  d  by  employing  the  appropriate 
one-to-one  transformation. 

The  present  section  presents  an  analytical  method  for  defining  optimal 
schemes  for  obtaining  the  estimate  and  error  signal  based  on  employing 
and  solving  a  system  of  equations  of  the  type  (4.38). 
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We  shall  assume  that  the  signal  contains  two  useful  parameters  a  and  3 
The  optimum  output  Y(a,3),  as  follows  from  (4.29)  and  (4.33),  can  be 
represented  as  the  series  (4.20) 


Y{a,p). 


(4.66) 


=  ‘I  S  (“  -  *0  )'(P  -  Po (/-!-'=  2«;  /./.  rt  =  0. 1 .2. . . .). 

V  ' 


where 


y"  (/+/)! 


d'+/ 


da'dfl 


(4.67) 


The  use  of  assumption  (4.33)  simplifies  the  illustration  of  the  method, 
but  is  not  obligatory. 

We  shall  assume  further  that  the  estimate  3§  is  obtained  in  some  way, 

and  that  it  is  required  to  find  the  best  estimate  of  the  parameter  a 
without  forming  an  output  for  all  values  of  that  parameter.  In  accordance 
with  the  above,  one  method  for  obtaining  the  best  estimate  of  one  parameter 
is  to  solve  a  system  of  two  equations.  The  choice  of  the  system  of  equations 
for  obtaining  the  estimate  can  be  made  in  various  ways.  We  shall 
ourselves  to  the  class  of  equations  in  which  one  of  the  equations  defines 
the  error  signal  (i.e.,  expresses  the  error  as  a  function  of  the  estimate), 
and  the  other  makes  it  possible  to  define  or  exclude  the  unknown  coefficie 
from  the  first  equation.  For  example,  we  can  use  the  following  as  our 
system  of  equations: 
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l.e.,  formation  of  the  optimum  error  signal  consists  of  multiplying  the 
received  oscillation  x(t)  by  the  function  3^(t,a^)  and  integrating  the 

product  obtained.  In  order  to  form  the  estimate  it  is  necessary  to  form 
a  second  number  S)  or  compose  a  second  equation 


(4.71) 


by  performing  analogous  operations  with  the  function  3"(t,a^). 

a  0 


The  problem  of  determining  an  optimum  scheme,  as  well  as  the  problem  of 
composing  a  system  of  equations  for  obtaining  the  estimate,  does  not  have 
a  single  solution.  Another  possibility  for  obtaining  the  optimum  error 
signal  and  best  estimate  is,  for  example,  to  form  two  values  of  the  output 
Y(a(j^+A,Bg)  and  Y(a^-A,8*)  for  the  points  a^+A,B*  and  a^-A,B*  offset 

S3mimetrically  with  respect  to  the  fixing  point.  The  optimum  error  signal 
is  then  obtained  by  composing  the  difference 


K(a^  A.  Po* )-  -  A,  pj  )*  2Ar;  (a  ). 


(4.72) 


while  an  additional  equation  for  obtaining  the  estimate  can  be  the  sum 


(4.73) 


Keeping  (4.66)  in  mind,  we  obtain  approximately 


(4.74) 


Consequently,  this  method  for  designing  optimal  systems  for  obtaining  the 
signal  parameter  estimate  consists  of  the  following.  Based  on  representing 
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the  optimum  output  as  a  function  of  the  maximum  likelihood  estimate 
(4.27)-(4.39)  a  system  of  equations  is  composed  whose  solution  is  the 
estimate.  The  number  of  equations  n  must  exceed  the  number  of  measured 
parameters  by  one.  When  obtaining  the  estimate  of  one  parameter,  the 
error  signal  as  a  function  of  the  estimate  can  be  employed  as  one  of  the 
equations.  Then,  returning  to  the  representation  of  the  optimum  output 
as  a  function  of  the  received  oscillation  x(t),  the  equations  constructed 
are  modeled  with  the  help  of  an  electrical  circuit.  The  circuit  forms 
n  numbers,  which  define  the  value  of  the  free  terms  in  the  system  of 
equations  and,  consequently,  its  solution.  The  numbers  formed  by  the 
circuit  are  input  to  the  decision  device.  The  algorithm  for  obtaining 
the  estimate  in  the  decision  device  is  determined  directly  by  the  system 
of  equations  selected,  and  may  be  very  simple,  especially  if  the  errors 
are  small  and  representation  (4.35)  can  be  used.  Practical  examples  of 
the  application  of  this  method  for  constructing  optimum  circuits  for 
obtaining  estimates  and  optimum  error  signal  formation  circuits  will  be 
presented  during  the  investigation  of  direction  finding  systems. 

4.5.  Potential  Sampling  Accuracy  of  Range  Coordinate  from  Stationary 
Targets 

As  our  first  application  of  the  theory  of  radar  signal  parameter  estimates 
we  shall  calculate  the  potential  accuracy  of  sampling  the  delay  of  a 
signal  reflected  from  a  stationary  target.  The  only  useful  parameter  is 
the  delay  time  T.  The  modulation  function  S(t-T)  in  this  case  satisfies 
condition  (4.33).  We  shall  assume  at  the  outset  that  the  parasitic  random 
parameters  (j)  and  e  are  absent.  Then  the  optimum  output 


is  the  high  frequency  function  with  slowly  varying  envelope  L[Y(t)]. 
Determination  of  the  maximum  point  of  function  Y(t)  by  solving  the 
likelihood  equation 

y'(<t)  =  0^  (•*,  — •<'*<‘‘*  +  0  (4.76) 

yields  ah  ambiguous  result  because  of  the  high  frequency  carrier.  The 
region  (Tq-E,  Tq+e)  in  this  case  designates  the  vicinity  of  the  point  Tq 

the  true  delay  time  of  the  reflected  signal.  Conversely,  we  can  assume  for 
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the  envelope  of  the  output  L[Y(t)]  that  in  the  region  (Tq-E,  Tq+e)  the 

probability  of  the  occurrence  of  more  than  one  maximum  is  near  zero.  For 
this  reason,  we  shall  initially  seek  the  maximum  point  of  the  function 
Y(t)  from  the  envelope  of  L[Y(t)],  'i.e.,  from  the  equation 


(4.77) 


after  which  it  will  be  possible  to  answer  the  question  as  to  whether 
the  high  frequency  carrier  can  be  utilized. 

Assuming  that 


x(t)  —  n  (t)  -{-  >5  —  \) 


(4.78) 


the  expression  for  the  envelope  of  L[Y(t)],  assuming  that  the  signal  is 
significantly  stronger  than  the  interference,  can  be  reduced  to  the  form 
of  (4.61)  on  the  basis  of  (1.98) 


i  [r  (t)l = /2, 6  (,)  = 
~V[l— 5- -<.)■]  +  / 2,  »(,),  . 


where  'I'(t)  is  the  correlation  function  of  the  modulation,  3^  is  a  parameter 

which  determines  the  signal  spectrum  width  by  means  of  formula  (1.23),  and 
6(t)  is  a  normal  random  process  with  correlation  function 

(0(T,)0(T,))  =  'r(t,  — T,).  (4.80) 


The  solution  of  likelihood  equation  (4.77)  is  determined  by  formula  (4.62) 


(4.81) 


or 
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(4.82) 


1  00 


{  |5'(0|W 

—00 


Formula  (4.81)  was  derived  by  Woodward  [8]  by  calculating  the  dispersion  of 
the  averaged  a  posteriori  distribution. 

If  the  mean  square  error  6^  is  significantly  smaller  than  the  period  of  the 
high  frequency  oscillations  of  the  signal 


i  ^  1  (4.83) 

^  /.  ‘ 


knowledge  of  the  envelope  estimate  makes  it  possible  to  reduce  the  extent 
of  the  vicinity  of  the  point  Tq  in  initial  equation  (4.76)  enough  to 

eliminate  the  ambiguity  in  the  solution.  In  this  case,  solving  initial 
equation  (4.76)  under  condition  (4.78)  yields 


_  1 


(4.84) 


Consequently,  if  condition  (4.83)  is  satisfied,  i.e.,  if  the  modulation 
principle  and  power  ratio  q^  are  such  that  the , envelope  coordinate  sample 

eliminates  the  ambiguity,  the  high  frequency  carrier  of  the  signal  can  be 
used  to  sample  the  range.  The  potential  sampling  accuracy  is  then 
expressed  by  formula  (4.84).  Condition  (4.83)  can  be  satisfied  for 
certain  radio  navigation  and  remote  control  systems.  The  reverse 
inequality 


I 


/. 


(4.85) 


holds  for  radar  systems,  and  the  potential  range  sampling  accuracy  is 
determined  by  formula  (4.81)  (4.82). 
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Thus,  even  if  the  initial  phase  of  the  received  signals  is  not  a 
statistically  independent  random  parameter,  and  contains  useful 
information,  this  information  cannot  be  used  in  real  radar  lines  to 
increase  the  accuracy  of  the  range  coordinate  sampling,  or  to  increase 
detection  reliability.  Therefore,  we  shall  assume  everjrwhere  below 
that  the  initial  phase  of  the  received  signals  is  a  random  quantity 
which  is  distributed  uniformly  over  the  Interval  0,  217.  This  assumption 
also  corresponds  to  actual  operating  conditions  of  radar  systems. 

The  potential  range  coordinate  sampling  accuracy  in  radar  systems,  as 
follows  from  the  above  formulas,  depends  in  the  final  analysis  only 
upon  the  spectral  intensity  of  the  interference  n^  and  the  total  "energy" 

00 

Jl  (0 contained  in  the  derivative  of  the  complex  envelope 

of  the  signal.  The  accuracy  is  determined  by  the  parameter  3^  when  the 

signal/noise  ratio  is  given.  Therefore,  in  contrast  to  the  threshold 
ratio,  range  coordinate  sampling  accuracy  depends  significantly  upon  the 
form  of  the  signals  and  is  the  higher,  the  wider  the  modulation  spectrum 
and  the  more  of  the  modulating  signal  energy  is  concentrated  in  the  region 
of  high  modulation  frequencies. 

We  can  pose  the  problem  of  determining  the  form  of  the  signal  which 
provides  the  highest  sampling  accuracy  for  the  parameter  x  for  a  given 
signal/noise  power  ratio  q^  and  assuming  that  the  signal  spectrum  stays 
practically  within  the  band  (f^-F  ,  f„+F  ).  Such  a  signal  is 

J(0=l(t,  r)S^  cos  cos  2n^.  (^-86) 


since  its  spectrum  approximates  the  form  6  (f-fQ+F^)-l-6  (f-f^-F^)  and  has  the 

p 

maximum  possible  "dispersion"  (ZirF^)'^  with  respect  to  f=fQ. 

However,  no  matter  how  high  the  frequencies  F^  function  (4.86)  is  not 
suitable  for  use  in  the  capacity  of  radar  signals. 

Some  limitation  of  the  area  of  application  of  formula  (4.81)  must  also 
be  emphasized.  This  formula  is  valid  only  for  signals  above  the 
threshold.  As  the  parameter  3^  increases  the  threshold  increases  slowly, 

but  monotonically.  Therefore,  for  a  given  Interference  level  it  is 

not  possible  to  increase  the  sampling  accuracy  without  limit  by  changing 
the  form  of  the  signal  without  changing  its  power.  Sampling  accuracy  can 
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be  increased  in  this  way  only  as  long  as  the  reliable  detection  condition 
2  2 

q  >q  ,  is  retained.  In  addition,  formulas  (4.81)  and  (4.82)  are  based 
^  —  thr 

on  representation  (1.24)  and,  consequently,  are  valid  only  for  those 
signals  and  those  intervals  T  for  which  the  approximate  equality 


holds. 


This  representation  is  not  applicable,  for  example,  to  the  frequently 
employed  Idealization  of  the  envelope  of  valid  signals  in  the  form  of 
square  pulses  with  correlation  function 


(4.87) 


In  this  case,  using  formula  (4.17)  for  the  mean  square  errors  and 
representation  (4.28),  assuming  e=e*=l,  we  obtain 


0  " 


dx 


dx 


(4.88) 


or 


2,* 


c4 

•^0 


where  is  the  amplitude  of  the  input  signal  pulse.  According  to  the 

above,  for  the  signals  in  question  the  dispersion  (4.88),  like  the 
dispersion  (4.81),  is  the  same  as  the  averaged  dispersion  of  the  estimates 
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The  potential  sampling  accuracy  for  square  pulsed  signals  depends  only 
upon  the  spectral  intensity  of  the  interference  and  the  signal  amplitude, 
but  not  upon  the  pulse  duration.  Expression  (4.88)  can  be  viewed  as  the 
limit  which  the  potential  sampling  accuracy  approaches  as  the  slope  of 
the  leading  edge  of  the  signal  pulses  increases  without  limit. 

In  order  to  obtain  the  final  .expressions  for  the  sampling  accuracy  of  the 
range  coordinate  of  stationary  targets,  we  must  also  allow  for  the 
occurrence  of  two  parasitic  parameters  in  the  received  signal:  the  random 
initial  phase  ip  and  the  random  intensity  £.  In  this  case,  the  optimal 
output  Y(t)  can  be  the  envelope  of  oscillation  (4.75).  Assuming 


x(t)  —  n  (t) (( —  Xj,  vfi) 


(4.89) 


and  performing  transformations  analogous  to  those  which  lead  to  (4.79), 
we  obtain 


y(T)  =  e.V^(x-x.)  +  /29&(x)  (4.90) 


and 


I 


(4.91) 


As  might  be  expected,  in  accordance  with  the  conclusions  above  the 
occurrence  of  the  random  initial  phase  of  the  received  signal  does  not 
increase  the  dispersion  of  the  samples  of  the  range  coordinate. 

As  regards  the  influence  of  random  intensity,  comparison  of  formulas 
(4.91)  and  (4.81)  permits  the  conclusion  that  the  potential  dispersion  of 
the  samples  does  not  depend  upon  whether  or  not  the  intensity  of  the 
received  signal  is  known,  but  rather  upon  the  intensity  of  the 

signal  which  is  actually  received. 

The  averaged  potential  dispersion  of  the  errors  <"6^^  is  obtained 
theoretically  by  statistically  averaging  with  respect  to  the  parameter 
Eq  the  dispersion  6^  obtained  assuming  that  the  parameter  is  fixed. 

However,  averaging  cannot  be  done  directly  with  respect  to  all  possible 


80 


values  of  of  the  obtained  quantity  (4.91),  since  all  of  our  formulas 

as  based  on  the  assumption  that  reliable  detection  occurs,  and  are 
consequently  valid  only  for  sufficiently  large  e^.  In  addition,  formula 

(4.91)  has  a  perfectly  clear  physical  meaning  —  it  represents  the 
potential  accuracy  for  received  signal  energy  of  eQ2Q2^  regardless  of 

whether  or  not  the  intensity  of  the  received  signal  is  known  while  the 
estimate  is  being  obtained.  This  formula  provides  a  sufficiently  good 
characterization  of  the  potential  capabilities  of  the  system.  Therefore, 
we  shall  not  concern  ourselves  here  and  below  with  seeking  intelligent 
methods  of  averaging  the  expressions  for  the  dispersion  of  the  estimates 
of  various  parameters  with  respect  to  e^. 

During  reception  against  a  background  of  random  interference  with 

spectral  intensity  depending  upon  frequency  (normal  correlated  interference) , 

the  analytical  expression  for  potential  accuracy  (4.91)  takes  on  the  form 


»Kfi  0 


(4.92) 


The  parameters  and  were  defined  in  §  3.10.  N_  is  the  spectral 

3K5  T  3KB  U 

intensity  of  the  interference  at  carrier  frequency  f^. 


Chapter  7.  Azimuth  Scanning  and  Target  Direction  Finding 

7.1.  Azimuth  Scanning  by  Flat  Rotation  of  Directivity  Pattern 

It  has  been  assumed  up  to  now  that  the  scanning  has  been  done  in  one 
fixed  direction.  We  shall  now  assume  that  scanning  is  also  done  by 
angular  coordinate  6  within  some  arbitrary  sector  ’  which  we 

shall  assume  for  definition  to  be  equal  to  (0,0  ).  We  shall  begin  by 

m 

studying  the  basic  scanning  method  —  that  of  flat  rotation  (or  rocking)  of 
the  pattern  of  the  antenna  system  within  the  working  sector  (0,  6jjj)  • 

The  (voltage)  directivity  patterns  of  the  transmitting  and  receiving  antennas 
will  be  assumed  to  be  real  functions  of  0  and  will  be  designated  Yj(6) 

Yjj(0).  The  direction  of  maximum  gain  of  both  antennas  corresponds  to  the 
angle  0=0. 


* [Chapter  3  not  included  in  this  translation.  Tr.] 
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The  objective  of  the  investigation  is  to  obtain  analytical  formulas  for 
the  threshold  signal  and  for  the  potential  accuracy  of  the  estimate  of  the 
angular  coordinate,  which  is  sampled  jointly  with  the  range  coordinate. 

In  order  to  do  this,  we  must  write  the  analytical  expression  for  the 
optimum  output  and  study  its  statistical  characteristics  in  the  vicinity 
of  the  estimate  (in  the  region  of  high  correlation)  in  a  manner  analogous 
to  that  employed  in  Chapter  5  for  signals  containing  the  two  useful 
parameters  T  and  The  problem  of  determining  the  optimum  azimuth  search 
mode  is  not  dealt  with.  It  is  assumed  that  the  trials  are  made  at  fixed 
Intervals  (the  signal  has  fixed  duration).  It  is  assumed  accordingly  that 
the  angular  scanning  rate  or  rotation  rate  of  the  directivity  pattern  is  a 
constant  quantity  equal  to  (o  deg/sec,  and  that  wt  is  the  current  angle  of 
revolution  (angular  position)  of  the  pattern.  We  shall  examine  a  single 
scan  cycle,  so  that  the  angular  position  of  the  target 


0  =  (!>/ 

At 


(7.1) 


where  t  is  the  instant  in  time  or  the  time  offset  corresponding  to 

A 

rotation  of  the  antenna  system  through  an  angle  equal  to  the  angular 
coordinate  of  the  target.  Then  the  function 


represents  the  signal  emitted  in  the  direction  0,  and  the  received 
signal  which  is  reflected  from  the  target  and  whose  coordinates  and 
velocity  correspond  to  the  parameters  x,  t^  and  $,  will  be  expressed 
by  the  function 


*Y,(“  -f  t  -  /)]f„  [<0  -  0)5  (/  -  X,  <p,  «>).  (7.2) 


In  writing  (7.2)  it  is  assumed  in  accordance  with  the  actual  operating 
conditions  of  radar  systems  that  the  target  irradiation  time  in  each  space 
scanning  cycle  does  not  exceed  the  coherence  integral,  and  that  the 
fluctuation  in  the  scanning  rate  u  caused  by  movement  of  the  target  is 
negligibly  small. 
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In  theoretical  studies  it  is  convenient  to  approximate  the  directivity 
pattern  of  radar  antennas  by  the  bell-shaped  function 


Here  and  are  the  respective  pattern  widths  of  the  transmitting 

and  receiving  antennas  taken  at  the  half-power  level,  and  k  is  the 
proportionality  coefficient. 

Then  the  generalized  directivity  pattern  of  the  transmitting  and 
receiving  antennas  also  represents  a  bell-shaped  function 


+  ^  — OITiiK^ 

jfe* 


-0^ 


exp 
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1  “II 


(7.4) 


with  pattern  width 


^1^11 


(7.5) 
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which  corresponds  to  an  effective  signal  duration  of 


T 


2 


B 

111  2 


(7.6) 


The  occurrence  of  a  time  offset  between  reception  and  transmission 
results  in  an  offset  of 


vt 


(7.7) 


between  the  moment  at  which  the  maximum  of  the  pattern  is  directed 
toward  the  target  and  the  moment  corresponding  to  the  maximum 
amplitude  of  the  received  signal,  as  well  as  a  reduction  in  the 
resultant  gain  of  the  antenna  system,  which  is  expressed  by  the  first 
exponential  factor  in  (7.4). 

In  order  to  make  it  easier  to  write  subsequent  formulas,  we  shall  assume 
that  the  time  offset  VT  in  (7.4)  is  small  and  can  be  disregarded.  Then 
the  reflected  signal  takes  on  the  form 

*j(/  — T,  <p,  (]>,  e)  =  «T[«>(/^-0]Re{S((— (7.8) 


where  y(0)  is  the  generalized  directivity  pattern  of  the  transmitting 
and  receiving  antennas.  The  simplification  in  writing  (7.8)  which 
comes  about  by  disregarding  the  quantity  VT  is  acceptable  in  most 
cases.  In  the  more  general  case  the  time  offset  VT  can  be  included  in 
the  parameter  t  ,  which  again  leads  to  (7.8). 

A 

The  reflected  signal  is  structurally  the  same  as  the  signals  studied  in 
the  preceding  chapters,  except  for  the  fact  that  the  time  position  of  the 
function  y,  which  plays  the  role  of  cutoff  function,  depends  upon  the 
azimuth  of  the  target.  For  this  reason,  further  analysis  duplicates  to  a 
significant  extent  the  analysis  presented  in  Chapters  5  and  6.  We  limit 
ourselves  here  to  reproducing  the  order  of  the  calculation,  dealing  in 
detail  only  with  singular  features  associated  with  the  occurrence  of  the 
additional  random  parameters  —  the  delay  time  of  the  cutoff  function. 
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The  signal  or  modulating  function  S(t)  figuring  in  (7.8)  can  be  any  of 
the  radar  signals  examined  above,  including  the  sum  of  some  number  of 
incoherent  elementary  signals.  If  S(t)  is  a  coherent  or  incoherent  sequence 
of  periodically  repeating  signals  (t),  the  modulating  function  in 

(7. 8)  will  appear  as 


CO  s 

S  (7.9) 

l=~00 


where  ^  is  a  quantity  uniformly  distributed  over  the  interval  |c| 

and  which  allows  for  the  absence  of  synchronization  between  the  rotation  of 
the  directivity  pattern  and  the  modulation.  However,  C  is  a  parameter  of 
the  emitted  signal  which  can  be  allowed  for  exactly  during  the  processing, 
and  therefore  must  be  viewed  as  a  known  non-random  parameter.  In  addition, 
for  coherent  signals  all  of  the  <f)^  in  (7.9)  must  be  assumed  to  be  equal 
non-random  quantities. 

The  likelihood  coefficient  for  the  parameters  T,0  and  $,  which  determines 
the  potential  capabilities  of  the  system  and  the  optimum  processing  circuit, 
is  expressed  by  the  same  formulas  which  were  derived  in  our  examination  of 
spatial  scanning  in  one  fixed  direction.  The  only  difference  is  that  the 
function  Y  which  figures  in  the  formulas  for  the  likelihood  coefficient 
presented  in  Chapters  5  and  6  are  now  in  the  general  case  functions  of 
three  parameters  and  have  the  following  appearance: 

—  for  coherent  processing  [formula  (5.12)] 

.  OQ 

Y{x,  0.  0)  =  z.|^  (7.10) 


and  for  incoherent  processing  of  a  small  number  of  repetitions  [formulas 
(6.11)  and  (6.15)] 


0. 

* 

k  —00 


(7.11) 
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We  shall  assume  below  that  a  fixed  signal  with  parameters  Tq,  SgCtot  q), 
Bq,  (j)Q  is  input  to  the  receiver.  Then  the  outputs  (7.10)  and 
(7.11)  can  be  represented  in  a  form  analogous  to  (5.26) 


>'(x,  0.<I>)  = 


Ve,|>P(-c_v  e-0o.'&-®.)|+W„(t.0.4.) 

—  within  high  correlation  region, 

LI/i,{t,  6,  »i),  <p)] 

—  outside  high  correlation  region, 


(7.12) 


where  N  (t,0,$)  and  n  (T,0, $,({)) 
cx  i 


are  random  processes  for  which 


0..  0.,  0..))  = 

=:2^*|'l^(t,--c..  O.-0,.a>,-<I>.)|i 

(^1  ('*11  Oj,  'Pi)rt, (Xj,  0,,  ?j)) — 

*»>  Oi  0»'  ^’a)|  [2Tt/j  (x,  *»)"f* 

+  *(*!  — •'t.  — ‘I'a.  ^  — U  +  ?.  — <P.K 


(7.13) 


(7.14) 


and  |'F(t,6,$)|  and  xC^jQ,^)  are  the  absolute  value  and  argument  of  the 
joint  (for  the  parameters  T,0,$)  complex  correlation  modulation  function 
^(t,0,$).  We  shall  provide  the  analytical  expressions  of  the  joint 
modulation  correlation  function  |’1'(t,6,$)|  entering  into  (7.12)  for  the 
most  interesting  practical  cases. 

In  the  case  of  a  coherent  system  for  processing  the  received  oscillations 


09 

I'F  (^.  0.  0)1  ^  (U)0  ^ (‘0/ -  0)  S(t)  I*  (7 . 15) 

—00 
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The  signal  S(t)  in  (7.15),  as  a  rule,  can  be  represented  as  the  sum  of 
harmonic' components.  For  example,  when  Doppler  signals  are  used 


S  (/)  =  cos  2nFt  = 


1  Q 

2  “Jo  c 


(7.16) 


and  direct  substitution  of  (7.16)  in  (7.15)  yields 

|»F(‘*,  0.  a>)|  =  |cos2«fxl|'F  (0,  ®)|-{- 
®-2F)l-fi|'F(0.  ®-i-2/=-)|. 


Here  and  below  'i'^(6,$)  is  determined  by  the  formula 


'r^(0,  4>)  = 


w 

j  Y(w/)  y  {«*>/  —  6) 

"00 _ 

00 

J 


(7.17) 


(7.18) 


and  is  the  joint  correlation  function  with  respect  to  the  parameters 
6  and  $  of  the  directivity  pattern  or  cutoff  function  Y(wt). 

For  a  bell-shaped  approximation  of  the  directivity  pattern 


Y(fl)  =  *,exp/— 2In2 
\ 


(7.19) 


Furthermore 
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I'V.  '»)i  =  <=>'p(-'»2f)exp(-^). 


(7.20) 


In  the  general  case,  when 


(7.21) 


employing  the  notation  from  §  5.9 


.  4’)  =  - 


wu 

— 00 

SlSr-(01*rf< 


(7.22) 


for  the  joint  correlation  function  of  a  signal  taken  over  a  single 
repetition  period,  as  well  as  the  notation 


(7.23) 


and  assuming  that  the  adjacent  spectralized  r(f-)  and  r(f+^ — ) 
practically  do  not  overlap,  we  find  M 


(7.24) 


In  the  working  region  T 
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(7.25) 


l'r(x.  0. 


For  an  incoherent  system  with  a  small  number  of  repetitions  the 
correlation  function  entering  in  (7.12)  is  expressed  by  the  formula 


I'lf  {-c,  0,  tt))l  =  I  J  If  (««/)  If  H  -  0)  SrJ.i  -  kTjX 

k  —00 


(7.26) 


Using  the  same  notation  and  assumptions  as  in  the  preceding  case, 
(7.26)  becomes 


a.  ®)|=pp  (0.  o)|)'Fr>,  i>)t: 


(7.27) 


Thus,  in  systems  which  search  (scan)  the  range  coordinate  T  and  the 
angular  coordinate  $  jointly,  the  joint  modulation  correlation  function 
|'1'(t,6,$)  1  ,  which  determines  the  potential  capabilities  of  the  system, 
represents  the  product  (or  the  sum  of  the  products)  of  the  two  functions 

and  'F  ,  one  of  which  Includes  the  parameter  T ,  and  the  other  —  the 
TM  Y 

parameter  0.  Furthermore,  the  function  depends  only  upon  the  shape  of 

the  system  antenna  pattern,  while  the  function  depends  upon  the  form  of 

the  signal.  Because  of  this,  the  potential  system  indicators  in  terms  of 
the  angular  coordinate  9  depend  only  upon  the  directivity  pattern,  and 
not  the  form  of  the  signal.  On  the  other  hand,  the  potential  range  coordinate 
indicators  depend  upon  the  form  of  the  signal,  but  not  the  directivity 
pattern. 

With  incoherent  processing  the  information  about  target  velocity  contained 
in  the  set  of  Initial  phases  of  the  high  frequency  oscillations  with 
different  repetition  periods  is  lost.  The  only  portion  of  the  velocity 
information  retained  is  that  which  is  included  in  one  repetition  period  of 
the  reflected  signal.  This  physically  obvious  fact  is  reflected  in  formula 
(7.25). 
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The  expressions  derived  for  the  joint  modulation  correlation  function 
make  it  possible  to  determine  the  region  of  high  correlation  in  the 
working  sector  of  space  (t,6,$),  after  which  the  threshold  ratio  can  be 
calculated  by  the  corresponding  formulas  presented  in  the  preceding 
chapters.  The  equivalent  number  of  orthogonal  signals  (or  number  of 
resolvable  signals)  can  be  assumed  to  be  approximately 

volume  of  region  of  variation  of  useful  parameters  t,6,$  2\ 

™  ~  volume  of  region  of  high  correlation  in  working  sector  of 
space  (t,9,$) 


or 


(7.29) 


where  m  ,  mg,  m^  —  equivalent  number  of  orthogonal  signals  (number  of 
resolvable  positions)  along  the  T,0,and$  axes. 

The  extent  of  the  region  of  high  correlation  A0,  which  is  a  measure  of 
azimuth  resolution,  can  be  estimated  approximately  by  the  correlation 
angle 

j  |»K(0.  a,  0)|*rf8.  (7. SO) 


For  a  bell-shaped  approximation  of  the  directivity  pattern  (7.19) 


and 


»r(0,  0,  0)  =  v  (e)  =  exp^-ln2-^^ 


(7.31) 


(7.32) 
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The  problem  of  determining  the  potential  accuracy  of  direction  finding 
by  flat  rotation  of  the  pattern  for  pulsed  incoherent  systems  is  solved 
in  [31]  with  a  number  of  additional  limitations,  specifically,  assuming 
that  the  range  coordinate  is  known.  The  determination  of  the  optimum 
output  in  the  region  of  high  correlation  which  is  provided  makes  it 
possible  to  estimate  the  potential  angular  coordinate  sampling  accuracy 
in  more  general  form.  j 

In  accordance  with  the  analytical  method  employed  we  shall  assume  that  the 
oscillation  input  to  the  receiver  contains  a  signal  with  fixed  parameters 
Tq,  9q,  Eq,  (})q  and  that  the  angular  coordinate  and  range  coordinate  are 

estimated  jointly.  We  shall  assume  that  the  signal  in  each  repetition 
period  does  not  have  velocity  resolution 


(7.33) 


We  shall  also  assume  that  for  coherent  systems  the  sampling  is  done  in  a 
channel  for  which  in  (7.25) 


'  M 


(7.34) 


or  in  (7.17) 


(7.35) 


Then  the  output  of  an  incoherent  system  with  a  small  number  of 
repetitions,  as  well  as  of  a  coherent  system  for  a  frequency  channel 
through  which  the  signal  passes,  will  have  the  same  appearance  (4.52)  in 
the  region  of  high  correlation 


Y  (X.  fj)  =  2q\  1  'F,  (0  -  8,)  'F^,  {x  -  x,)  1  +  ff,  (x,  0). 


(7.36) 


where  the  random  process  N^(t,0)  is  such  that 


(W.  (X..  0.)  N,  (X..  e,'))’  =  I  'F  (8.  -  0,)  «F^^  (x._x,)  j.  (7.37) 
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The  matrix  of  second  moments  of  the  best  estimate  6*,  T*  is  defined  by 
formulas  (4. 57)- (4. 60) 


t? 

' 

t(  a 

1 


K 

0 


29**o 


where 


p; 


—  'y  h) 

at*  Tvi  vv 


«=0. 


* 


de*  7^ ' 


lt=:0 


(7.38) 


(7.39) 


Thus,  when  the  angular  coordinate  and  range  coordinate  are  sampled 
jointly,  the  potential  dispersions  of  the  estimates  or  the  dispersions 
of  the  estimates  of  the  maximum  likelihood  6^  are  6g  the  same  as  when 

one  of  these  parameters  is  estimated  with  the  second  parameter  known 
exactly.  Accordingly,  when  the  range  coordinate  is  known  exactly  (and 
when  the  parameters  are  estimated  jointly  (t  and  0))  the  potential 
sampling  accuracy  of  the  angular  coordinate  for  coherent  systems  and 
for  incoherent  systems  with  a  large  number  of  repetitions  is 


(7.40) 


When  the  generalized  directivity  pattern  of  the  receiving  and  transmitting 
antennas  are  approximated  in  the  form  of  bell— shaped  function  (7.19),  we 
find 


21n2 


(7.41) 
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The  formulas  presented  in  the  present  section  can  also  be  used  to  obtain 
more  complex  formulas  which  determine  the  dispersion  of  the  joint  estimate 
of  the  parameters  T,9  and 

7.2.  Optimal  Processing  During  Azimuth  Scanning 

We  have  not  yet  touched  upon  the  question  of  optimum  processing  of  the 
received  oscillations  during  spatial  scanning.  Our  investigation  of 
potential  capabilities  has  been  based  on  the  assumption  that  an  optimum 
output  Y(t,9,$)  or  Y(t,9)  is  formed  if  the  system  does  not  have  velocity 
resolution. 

We  shall  now  discuss  the  simplest,  and  near-optimum,  method  for  forming 
the  output.  This  method  is  based  on  the  fact  that  two  parameters  —  range 
and  azimuth  —  are  determined  by  the  delay  time  of  the  modulating  and  cutoff 
functions,  respectively.  A  delay  search  can  be  made  by  means  of  optimum 
filtering  or  some  equivalent  method  of  carrying  out  an  optimum  linear 
operation. 

As  a  simple  example  we  shall  examine  a  system  designed  to  receive  signals 
reflected  from  stationary  targets.  The  optimum  output  of  such  a  system  is 
determined  by  expression  (7.10),  assuming  that  $  =  0 


6)  =  L 


J  a:  (/)  Y  («./  -6)2  (/_x  _  kT^)  dt 

—00  k 


(7A2) 


We  shall  assume  that  the  receiver  consists  of  a  linear  system  with 
impulse  response 


h it)  =  T (“0  S  (-  /  -  kTJ  (7.43) 

* 


and  an  envelope  detector.  In  order  to  make  it  easier  to  write  the 
formulas,  the  time  offset  in  (7.43)  is  omitted,  which  makes  it  possible 
to  satisfy  the  conditions  for  implementation  of  the  system.  The  voltage 
at  the  output  of  such  a  receiver  will  appear  as 


(7.44) 
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and,  if  we  introduce  the  notation 


and 


/  = /iT’jj -j- ^  ^ 


(7.45) 


(7.46) 


Comparison  of  (7.46)  and  (7.42)  shows  that  Z(nTj^+T)  is  the  optimum  output 
for  discrete  values  of  the  angular  coordinate  (equal  to  wnT^+wx) 


Z(0  =  2(rt7„  +  't)  =  l'('c.  (7.47) 

The  angular  offset  WT,  when  significant,  can  be  considered  exactly, 
since  it  is  a  known  quantity  which  is  determined  during  independent 
sampling  of  the  range  coordinate. 

An  analogous  result  is  obtained  in  systems  designed  for  receiving  signals 
reflected  from  moving  targets. 

The  optimum  filtering  method,  or  other  equivalent  methods  for  accomplishing 
an  optimum  linear  operation,  thus  make  it  possible  to  obtain  in  a  single¬ 
channel  system  an  optimum  output  which  is  continuous  with  respect  to  the 
parameter  x  and  discrete  with  respect  to  the  parameter  6  with  a  digitization 
interval  of  co  T 

In  accordance  with  the  general  principles  for  practical  realization  of  the 
potential  capabilities  of  a  system  in  terms  of  detection  it  is  necessary 
that  the  extent  of  the  region  of  high  correlation  be  greater  than  the 
digitization  interval.  In  order  to  realize  potential  capabilities  in 
terms  of  accuracy  in  estimating  the  angular  coordinate  it  is  necessary 
to  have  at  least  two  values  of  the  output  in  each  interval,  equal  to  the 
length  of  the  region  of  high  correlation.  In  other  words,  in  order  to 
realize  the  potential  capabilities  of  the  system  when  a  single-channel 
arrangement  is  used  to  form  the  output  Z(t),  the  duration  of  the  reflected 
signal  must  be  at  least  twice  the  repetition  period 
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(7.48) 


T>2T^. 


If  condition  (7.48)  is  not  satisfied,  the  potential  capabilities  can  be 
realized  in  a  multi-channel  arrangement  for  forming  the  output  Z(t)  with 
the  appropriate  time  offset  between  the  channels. 


7.3.  Phase  Direction  Finding  Method 

The  phase  and  amplitude  direction  finding  methods  are  used  extensively 
in  radar  to  estimate  the  angular  coordinate  of  a  detected  target.  These 
methods  can  also  be  used  to  form  the  error  signal  and,  in  principle,  for 
azimuth  scanning  as  well.  From  the  viewpoint  of  antenna  engineering, 
phase  and  amplitude  direction  finding  systems  can  be  designed  in  many 
different  ways.  We  shall  not  be  dealing  with  the  design  singularities 
of  antenna  systems  or  other  details.  In  our  explanation  of  the  methods, 
we  shall  present  the  most  obvious  ways  of  implementing  them. 

We  shall  begin  by  examining  the  potential  angular  coordinate  sampling 
accuracy  in  the  phase  direction  finding  method.  The  antenna  system 
consists  of  two  reflectors,  with  feeds  and  Vi.2  the  foci  (Fig.  7.1) 

and  separated  by  distance  D.  The  patterns  of  both  antennas  are  the  same, 
and  have  maximum  gain  in  the  direction  of  the  geometric  axis  8  . 


Fig.  7.1.  Antenna  arrangement  for 
phase  direction  finding 


In  order  to  simplify  the  formulas  below,  it  is  assumed  that  the 
angular  position  of  the  geometric  axis  0^  coincides  with  the  direction 

corresponding  to  the  origin  of  the  angular  coordinate  9,  i.e.,  6^  =  0. 

The  angular  coordinate  is  sampled  within  a  relatively  small  sector  |0|_<  0^, 

for  which  the  antenna  gain  can  be  assumed  constant. 

The  form  of  the  signal  3(t)  or  the  form  of  the  modulating  function  S(t) 
is  of  no  importance.  It  is  assumed  here  and  below  that  the  reflected 
signal  has  fixed  intensity  and  contains  only  a  single  repetition  period 
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with  duration  T.  The  influence  of  the  fact  that  the  reflected  signal 
contains  a  greater  number  of  repetition  periods  and  has  random  intensity 
can  be  allowed  for  easily.  This  issue  was  examined  in  sufficient  detail 
in  our  study  of  coherent  and  incoherent  reception  methods. 

The  signals  reflected  from  the  target  and  received  by  the  first  and  second 
antennas  can  be  written  in  the  form 


{t  —  0.  ?)  =  Re  |5(/— x) exp  (f~x) -f-/ 0-}-/<p  j | . 

5,  (/  —  -c.  0,  <p)  =  (/— t)  exp  (t—x)  —}  ^  O-f /Yj|. 


In  (7.49)  sin  9  is  replaced  with  6  in  accordance  with  the  assumption  of  the 
smallness  of  angle  0,  and  the  initial  phase  cj)  can  contain  a  term  which 
determines  the  relationship  between  the  received  oscillation  and  the 
target  velocity. 

The  investigation  is  done  for  internal  interference  in  the  receiving  channel  . 
Therefore,  interference  in  the  channels  of  the  first  and  second  antennas 
are  mutually  independent  of  random  functions.  Considering  that  interference 
is  added  mainly  during  the  amplification  and  frequency  conversion  of  the 
received  signals,  we  can  imagine  two  processing  circuits  or  methods.  In 
the  first  method  the  received  signals  3^^  and  82  are  amplified  directly, 

converted,  etc.  In  the  second  system  the  sum  and  difference  signals 


3^,(1  — X,  0.  <f)  —  S^{t  —  x,  0,  <p)  +  J,(/  — t,  0,  9)  = 

2  cos  0)  Re  {S{t  -X)  exp  [M,  {t  -  x)-\-}<f]}, 

—  8.  ?)  =  •»,(<  —  ■*.  8.  ?)  — — 8,  9)  = 

=  2  si  n  0  j  I  m  {5  (/  —  x)  exp  I/2ir/,  (/  —  x)  -j-  /9]}, 


This  limitation  has  not  been  introduced  previously. 
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are  first  formed  in  the  antenna-feed  system,  after  which  the  signals 
3^^  and  3^2  amplified,  converted,  etc. 


In  accordance  with  the  definition  of  the  parameter  estimation  problem, 
we  shall  assume  that  the  first  and  second  channels  receive  as  input 
the  oscillations  x^(t)  and  X2(t),  consisting  of  statistically  independent 

interference  n^^Ct)  and  n2(t)  and  the  signals  3^  and  32  (or  3^1  and  3^2> 

with  unknown  useful  parameters  T  and  6.  We  shall  assume  that  the  a  priori 
distribution  function  of  the  useful  parameters  is  uniform  over  the 
required  working  region  of  variation  of  these  parameters,  and  is  equal 
to  zero  outside  the  working  region. 


The  likelihood  coefficient.  y\(T,e),  which  differs  by  a  constant  coefficient 
from  the  function  of  a  posteriori  probability  of  occurrence  in  the 
received  oscillations  x^^  and  X2  of  signals  with  the  parameters  (t,6)  is, 

for  the  first  method. 


where 


A(x,  0)  =  exp(— 2<7*)/, 


J  (/—■«.  6 •  (7.51) 


00  ' 

=  i  J  s]{i.  0.  J  ((=1. 2). 


(7.52) 


In  writing  (7.51)  it  is  assumed  that 


_L 

2n 


■  ^  n  CO 

Lj=o  —00 


d<f  = 


2n 


j  exp  |c 


n  oo  , 

cos  [*{x)  —  fli  SI  0)dt 

-<-1  -a 


df. 
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where  is  the  argument  of  the  complex  quantity 


2  ?  X^(t)S^(t~x)dt. 
<-i-4o 


We  can  use  as  the  optimum  output  Y(t,9)  the  function 


5  —  0.  + 


-f  J  x,(/)s,(/  —  t,  0,  <{•)<//]  , 


(7.53) 


which,  assuming  that  we  are  dealing  with  an  ensemble  of  cases  in  which 
a  signal  with  parameters  (Tq,0q)  is  input,  transformed  into 


r  (X.  0)  =  4^;  |,'K  -  T,)l  I  cos  ^^(0  -  oj  1+ 


(7.54) 


where 


<l»c,('‘*.  0.)  +  »c.("..  0.)]l»c.("..  0,)  +  »c2(-'..  ».)])  = 


kD 


2»ir(x,-x.)cos^(6,-e.), 


(7.55) 


The  joint  modulation  correlation  function  for  the  parameters  T  and  6, 
like  in  the  other  cases  of  joint  sampling  of  the  range  coordinate  and 
angular  coordinate,  is  the  product  of  the  corresponding  correlation 
function 


W(x,  6)  =  'r(,)»F^(0). 


(7.56) 
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Furthermore,  as  was  shown  above,  in  order  to  calculate  the  potential 
sampling  accuracy  of  the  angular  coordinate  6  when  T  and  6  are 
measured  jointly,  it  can  be  assumed  that  the  range  coordinate  T  is 
known  exactly. 

Considering  that  the  correlation  function  of  the  signal  for  the  angular 
coordinate 


(7.57) 


has  a  repetition  period  of  A/D,  elimination  of  ambiguity  in  the  sample 
requires  that  the  working  sector  be  selected  on  the  basis  of  the 
condition 

(7.58) 

m  p 


The  optimum  output  Y(t,9)  has  the  form  of  (4.52).  Accordingly,  we  use 
formula  (4.57)  or  (4.62)  to  find  the  potential  dispersion  of  the  angular 
coordinate  samples 


(7.59) 


Using  the  second  processing  method,  in  which  the  sum  and  difference  signals 
are  formed  in  the  antenna  feed  system,  the  likelihood  coefficient  A(t,0)  is 
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A('e,  0)  =  exp(— 

'  X  Re  {5  (/-■ t)  e'"*'* d/  -f  si n  ^ 


0  J  .x,(OX 


—00 

00 


fl  J.^,(0X 


(7.60) 


If  we  examine  an  ensemble  of  cases  with  fixed  parameters  (x^,  6^)  of  the 

received  signal,  the  expression  for  the  corresponding  optimum  output 
Y(t,6)  becomes 


r(t,  o)=:89‘1'F(x-x.)| 


cos^  (0  — oj 


+  2<7lO„(x,  +  0)). 


(7.61) 


The  random  process  defined  by  (7.55).  Output  (7.61)  is 

the  same  as  (7.54),  to  within  the  coefficient  2  of  the  quantity  q^. 

Therefore,  the  dispersion  of  the  estimates  of  the  angular  coordinates 
will  be 


(7.62) 


Thus,  if  we  can  only  disregard  noise  in  the  antenna  feed  system,  in 
contrast  to  noise  in  the  other  components  of  the  receiving  circuit,  the 
dispersion  of  the  estimates  of  the  coordinate  9  is  half  as  large  in  the 
second  circuit. 
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7.4.  Amplitude  Direction  Finding  Method 

We  shall  examine  a  simple  circuit  which  implements  the  method.  We  shall 
asstime  that  the  antenna  system  consists  of  one  reflector  and  two  feeds 
which  are  offset  symmetrically  with  respect  to  the  focus  in  the  focal 
plane  such  that  the  directions  of  the  maxima  of  the  directivity  patterns 
Y  (0)  and  first  and  second  feeds  are  offset  accordingly  by 

an  angle  of  -a  and  +a  with  respect  to  the  geometric  axis  of  the  antenna 
system,  as  is  shown  in  Fig.  7.2  and  expressed  by  the  replationships 


Y,(0)  =  *f(O-fa).  Ya(0)  =  lf(9-“)- 


It  is  assumed  that  the  directivity  pattern  yO)  is  an  even  function  with 
its  maximum  at  6=0. 

Within  the  working  sector,  which  approximately  surrounds  the  region  |6|<a, 
the  signals  3^  and  82  received  by  the  first  and  second  radiators  have  the 

same  phase  and  differ  only  in  amplitude 


Siit  —  t,  0,  <p)  =  Y(®+a)'»(^  — ■«.  f)> 

— t,  0,  (f)  =  Y(0  — a) J(/  — X,  ip). 


Fig.  7.2.  Directivity  pattern  of 
radiators  and  H2  for 

amplitude  direction  finding 
method 


The  likelihood  coefficient  for  the  parameters  (T,e)  using  the  first 
processing  method  (without  forming  the  sum  and  difference  signals)  is 
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A  (t.  0)  =  exp  lY*  (0.+  a) -f  Y*  (O-a)jJx 

^  —00 

+  Y(0-a)  J  X,it)3it-x. 


where  q  denotes 


(7.66) 


We  then  compose  the  likelihood  equation,  assuming  that  a  signal  is 
input  with  parameters  (t^,  0  )  and  that  the  range  coordinate  x.  is 
known  in  advance  ^ 

'  ~i^I“Y'('Jo+“)T(&;Tf  «)- 

“  Y'  a)  Y  -1-Y'  (O^N  Y  {«o+“)  +Y'  Y  (0,-a)l  + 

+  Y^if  Ilf '  (^0 + ®)  <'ci  4- Y' (®o  —  «)  ^cjl  =  0- 

We  expand  the  functions  yCTq-Kx)  and  YCT^-a)  in  the  vicinity  of  the 

estimate  by  powers  of  the  small  parameter  (e*-0  ) .  Then  the  likelihood 
equation  becomes  u  ,  u 

{®o  ~f*B)  {(tT  {®o  +  ®)J’  -f-  Iy"  (®o  ~  ?)]’}  + 

: :  +»75JKiT''(9;+«)i-+ir' (»;-.)]•  =0, 


where  0  is  a  normal  normalized  random  quantity.  Assuming  that 


[f  (<>0  +  a)]*  +  (Y'  (00  -  a)]*  2  (f  (a)]*. 


(7.68) 


102 


we  obtain  this  expression  for  the  error  dispersion: 


I  fY{“) 


V  h'w 


(7.69) 


With  the  second  processing  method,  in  which  the  antenna-feed  system  forms 
the  sum  and  difference  signals, 

•>11  —  ¥)  =  (Y(0+®)+V(®— f). 

'*12  ((  —  '*<  ¥)  =  [Y(®  +  «)  —  Y(®—®)]J  ((—'*.  9)- 

The  likelihood  coefficient  is 

A  (t .  8)  =  exp  ( lY’  (0 + «)  +  Y’  (0  -.4] )  X 

+  {Y(0-fa)-W-a)}  I  jc.COjC/--*,  'P)d/]). 

—09 

If  we  assume  that  the  received  oscillations  x^  and  contain  a  signal 
with  parameters  Tq  and  0q,  the  argument  of  the  function  in  (7.71)  can 
be  transformed  as  follows; 

{lY  (0  +  «)  +  Y  (0  -  ®)1  lY  (0.  +  «) + Y  (0.  -  «)1 1  ('*— ‘.)1  + 

+  lY  (0  4- a)  -  T  (0  -  a)]  (t  (0, -f  a)  - 1  (0, -- a)l  1 'F  (t-t.)  1 +• 

-I-  J^{[Y  (0+») + Y  (07<»)1  <>.  (^) + Iy  (0+“)  -Y  (0-“)]  ('')}  == 

'  +-^Iy(0+“)Oc.(-)+y(0-«)»c2(")1. 


(7.70) 


(7.71) 
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which  again,  to  within  the  coefficient  /l  of  the  first  term,  leads  to  a 
likelihood  equation  of  the  type  (7.67).  Accordingly,  for  processing  in 
which  the  sum  and  difference  signals  are  formed 


2 _  1  flifLl*  (7.72) 

84*  [y'H  J  ■  ■ 


Like  with  the  phase  direction  finding  method,  if  the  noise  in  the  antenna 
feed  system  can  be  disregarded  as  compared  with  the  noise  in  the  other 
elements  of  the  receiving  channel,  the  use  of  the  second  arrangement  in 
which  the  sum  and  difference  signals  are  formed  makes  it  possible  to 
reduce  the  dispersion  of  the  estimates  by  a  factor  of  2. 

In  order  to  illustrate  the  formulas  derived  above,  we  shall  calculate  the 
potential  angular  coordinate  estimation  accuracy  for  the  first  processing 
method  for  an  antenna  system  with  a  bell-shaped  directivity  pattern 
(7.19)  and  with  a  pattern  width  (at  the  half-power  level)  of  A.  Here 


,2.  (7.73) 

I  ~  •I9»(4aln2j*' 

and  assuming  that  the  patterns  intersect  at  the  half-power  level  (a  =  4) , 
we  find 


X2_  A» 

^  V(2ln2)»  * 


(7.74) 


It  is  interesting  to  compare  the  accuracy  with  which  the  angular  coordinates 
are  sampled  in  the  phase  and  amplitude  direction  finding  methods.  We  shall 
assume  in  both  cases  that  the  field  intensity  of  the  reflected  signal,  as 
well  as  the  intensity  of  the  interference  in  the  receiving  channels,  are  the 
same  and,  as  Fig.  7.3  shows,  the  total  aperture  area  of  the  antenna  system 
employed  in  the  phase  method  ($)  is  equal  to  the  aperture  area  of  the  antenna 
employed  in  the  amplitude  method  (A).  Then  the  total  power  of  the  signals 
received  by  radiators  and  H2  are  the  same  in  both  systems.  Accordingly, 

the  power  ratios  q^  which  figure  in  the  expressions  for  potential  sampling 
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accuracy  are  the  same  for  both  direction  finding  methods.  Employing  the 
approximate  relationship  [29] 


ii^l,22^Jpad\, 


(7.75) 


we  can  write  instead  of  (7.74) 


(7.76) 


Fig.  7.3.  Comparison  of  potential 

capabilities  of  direction 
finding  methods.  Aperture 
areas  of  antenna  systems 
for  phase  ($)  and  amplitude 
(A)  direction  finding  methods. 


On  the  basis  of  (7.59)  and  Fig.  7.3  we  have 


(7.77) 


In  the  phase  direction  finding  method. 

The  comparison  above  is  approximate;  even  so,  it  indicates  definitely 
that  the  potential  angular  coordinate  sampling  accuracy  for  both  methods 
is  of  the  same  order  and  that,  for  a  given  reflected  signal  field  intensity, 
it  depends  basically  on  the  size  of  the  aperture  of  the  antenna  system. 
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7.5.  Optimal  Processing  Schemes  for  Amplitude  and  Phase  Direction 
Finding  Methods 


It  was  pointed  out  in  Chapter  4  that  in  order  to  obtain  the  best  estimate 
of  a  signal  parameter  there  is  no  need  to  form  an  optimum  output  for 
each  value  of  the  measured  parameter.  Considering  the  random  nature  of 
signal  intensity,  in  order  to  sample  the  parameter  6  it  is  sufficient  to 
form  two  values  of  the  optimum  output,  after  which  the  best  estimate  is 
obtained  by  solving  a  system  of  equations.  In  defining  the  optimum  scheme 
we  shall  assume  that,  in  accordance  with  the  proofs  in  §  4.2,  the 
likelihood  coefficient 


A(t.  0.  6)  =  exp|-~r  J  j;  (/.  0.  9)^/4- 

—OO  ,  —00 

— QQ  — 


(7.78) 


in  the  vicinity  of  the  estimate  (t*,  6*,  e*)  can  be  represented  as 

A  (X,  0.  *)  =  exp|— J  0,  <p)rf^+  J<5®(/,  e,9)rf/jjx 

0,  — xj,  o;, 

'  *  ‘--00  ' 

CO 

+  j  0,  <?)«»,(/  — Xo,  Oo, 

—00 


(7.79) 


When  a  processing  scheme  is  employed  in  which  sum  and  difference  signals 

are  formed,  3„>]  and  must  be  substituted  for  3^  and  3.,  in  (7.78)  and 

(7.79).  2.1  2.^  12 

The  method  for  constructing  the  error  signal  and  estimate  formation  circuit 
was  described  in  §  4.4.  For  direction  finding  systems  this  method  consists 
of  the  following.  On  the  basis  of  (7.79)  operations  are  selected  which 
convert  the  function  A(T,0,e)  for  some  fixed  values  of  6  and  e  to  the  output 
d(T)  which,  for  a  value  of  x  equal  to  the  true  delay  time  Tq,  represents  a 

monotonic  function  of  the  deviation  of  the  direction  of  the  geometric  axis 
from  the  best  estimate  0*.  The  error  signal  d(TQ),  in  particular,  can  be 

proportional  to  the  deviation  of  the  geometric  axis  from  the  best  estimate. 
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In  order  to  obtain  the  estimate  6*  it  is  necessary  to  define  a  second 

conversion  of  the  function  A(r,0,  e)  to  an  additional  output  2) (t) , 
whose  knowledge  for  would  make  it  possible  to  eliminate  the  random 

intensity  parameter  from  the  error  signal.  We  shall  designate  the 
conversions  of  the  function  A  (t,6,£)  which  we  have  selected  by  the 
symbols  and  ^2’  the  fixed  values  of  the  parameters  0  and  £  in  the 

first  and  second  conversions  by  0^,  £^  and  02,  ^2*  respectively.  Then 

solving  the  system  of  equations 


1 ,  ( A  (t,  0, 


(7.80) 


makes  it  possible  to  obtain  the  estimate  0*  in  the  form  of  some  known 
function  ^  of  the  variables  ^(Tq)  and  £5  (t^) 

0o  =  9(‘fK).  ^K)l-  (7.81) 


The  receiver  must,  consequently,  normalize  the  output  d(T)  and  S3(t). 

The  operations  over  the  received  oscillations  x^(t)  and  X2(t)  and  the 

corresponding  radio  circuits  needed  to  do  this  _are  defined  on  the  basis 
of  representation  (7.78).  The  outputs  are  then  input  to  a  decision  device 
whose  functions  are  determined  through  (7.81). 

It  was  noted  in  §  A. 4  that  there  is  no  single  solution  to  the  problem  of 
composing  the  system  of  equations  for  obtaining  the  best  estimate.  A 
large  number  of  systems  of  equations,  and  accordingly  of  radio  circlets, ^ 
can  be  proposed  which  provide  potential  sampling  accuracy .  The  designer  s 
problem  is  obviously  to  select  the  systemof  equations  (7.80)  such  that 
the  corresponding  radio  circuit  has  certain  practical  advantages ,  such  as 
simplicity  of  implementation.  The  choice  of  systems  of  equations  (7.80) 
in  the  present  section  is  made  so  that  the  circuits  which  are  obtained  are 
as  close  as  possible  to  those  actually  used  in  practice. 
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In  the  phase  direction  finding  method,  employing  a  processing  scheme 
without  forming  some  indifferent  signals,  we  accordingly  use  the 
following  as  the  output  of  the  channel  which  forms  the  error  signal 
d(T)  and  additional  output  for  obtaining  the  estimate  (x) : 


d{t)  = 


^  (In  A  (IS.  0,  e))» 


!.=),  8-0' 


«=1.  8=0 


(7.82) 


We  shall  demonstrate  that  the  relationships  stipulated  above  between 
the  quantities  dCx^)  and^I^  (Tq)  and  the  estimate  are  obtained  in  this 

case.  If  we  consider  that  the  square  of  the  envelope  of  the  sum  of 
two  oscillating  processes  can  be  represented  as 

{L  [/,  (t)  cos  (mj  4-  cp,)  -f-  /,  (/)  COS  ((0,^4-  ?,)]}»  = 

=  /;  ii)  +  fl  (0  +  2f ,  it)  (t)  cos  -  <p,), 


after  substituting  (7.79)  in  (7.82),  we  find 

and  the  expression  for  obtaining  the  estimate 


arc  tg 


2nD  d  {■<.,) 

,  X  S5(x.)  ‘ 


(7.83) 


(7.84) 
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We  now  define  the  circuit  for  forming  the  outputs  (7.82),  for  which  we 
substitute  in  these  expressions  the  likelihood  coefficient  (7.78) 


A  jx,(0Re{S(/-t)e  -  "  }dt+ 

-<» 

2 


‘'W=5rj 

00 

p,(0Re{5((-t)e"' . "  }dt 

—00 

|2i',(.)y.wcos  W-«.(')+^  o]}^ 


}  = 


d 

d9 


(7.85) 


where  Y^(t),  Y2(t)  and  Xj^(t)  and  X2(t)  are  the  envelopes  and  phases  of 
the  high  frequency  processes 


'  00  ,  00 
^  ^x^{t)3{t  —  'z)dt  and  ^  x^{t)  3 {t  — 'i)dt,  (7.86) 


Thus, 


d  (x)  =  2  r,  (x)  Y,  (x)  sin  [Xj  (x)  -  *,  (t)]. 


(7.87) 


Analogously, 


gj  (x)  =  2  (^)  Y^  (x)  cos  [x,  (x)  —  X,  (x)l 


(7.88) 


and  on  the  basis  of  (7.84) 


(7.89) 
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The  optimum  processing  scheme  consists  of  optimal  filtering  of  the 
received  oscillations  x^(t)  and  X2(t)  in  linear  systems  with  identical 
impulse  responses 


h^{t)  =  h^{l)=zc3{—t,  9)  (7.90) 


and  measuring  the  phase  difference  between  the  oscillations  at  the 
output  of  the  optimal  filters.  Figure  7.4  shows  a  processing  circuit 
which  is  a  direct  interpretation  of  these  expressions.  The  notation  in 
the  diagram  is  as  follows:  FD  —  phase  detector,  FI  —  phase  inverting 
network,  which  multiplies  the  phase  by  Tr/2,  DD  —  decision  device,  which 
forms  the  estimate  0*  on  the  basis  of  (7.84). 


Hi 


Fig.  7.4.  Optimum  processing  circuit  for  phase  direction  finding 
method 


As  a  check,  we  shall  calculate  the  dispersion  of  the  estimate  (7.89) 
obtained  in  the  optimum  processing  circuit.  We  shall  assume  that  the 
input  is  a  signal  with  the  parameters  (x^,  9^,  e^,  (p^) .  Then  the 

oscillation  at  the  input  of  the  optimal  filter  will  be 
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(7.91) 


N» 


j  JC,  (/)  3  (t  —  2q\  1  vr  (t  —  t.)  1  X 

—00 


X  cospit/,  (t  —  -|-  <Po  +  X  ®»] 

+  ^cl  (■*)  +  ‘Po  +  X  ®*]"^ 

+  N^^  (t)  sin  2ic/.  (X  —  X,)  +  <p,H-  X  ®»]  • 


Since  the  signal  is  significantly  stronger  than  the  interference 
vicinity  of  Tq,  we  can  write 

in  the 

*1  ('*o)  —  ?o  +  X  .  2q‘u  ' 

(7.92) 

Analogously, 

»(  «)  T»  X  2q\  • 

(7.93) 

whence  according  to 

(7.89) 

-  (.'—A  1  X 

“  «  2v’«,  2nD 

(7.94) 

and 

®(  \(®o  ^o)  )  .,2/«0N*' 

(7.95) 

As  might  be  expected,  the  dispersion  of  the  estimates  obtained  in  the 
optimum  processing  circuit  is  the  same  as  the  potential  dispersion  of 
the  estimates  (7.59),  if  we  only  allow  for  the  random  nature  of  the 
intensity  of  the  received  signals  in  (7.59).  The  coincidence  of  these 
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dispersions  need  not  be  checked  below. 

In  the  amplitude  direction  finding  method,  employing  a  processing  circuit 
which  does  not  form  sum  and  difference  signals,  the  output  of  the  channel 
which  forms  the  error  signal  d(T)  and  the  additional  output  (t)  can  be 


and 


‘'W=?w 


,lnA('t,  0,  s) 


.-I.  8=0 


■^irci+oj+T’ (»-«))} 


»=l,l=0 


(7.96) 


(7.97) 


According  to  the  expression  for  the  likelihood  coefficient  in  the 
vicinity  of  the  estimate  (7.79) 


<i  w = -^  (T  (K +«)-•(('<:-«)).■ 


(7.98) 

(7.99) 


The  estimate  0*  can  be  found  by  solving  (7.98)  and  (7.99)  simultaneously. 

In  order  to  define  the  optimum  processing  circuit,  we  substitute  (7.78) 
in  (7.96)  and  (7.98).  Since  the  signals  in  both  channels  are  in  phase 


+  7(0  — a)  J  .v,(/)  — 


T(0+«)l',(^)+r(0-a)K.(x). 


(7.100) 
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where  Y^(t)  and  Y2(t),  as  above,  are  the  envelopes  of  processes  (7.86). 
Therefore, 


(7.101) 

(7.102) 


Figure  7.5  shows  the  corresponding  processing  circuit.  We  note  that  this 
circuit  does  not  respond  to  variations  in  the  initial  phase  of  the 
oscillations.  The  initial  phases  of  the  oscillations  of  the  received 
signals  Sj^  and  can  be  different,  or  even  statistically  independent. 

For  this  reason,  it  is  possible  for  this  circuit  to  employ  a  single 
processing  channel,  with  the  first  and  second  radiators  connected 
to  it  alternately  for  a  single  repetition  period. 


Fig.  7.5.  Optimum  processing  circuit  for  amplitude  direction  finding 
method 


In  the  case  of  processing  in  which  sum  and  difference  signals  are  formed, 
the  question  of  the  structure  of  the  optimal  receiving  circuit  will  be 
considered  in  more  general  form,  for  phase  and  amplitude  direction  finding 
simultaneously.  We  note  that  for  both  direction  finding  methods 


((  -'*.  0,  =  -0.  0,  «p), 

ii™  4r «■  = I  w  *  <'— •  *■  1„- 

i.*0  I  '•=<> 


0,  ?). 


(7.103) 

(7.104) 
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By  analogy  with  the  preceding  case,  the  output  of  the  error  signal  formation 
channel  can  be  the  function 


d{i)  =  c 


-^lnA('c,  (I,  e)  = 

•=|,J=0 

.“00 


“00 


0,  <p)rf/ 


(7.105) 


Further  transformation  of  (7.105)  considering  (7.104)  allows  us  to 
obtain 


c  lim 
a-«o 


20 


.f  (0*^11  ((  ■*>  9)dt 


00 

~f~  f -^a  (^)  ^)dt  L  J  — '*> — ^>v)d(-\~  • 

—00  J  _!oo  , 

00  •  '  li  ■ 

+  I  -e.  <f)dt  = 

“OO  J ) 

{00 

L  Umi  ^x^{t)s{t  —  x,  0,  i)dt-\- 
—00 

00 

+  j  x^{l)3\{t  —  x,  0.  ^)dt 
-00 

00 

ini  ^  r  x^{t)s{t  —  x,  0,,^)dt 
-  —00 

00  .  • 

J  X,{t)3\{,t  —  x,  0,  ^)dt 


—  L 


(7.106) 
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If 


CQ  W 

j  a:,  (0'S  (<  —  '*.  0.  (0 -sj  (<  —  ■'.  0.  ‘f)dt 


(7.107) 


(which,  of  course,  occurs  with  both  direction  finding  methods  near  1=10) > 

the  occurrence  of  a  large  coefficient  ahead  of  the  first  terms  in  the 
square  brackets  has  practically  no  influence  on  the  result,  since  these 
terms  compensate  one  another  after  the  operation  L  has  been  executed.  We 
can  therefore  write 


</(0  =  i  J  JC, (0'S (<  —  '*.  0.  + 

J  x^{t)3\{t  —  t.  0,  ?) d/  — 

J 

Ct  J Jf»{0J (/—'*.  0,  <i)dt 
— «> 
v9 

■Ct  f^*(0'»i(^  — 0,  <p)d( 


(7.108) 


where  c^  and  are  arbitrary  constants,  but  are  such  that  adding  them  as 

factors  to  the  right  and  left  parts  of  (7.107),  respectively,  does  not 
violate  the  inequality. 

We  assume  the  additional  output  ^  (x)  to  be 


^  (0  -r-Cjhn  A  (x,  0,  •)  -f-  ^  'Sji  -{-  (it  j 

-  00  .  - 

~  J  X^{t)3{t  —  x,  Oi  <p)d( 


9ss0 


1.  '.—00 


which  can  be  written  as 


(7.109) 
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«J(x)  =  Z.  L,  J  0,  <{))rf/+  ' 

L  —00 

”  1 


-\-L 


(7.110) 


—  c,  J  Xi(/)s'  (i~x,  0,  f) 


dl 


Figure  7.6  shows  an  optimum  processing  circuit  which  forms  (7.108)  and 
(7.110).  The  component  and  E^  shown  in  the  figure  form  the  sum  (+)  and 

difference  (-)  of  the  input  oscillations.  The  impulse  responses  hj^(t)  and 

h2(t)  of  the  linear  systems  must  be,  to  within  a  constant  coefficient: 

—  for  the  amplitude  method 

/«.  (/)  = /»,(()  =  3  (/.-/.  0. ,  ¥) 

—  for  the  phase  method 

/j,  (/)  =  Re  { S  (/,  —  t)  exp  {/  +  f)] } . 

/i,  (0  =  Im  { 5  (/,  —  0  exp  {/  (p))  | . 

The  functions  Y^(t)  and  Y2(t)  input  to  the  decision  device  (DD)  in 
Fig.  7.6  are  connected  with  the  outputs  d(T)  and£>(T)  as 


(7.111) 


(7.112) 


rf(x)  =  K.(x)-r.(x). 
23(x)  =  K,(x)  +  r,(x). 
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Fig.  7.6.  Optimum  processing  circuit  for  phase  and  amplitude  direction 
finding  forming  sum  and  difference  signals 


The  equations  for  obtaining  the  estimates  are  analogous  to  the 
corresponding  equations  in  circuits  which  do  not  form  the  sum  and 
difference  signals.  If  we  again  use  representation  (7.79)  for  the 
likelihood  coefficient  in  the  vicinity  of  the  estimate,  the  equations 
can  be  transformed  as  follows: 

—  for  the  phase  direction  finding  method 

i>(t,)=ccos^0‘. 


(7.113) 

(7.114) 


—  for  the  amplitude  direction  finding  method 


(7.115) 

(7.116) 


The  circuits  obtained  in  this  section,  with  the  possible  exception  of  that 
in  Fig.  7.6,  are  not  new.  But  even  so,  it  is  interesting,  first,  that  these 
circuits  were  obtained  by  applying  the  analytical  apparatus  for  defining 
optimal  circuits  and,  second,  that  it  has  been  proved  that  the  circuits  in 
Figs.  7.4,  7.5  and  7.6  exhaust  the  potential  capabilities  of  amplitude  and 
phase  direction  finding  systems. 
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7.6.  Multichannel  Systems 

Amplitude  and  phase  systems  are  designed  primarily  for  forming  an  error 
signal  and  sampling  the  angular  coordinate.  However,  they  can  also  be 
used  for  detection  as  well.  One  advantage  is  their  capability  of 
operating  with  stationary  antenna  systems.  Even  so,  they  have  serious 
deficiencies:  the  received  signals,  which  differ  only  in  angular 
coordinate,  are  strongly  correlated.  In  other  words,  amplitude  and 
phase  systems  have  no  azimuth  resolution.  Another  deficiency  is  the 
extremely  small  size  of  the  working  sector,  which  is  determined  for 
the  amplitude  method  by  the  width  of  the  directivity  pattern,  and  for 
the  phase  method  by  the  sector  of  indeterminacy,  the  angular  extent  of 
which  is  X/D  rad.  These  deficiencies  can  be  reduced  in  principle  to  some 
extent  by  employing  multi-element  and  multichannel  systems  [73,  74]. 


Fig.  7.7.  Directivity  patterns  of  radiators 

multichannel  amplitude  direction  finding  method 


With  the  amplitude  method,  a  multichannel  system  consists  of  some  number, 
such  as  2n,  of  radiators,  the  relative  positioning  of  whose  directivity 
patterns  y  (6) >  (6), . .Y2^(0)  shown  in  Fig.  7.7.  The  signals 

received  by  the  radiators  are  input  to  2n  independent  processing  channels. 
Each  received  signal  is  passed  to  the  output  of  two  channels.  The  angular 
coordinate  is  determined  by  the  numbers  of  these  channels  and  the  ratio  of 
the  amplitudes  of  the  signals  in  them. 

In  studying  the  phase  method  we  shall  again  examine  the  simple  example 
in  which  the  antenna  system  is  a  linear  array  (Fig.  7.8)  consisting  of 
2n  reflectors  with  radiators  H^,  H2,  their  foci.  The 

signal  with  parameters  (t,6,(())  received  by  radiators  ^2i-i  ^ 

be  written  as 
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(7.117) 


'*•  (if  ?) - 

=  Re  [5  (/  -  ti  exp  f  /2zf,  it  -  -)-f-  /  (2/  -  J )  -{-  ,>  j  j , 

— '*•  (*.  ?)  = 

=  Re  j5  (/  -  x)  exp  [ /2«/,  (/  -  x)  _  /  (2,-  _  | )  4. j  j . 


Fig,  7.8.  Antenna  ari'angement  for  multichannel  phase  direction  finding 
method 


If  the  received  signals  3^,  32,...,  82^^  are  processed  directly  in  2n 
independent  receiving  channels,  the  likelihood  coefficient 


A(t,  0): 


r  2/1 

’exp(— 2(?*/i)/,  — t.  0,  't)dt 

l.<=i  .  ' 


(7.118) 


assuming  that  a  signal  with  parameters  (x^,  0q)_ is  input,  is 


A  (^,  0)  =  exp  (-  2<7*/i)  /,  [2n2q*  1 V  (x  -  t.)  1 1  (0  -  0.)  1  + 

.  6)1. 


where  [30,74] 


(7.119) 


(7.120) 
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and  the  random  process  N(t,6)  Is  such  that 


(7.121) 


The  angular  coordinate  correlation  function  (7.120)  is  periodic,  with  a 
repetition  period  of  X/D.  Figure  7.9  shows  the  function  in  one 

repetition  period  for  n=l  and  n=4.  The  width  of  the  base  of  the  main 
lobe  is  X/nD.  The  side  lobes  of  the  function  are  relatively  small,  and 
drop  off  as  n  becomes  smaller. 


The  unambiguous  detection  angle  in  multichannel  phase  systems,  which 
determines  the  working  sector,  thus  does  not  depend  upon  the  number  of 
antennas,  but  rather  the  distance  between  adjacent  antennas.  On  the 
other  hand,  the  resolution  depends  upon  the  number  of  antennas  and  is 
determined  by  the  overall  extent  of  the  antenna  system  2nD.  The  number 
of  resolvable  angular  coordinate  positions  can  be  assumed  to  be  2n. 


Fig.  7.9.  Angular  coordinate  correlation  function  of  signal  for  phase 
direction  finding  method  (2n  —  number  of  antennas) . 


We  note  that  when  the  phase  direction  finding  method  is  employed  the 
analytical  formulas  for  the  angular  coordinate  performance  indicators 
are  similar  to  the  formulas  for  the  performance  indicators  for  the  range 
coordinates  in  Doppler  systems.  For  example,  the  correlation  functions 
'I'(t)  of  a  Doppler  system  with  2n  frequency  components  spaced  2F  apart  is 
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the  same  as  the  correlation  function  of  a  phase  direction  finding 

system  with  2n  antenna  systems  spaced  distance  D  apart  (Fig.  7.8).  The 
angle  D/X  rad  is  the  analog  of  the  frequency  interval  of  2F. 

The  potential  sampling  accuracy  of  the  angular  coordinate  is  determined  by 
formula  (4.62),  and  is  equal  to 


(7.122) 


where  [30] 


As  the  number  of  antennas,  and  accordingly  the  number  of  independent 
receiving  channels,  is  increased  the  accuracy  with  which  the  angular 
coordinate  is  estimated  increases  monotonically.  However,  when  the 
signal  energy  is  fixed,  like  in  the  case  of  expanding  the  signal 
spectrum,  the  accuracy  does  not  increase  without  limit,  but  only  until 
the  power  ratio  exceeds  the  threshold. 


If  the  noise  in  the  antenna  feed  system  is  significantly  smaller  than 
the  total  noise  in  the  receiving  channel,  processing  without  forming 
sum  and  difference  signals  does  not  fully  utilize  the  system  capabilities. 
The  loss  increases  as  the  number  of  channels  2n.  Let  us  suppose  that  the 
signals  3^^  2^  and  32i*  received  by  the  antennas' in  Fig.  7.8,  are  not 

input  directly  to  2n  processing  channels,  but  rather  form  2n  sum-and- 
difference  outputs  order  to  find  the  optimum  nature 

of  summation  of  the  signals  3^^,  •  -  - ,  shall  examine  the  result 


of  simple  summation  of  the  outputs  of  all  of  the  antennas  assuming  that 
the  signal  has  the  parameters  (Tq,  0q,  (JJq) 
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2n 

y]  *»,(/- s.  ^.  ¥.)= 

/=i 

ft  r 

=  2  Re  {5(/— t*)exp(/2K/,(/— x,)4-/>,)}  cos^yC!^/— l|f>. 


^2/m(/  — X.,  0,  f.) 


-  kD 

Sin  2/1  Y®» 

iiD&( 

2/1  tin  Y 


(7.124) 


The  summary. output  (7.124),  which  is  applied  to  the  input  of  one  of  the 
2n  processing  channels,  has  the  maximum  possible  intensity  when  9q=0. 

It  can  be  assumed  that  for  any  6^  belonging  to  the  working  sector  |6q|  _<  X/2D 

the  signal  intensity  will  be  close  to  the  maximum  possible  value  at  the 
input  of  one  of  the  2n  channels  if  the  angular  offset  between  the  directions 
of  maximum  intensity  in  the  adjacent  channels  is  X/2nD.  In  order  to  achieve 
this,  the  input  to  an  arbitrary  kth  channel  must  be  the  signal 


(7.125) 


where  k  —  —  n,  —  /t  -}-  1 2,  n  —  1 . 


Considering  (7.120)  and  (7.124),  the  expression  for  the  sum-and-difference 
signal  at  the  input  to  the  kth  channel  can  be  changed  as  follows: 
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(7.126) 


-5i*  (/  -  “to.  (»,.¥,)=  ; 

'*  , 

=  cos  j^ic  ^  (2/  —  1)  j  Co,  Oo.  <Po)+J2,_|  {/-  ‘'o.  ?o)l-H 

1=1  ■  ■  '  ,  ^  ■ 

+  ^‘"[’'2^(2'“ ')]K/(^-V  0..  ?o)l. 


which  indicates  the  manner  in  which  the  optimum  sun-and-dif ference  outputs 
are  formed  at  the  input  of  the  2n  processing  channels.  If  we  disregard  the 
occurrence  of  side  lobes  in  function  (7.120),  the  system  becomes  identical 
to  the  multichannel  amplitude  system  whose  directivity  patterns  are  shown 
in  Fig.  7.7.  The  width  of  the  base  of  a  directivity  characteristic  of  each 
channel  is  X/nD,  and  the  angular  offset  between  the  characteristics  of 
adjacent  channels  X/2nD.  The  diagrams  intersect  at  the  0.405  power  level. 

Finally,  if  a  single  processing  channel  is  employed  [73,  74],  to  which 
signals  are  applied  alternately  from  different  radiators  with  the  amplitude 
method,  or  different  sum-and-difference  signals  for  the  phase  method,  the 
system  becomes  equivalent  to  one  which  employs  flat  rocking  of  the  directivity 
pattern.  This  again  indicates  that  the  potential  capability  of  systems 
in  scanning  angular  coordinates  and  direction  finding  are  determined  basically 
by  the  geometric  dimensions  of  the  antenna  system,  as  well  as  the  number  of 
independent  processing  channels. 

In  order  to  simplify  the  task  in  this  chapter  we  have  excluded  certain 
capabilities  of  antenna  engineering.  For  example,  we  have  assumed  that 
the  (voltage)  directivity  pattern  y(9)  of  antenna  system  is  a  real 
function  of  0.  Furthermore,  the  correlation  function  'i’^(0),  which 

determines  the  potential  capabilities  of  the  system  with  respect  to  the 
angular  coordinate  0,  differs  little  from  the  function  y(0).  It  obtains 
accordingly  that  the  potential  capabilities  of  the  system  are  almost 
fully  determined  by  the  directivity  pattern  y(d)  itself  —  its  width  A. 
However,  in  the  general  case  the  function  y(®)  ^^7  complex,  and  the 
correlation  function  ’i'(0)  may  differ  significantly  in  width  from  y(0)* 

The  relationship  here  is  precisely  the  same  as  between  the  signal  S(t)  and 
its  correlation  function  ’i'(T)  [56]. 

Further,  during  our  examination  of  phase  (and  other)  systems  we  disregarded 
the  time  offset  between  the  complex  envelopes  of  the  signal  S(t-T) 
received  by  different  antenna  elements.  In  accordance  with  the  operating 
conditions  which  usually  occur,  the  latter  assumes  that  the  signal  correlation 
time  T  is  significantly  greater  than  the  time  intervals 

(~  D_^in_0j  corresponding  to  the  difference  in  travel  of  the  beam  to 
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different  antenna  elements.  However,  we  can  imagine  systems  in  which 
this  condition  is  not  met.  Then  the  joint  correlation  function  '}'(t;6) 
for  the  range  coordinate  T  and  angular  coordinate  9  cannot  be  represented 
as  the  product  ’I'(t)  'i'^(0),  and  the  performance  indicators  of  the  system 

in  terms  of  the  coordinates  T  and  0  become  interrelated  [56]. 

Finally,  it  would  be  helpful  to  examine  some  questions  involved  in  the 
practical  use  of  the  phenomenon  of  superdirectivity  [69]. 


Chapter  9.  Supplements 

9.1.  Signal  Reception  Against  Background  of  Interference  with 
Unknown  Intensity 

Up  to  now  we  have  been  studying  reception  against  the  background  of 
interference  with  fixed,  a  priori  known  intensity.  Reception  conditions 
such  as  these  occur  when  the  interference  is  caused  by  internal  radio 
equipment  noise.  We  shall  now  assume  that  the  interference  can  be 
approximated  by  Gaussian  white  noise  (§  1.2)  whose  spectral  intensity 
Nq  is  unknown  and  is  a  random,  but  constant  over  each  observation 

interval,  quantity 


(9.1) 


The  interference  intensity  parameter  Tl»  according  to  definition  (9.1), 
is  the  ratio  of  —  the  spectral  intensity  of  the  realization  of  the 

interference  in  a  given  observation  interval  —  to  the  mathematical 
expectation  of  that  quantity < Nq> .  We  shall  also  assume  that  the  received 

signal  e3(t-T,  (p,  $)  is  defined  by  formula  (5.10).  Then  the  likelihood 
coefficient  for  the  ensemble  of  all  possible  values  of  the  parameters  T 
and  $  can  be  expressed  as 


A(x,4»)=®- 


00  00 


Z*m  +  c*u*  \ 

r  r  •  .  r  1 

0 

■  1  )' 

CP(n)^  f  . 


T 


(9.2) 


where  p(e)  and  p(ri)  are  the  distributions  of  the  intensity  parameters  of 
the  signal  and  interference. 


1+T 


(9.3) 

t+r 

^  //I  (AT.)  1  ^ 

(9.4) 

% 


Integration  with  respect  to  t  in  (9.3)  and  (9.4)  must  be  done  over  the 
arbitrary  interval  which  includes  the  reflected  signal  3(t-T,(J),$) . 
Assuming  that  the  duration  of  the  reflected  signal  is  T,  the  integration 
interval  will  be  assumed  to  be  (t,  t  +  T) .  2m  is  the  number  of  measure¬ 
ments  of  the  interference  space,  which  is  determined  by  the  interference 
spectrum  width  f^(2m=2f^T)  and,  according  to  the  approximation 

employed  for  the  noise  background,  is  an  unbounded  large  (practically 
very  high)  number.  The  energy  signal/noise  ratio  is  the  ratio  of 
the  mathematical  expectation  of  the  signal  energy  to  the  mathematical 
expectation  of  the  spectral  intensity  of  the  interference 


(9.5) 

{NT)  * 


2 

The  number  Z  ,  as  will  be  shown  later,  is  defined  by  (9.4)  such  that  as 
m  ->  00  it  has  a  finite  limit. 

We  now  re-write  expression  (9.2),  transforming  it  somewhat: 


A(t.  4>)  = 


(9.6) 
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Asymptotic  estimates  of  the  integrals  in  the  numerator  and  denominator 
of  (9.6)  can  be  obtained  for  large  values  of  m  by  the  saddle-point 
method  [20].  The  idea  of  this  method  in  our  case  is  based  on  the  fact 
that  the  function 


(9.7) 


which  enters  into  both  expressions  beneath  the  integral,  has  a  single 
clearly  defined  maximum  at  the  point  n=z2.  The  larger  the  value  of 
the  parameter  m,  the  more  clearly  expressed  is  the  maximum  of  the  function 
[f(n)]™.  Therefore,  for  large  m  the  main  contribution  to  the  value  of 
the  integrals  comes  from  the  vicinity  of  the  point  of  the  maximum 
(Z^-h,  Z^+h) .  The  use  of  this  method  makes  it  possible  to  represent 
expression  (9.6)  as  the  ratio  of  the  functions  beneath  the  integral 
at  the  maximum  point  ir\=Z^) 


A (x.  •  J/> (e) exp (— j ^  (9.8) 

0 


As  m  increases  without  limit  the  sign  ~  of  the  asymptotic  estimate 
becomes  an  equals  sign. 

By  examining  formula  (9.8)  we  can  make  the  following  practical 
conclusions.  The  definition  of  an  optimum  system  does  not  require 
knowledge  of  the  a  priori  distribution  of  the  Interference  intensity 
p(ri).  The  likelihood  coefficient  A(T,<i’),  regardless  of  the  distribution 
of  the  signal  intensity  p(e)  for  every  given  value  of  the  number  Z^, 
is  a  monotonically  increasing  function  of  the  output  Y(t,'I>).  Therefore, 
the  optimum  decision  making  rule,  which  in  the  general  case  of  the  use 
of  the  Neumann-Pearson  criterion  is  expressed  by  formulas  (2.43)  and 
(2.44),  can  be  formulated  as  follows  in  our  case.  It  is  decided  that 
the  target  is  absent  if  for  all  possible  T  and  f 


(9.9) 
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It  is  decided  that  a  signal  is  present  with  parameters 


(t  ,  $  )  if 

X  X 


‘I',)  =  sup  IK  (t, 

t.  ♦ 


(9.10) 


2 

The  appearance  of  the  function  Y_(Z  ),  which  expresses  the  optimum 

^  2 

relationship  between  the  threshold  level  and  Z  ,  is  found  from  the 
equation 


=  jp(e)  exp(— 


(9.11) 


in  which  is  assigned  by  the  selected  false  alarm  probability.  On 

the  basis  of  (9.11)  we  can  write  the  optimum  form  of  function  Yq(Z^)  for 

the  distributions  p(e)  ordinarily  used.  For  a  fixed  signal  intensity, 
when  p(e)=6(e-l), 

A,  =  exp^— 

2  2 

and,  considering  that  q  »  Z  , 

K,^<7*  +  cZ*.  ...  '  (9-^2) 

In  the  case  of  Rayleigh  distribution  of  the  signal  intensities 

••••,.  '4' 

A  _  yg 

^0  —  z*  +  <7*  ^  P  [  4Z*  (Z‘  g') 

2  2 

and  considering  that  q»  Z  , 
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(9.13) 


The  constant  coefficients  c  figuring  in  (9.12)  and  (9.13)  are  determined 
by  the  given  value  of  Aq  or  the  given  false  alarm  level  F.  The 

formula  which  connects  the  coefficient  c  in  (9.13)  with  the  probability  F 
will  be  presented  below. 


One  feature  of  decision  making  rule  (9.9),  (9.10)  for  systems  with 
unknovm  interference  intensity  is  the  fact  that  the  threshold  level 
Yq,  like  the  output  Y(t,(})),  is  a  function  of  the  input  data  x(t).  Figure 

9.1  shows  an  optimum  processing  circuit  which  follows  directly  from  this 
decision  rule.  The  circuit  consists  of  two  channels,  the  first  of  which 
is  used  to  form  the  output  Y(t.'5).  The  structure  of  this  channel  was 
examined  in  sufficient  detail  in  the  preceding  chapters.  The  second 
channel,  which  forms  the  output  Z^,  incorporates  linear  system  LS2,  which 


is  a  wideband  (by  comparison  with  the  signal  spectrum)  amplifier,  a 
square-law  detector  and  a  low-pass  filter  (LPF)  with  cutoff  frequency 
of  approximately  1/T.  The  filter  provides  integration  over  the  interval 
(t,  t+T) .  The  outputs  Y(t,$)  and  are  input  to  the  decision  device, 
which  forms  the  optimum  threshold  level  Yq(Z'^)  and  makes  the  decision. 


There  are  some  possible  simplifications  to  the  structure  of  the  channel 
which  forms  the  threshold  level  Y-(z2).  In  order  to  clarify  these 

u  2 

capabilities  let  us  examine  the  analytical  expression  for  the  output  Z  . 


The  number  Z^  represents  the  result  of  extended  Integration  of  the  process 

n  , 

x^(t),  and  is  consequently  approximately  a  normal  random  quantity. 


Fig.  9.1.  Optimum  processing  circuit  for  signals  received  against 
background  of  interference  with  unknown  intensity 
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The  mathematical  expectation  and  dispersion  of  the  quantity  Z  can  be 
calculated  by  the  formulas  presented  in  §  1.2.  In  the  case  of  the  absence 
of  a  signal,  when  x(t)  =  n(t),  we  obtain 


where  6  is  a  normalized  normal  random  quantity,  and  tIq  is  the  realization 

of  the  noise  intensity  parameter  n  in  the  observation  interval  in  question 
In  case  the  signal  is  present,  when  x(t)=n(t)+eQ3(t-T, (}>,$) , 


2* 


(9.15) 


Assuming  that  the  number 


(9.16) 


is  extremely  large  when  the  noise  spectrum  f  is  extremely  wide, 
regardless  of  whether  the  reflected  signal  is  present  at  the  input, 
we  can  assume 


(9.17) 


We  now  note  that  the  output  Z  of  each  channel  in  Figi  9.1  retains  its 
structure  (9.14)-(9.17)  regardless  of  the  bandwidth  Af  of  the  amplifier 
in  each  channel,  if  only  TAf  is  large  enough,  e.g.,  of  the  order  of  10 
or  more.  In  addition,  the  square-law  detector  can  be  replaced  by  a 
linear  one,  since 


(9.18) 
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Keeping  these  remarks  in  mind,  the  optimum  processing  circuit  can  be 
modified  as  shown  in  Fig.  9.2.  Either  a  square-law  or  linear  detector 
can  be  employed.  This  circuit  assumes  that  a  linear  detector  is  used, 
and  that  the  intensity  of  the  reflected  signals  is  described  by  a 
Rayleigh  distribution.  Furthermore,  the  appropriate  choice  of  amplifier 
gain  K  ensures,  in  accordance  with  (9.13),  that  an  optimum  threshold 
level  is  formed. 


Fig.  9.2.  Simplified  optimum  circuit  for  processing  signals  received 
against  background  of  interference  with  unknown  intensity 
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We  shall  now  calculate  the  threshold  ratio  q  .  We  shall  assume  for 

nop 

simplicity  that  the  useful  signal  parameters  (t,$)  can  take  on  only  a 
single  value  (Tq,  ^g)>  i.e.,  we  shall  consider  a  two-alternative  system. 

The  false  alarm  probability  F  and  detection  probability  D,  keeping  (9.17) 
in  mind,  are  expressed  as 


?  PymdY, 
0  r/z’} 


D  =  lp(-n)drt  5  p,^{Yh)dY. 

6  rjv) 


(9.19) 


(9.20) 


in  which  Pj^(Y/ri)  and  Pgj^(Y/o)  are  the  distribution  functions  of  the 

output  Y=Y(Tq,  ii>g)  for  the  cases  of  the  absence  (N)  and  presence  (SN) 

of  the  reflected  signal,  assuming  that  the  spectral  intensity  of  the 
interference  is  fixed  and  equal  to  We  use  the  symbol  to  denote 

the  dispersion  of  the  interference  at  the  output  of  linear  system  LS^ 


130 


(Fig.  9.1  and  9.2),  assuming  that  ri=l»  Then 


(9.21) 


and  the  false  alarm  probability  for  a  given  fixed  value  of  the  intensity 
parameter  ri  is 

''(1)=  ?  (9.22) 


The  analysis  below  is  limited  to  the  most  interesting  case  in  which  the 
signal  intensity  parameter  e  is  described  by  a  Rayleigh  distribution, 
and  the  optimum  form  of  the  function  Yq(Z^)  is  determined  by  formula 
(9.13).  In  this  case 


(9.23) 


whence  it  follows  that  the  choice  of  the  optimal  form  (9.13)  of  the 
function  Yq(z2)  ensures  that  the  false  alarm  probability  remains  constant 

over  observation  intervals  with  different  values  of  the  realization  of 
the  quantity  n*  addition,  formula  (9.23)  makes  it  possible  to 
use  the  required  probability  F  to  select  the  quantity  c/a,  which  defines 
the  ratio  of  the  gains  in  the  channels  which  form  Y(t,$)  and  Yq. 

When  the  distribution  of  the  signal  intensity  parameter  is  described  by 
a  Rayleigh  approximation 


Ky  (^/v)  ,1  ^  exp  1^  ] 


(9.24) 


and  the  correct  detection  probability 

«0  n 

0 


(9.25) 
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or 


D-, 


exp[^ 


In  F 


*=2 


(9.26) 


We  can  conclude  from  examining  this  expression  that  when  the  detection 
probabilities  are  high,  on  the  order  of  0.9  or  more,  regardless  of  the 
distribution  of  p(Tl)»  and  especially  for  a  fixed  interference  intensity. 


(9.27) 


and  the  threshold  ratio  is 


(9.28) 


For  comparatively  small  values  (of  the  order  of  0.5  or  less)  the 
correct  detection  probability  D,  conversely,  depends  upon  the  distribution 
of  the  interference  intensity  p(ri). 

We  have  investigated  the  characteristics  of  reception  against  the 
background  of  random  interference  with  unknown  intensity  for  signals  of 
the  type  (5.10).*  Obviously,  the  method  used  to  calculate  the  likelihood 
coefficient  can  be  extended  to  the  case  of  receiving  signals  with  other 
sets  of  useful  and  parasitic  parameters.  Like  in  the  present  case,  the 
analytical  expression  for  the  likelihood  coefficient  for  systems  with 
unknown  interference  intensity  is  the  same  as  the  analytical  expression 
for  the  likelihood  coefficient  in  systems  with  fixed  interference 
intensity,  if  we  only  replace  in  the  latter  with  Z^<N3  .  Accordingly, 

by  substituting  Z  for  we  can  extend  all  of  the  results  and  formulas 

of  the  theory  of  detection  against  the  background  of  random  interference  with 
known  intensity  to  the  case  of  signal  reception  against  the  background  of 
interference  with  unknown  intensity. 


*[Chapter  5  not  included  in  this  translation.  Tr.] 
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9.2.  Use  of  Target  Information  Obtained  During  Previous  Observation 
Intervals 

All  of  our  investigations  have  assumed  that  the  duration  of  the  reflected 
signals  is  small  enough  that  the  received  oscillations  contain  information 
only  about  the  instantaneous  values  of  the  current  signal  parameters.  These 
signals,  of  comparatively  short  duration,  were  examined  in  isolation,  i.e., 
disregarding  the  information  contained  in  the  signals  reflected  from  the 
target  during  preceding  time  intervals.  This  statement  of  the  problem 
corresponds  to  the  actual  operating  conditions  of  modern  radars.  In  this 
regard,  Zibert  [15]  wrote,  "We  must  go  beyond  the  circle  of  ideal  associated 
with  finite  (and  consequently  short)  observation  time....  Every  radar 
target,  generally  speaking,  has  its  own  temporal  and  trajectory  history, 
which  is  a  more  thorough  expression  of  the  distinguishing  features  of  the 
target  than  the  instantaneous  values  of  its  current  coordinates  and^ 
velocity.  We  must  learn  to  use  this  target  characteristic.  More  difficult 
theoretical,  computational  and  practical  problems  arise  in  connection  with 
this  issue,  but  it  is  this  area  more  than  any  other  which  holds  hopes  for 
future  achievements  in  the  development  of  radar  engineering  . 

The  practical  solution  of  the  problem  of  more  or  less  total  utilization 
of  the  target  trajectory  in  order  to  increase  radar  operating  range, 
accuracy  and  resolution  is  an  extremely  difficult  matter.  In  the  present 
section  we  shall  suggest  that  a  decision  at  an  arbitrary  moment  in  time  T 
is  made  with  allowance  for  the  target  information  contained  in  the 
reflected  signal  for  all  of  the  preceding  observation  intervals,  and  we 
shall  discuss  in  general  features  the  extent  to  which  this  may  influence 
the  potential  capabilities  of  the  system  in  terms  of  range  and  detection 
accuracy.  The  extremely  important  and  interesting  question  of  the 
influence  of  this  processing  of  the  received  signals  on  the  resolution 
capability  is  not  dealt  with. 

For  definition,  we  shall  be  examining  a  periodic  scanning  detection  and 
tracking  system.  After  each  scan  cycle  (or  during  each  scan  cycle)  a 
statistical  trial  is  made,  as  is  a  decision  as  to  whether  or  not  a  target 
is  present  and  concerning  the  instantaneous  value  of  its  current 
coordinates . 

In  the  general  case,  target  coordinates  can  be  taken  to  mean  an  arbitrary 
set  of  useful  parameters  to  be  estimated.  Information  about  the  current 
coordinators  of  a  detected  target  is  produced  after  every  scan  cycle,  so 
that  the  time  variation  of  the  set  of  useful  parameters  —  the  target 
trajectory  —  is  determined  practically  unambiguously.  The  signal  duration 
in  each  scan  cycle  is  equal  to  the  target  irradiation  time.  The  only 
difference  of  this  system  is  the  fact  that  a  decision  is  made  during 
each  scan  cycle  on  the  basis  of  signals  reflected  from  the  target  within 
the  present  cycle  as  well  as  all  preceding  cycles. 
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The  signal  reflected  from  a  target,  considering  all  of  the  scan  cycles 
preceding  the  one  in  question,  called  the  kth  scan  cycle,  can  in  general 
be  written  as 


(9.29) 


where  s^  and  are  the  instantaneous  values  of  the  current  useful  and 

parasitic  parameters,  respectively,  during  the  ith  scan  cycle.  When  a 
decision  is  made  that  an  approaching  target  is  present,  the  number  k-1 
of  previous  scan  cycles  taken  into  account,  as  will  be  shown  below,  is 
not  important,  since  the  energy  (an  accordingly  the  information)  contained 
in  the  reflected  signal  for  distant  scan  cycles  is  negligibly  small.  The 
multidimensional  parasitic  parameters  u^  in  different  scan  cycles  will  be 

considered  independent  random  quantities  which  belong  respectively  to  the 
regions  U^.  The  functions  3^  with  different  values  of  the  subscript  i 

each  have  a  different  time  offset  by  the  corresponding  number  of  scan 
periods.  In  addition,  the  average  signal  intensity  varies  in  accordance 
with  variation  in  the  distance  to  the  target  during  different  scan 
periods. 

Since  we  agreed  not  to  deal  with  questions  of  resolution,  the  only 
useful  effect  provided  by  considering  previous  scan  cycles  is  the 
increased  signal  energy  realized  by  the  receiving  device  when  making 
a  decision.  We  shall  provide  a  quantitative  estimate  of  this  increase 
in  the  realized  energy.  We  shall  be  using  the  following  notation: 

^ —  range,  corresponding  to  the  moment  at  which  the  decision  is  made, 

Q  energy  of  reflected  signal  for  one  scan  cycle  with  target  range  of 
R»  AR  radial  component  of  length  of  path  covered  by  target  during  one 
scan  cycle,  a  —  relative  path  length 


(9.30) 


For  near-maximum  ranges,  the  quantity  a  is  small  —  of  the  order  of  0.01. 

Qj  total  energy  of  the  reflected  signal,  considering  all  preceding 

scan  cycles.  The  increase  in  realized  energy  can  be  estimated  by  the  value 
of  the  ratio  Qj,/q2,  Then,  assuming  that  the  energy  of  the  received  signal 
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is  inversely  proportional  to  distance  raised  to  the  fourth  power,  we 
find 


rfl  .  00 

^  ^  C  dz  _  \ 

Q*  2j  0  “  J  (<  3a  • 

*-0  .  0 


(9.31) 


The  energy  contained  in  sufficiently  distant  scan  cycles,  removed  from 
the  cycle  in  question,  e.g.,  by  more  than  1/ct  cycles,  is  negligibly  small. 
In  fact. 


I 


« 

s 

*->0 


I 

(1  +/!a)‘ 


1  <•  dz  _0,8]5, 
Tj  (1+2)*  3a  » 


(9.32) 


which  differs  little  from  (9.31) 

The  energy  gain  obtained  by  considering  previous  scan  cycles  thus 
increases  as  the  distance  to  the  target.  For  distances  close  to 
the  maximum  operating  range,  the  total  energy  of  the  received  signal 
over  all  of  the  scan  cycles,  or  over  the  closest  1/a  cycles,  is  several 
tens  of  times  greater  than  the  signal  energy  during  the  last  scan 
cycle.  It  can  be  expected  that  a  sharp  Increase  in  energy  will  lead 
to  a  significant  increase  in  range  (or  correct  detection  probability) 
and  parameter  estimation  accuracy. 

The  second  effect  which  must  be  dealt  with  when  attempting  to  employ 
previous  scan  cycles  in  the  decision  making  process  is  that  the  signal 
taken  over  several  scan  cycles  contains  more  information,  or  a  greater^ 
number  of  random  parameters,  than  the  signal  over  one  cycle.  Destruction 
or  extraction  of  additional  information  requires  the  consumption  of 
additional  signal  energy  and,  what  is  most  important,  it  complicates 
significantly  the  optimum  processing  circuit.  This  aspect  of  the  problem 
is  dealt  with  in  more  detail  below. 

The  a  posteriori  probability  for  a  signal  3  with  the  set  of  useful 

parameters  {s.}  =  s^,  . Sj^,  on  the  basis  of  general  principles, 

is 

(Q  33) 

P  •  •  »^*/^)  ^  ^xP  (®i»  •  •  •  >®*)  ^  (®>’  •  • 
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where  the  likelihood  coefficient 


(^if  •  •  S^)  —  ^  .  .J*  ^  (//,)  .  ,  ,  ^ 

u^  y* 

A  00  ,  . 

H  j*  ^  I 


(9.34) 


dUi.. .  dii^. 


/—I  —00 


The  coefficient  k  denotes  a  number  which  depends  only  upon  the  received 

oscillation  x(t),  and  which  can  have  a  different  value  in  different 
formulas.  The  likelihood  coefficient  taken  for  one  arbitrary  ith  scan 
cycle  will  be  designated  ■A^(s^) 


f  />(«<)  exp 

Ol 


((»  s^.  u^dt-\~ 


+4,  j 


s,,  «j)  dt  I  <//, 


(9.35) 


Then,  considering  that  the  functions  3^  and  3^  <j?^i)  do  not  overlap  in 

time,  and  also  that  the  a  priori  probability  p(s.,...,  s.)  can  be 
written  in  the  form  ^  ^ 


^  . h)  =  Z'  (Si . S*_,)  p  {sJSt,  .  . 


(9.36) 


we  find 

,  .  A(Sii  •  • '(S*)  =ri^/(s<)  —  A(Sif.v»s*_i)Aj^(s^)  (9.37) 

'  .  /=!  . 

I 

and 
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(9.38) 


y7  (S]  •  • .  |S^/x)  ^jiP  (®i  • 


•»S;^_l/x)  p  (sjh*  •  *  *»  ®*_l)  (®*)’ 


The  a  posteriori  probability  p(Sj^,...,  Sj^/x)  is  formed,  consequently,  by 

multiplying  three  functions:  p(s^,...,  the  a  posteriori 

probability  obtained  at  the  end  of  the  preceding  k-lth  scan  cycle; 
d(s  /s  ...  s  )  —  the  a  priori  conditional  probability  of  transition 
*  ^  Ic*  1*  *  *  *  lC“l 

of  the  object  to  the  point  s^^  after  it  has  covered  the  trajectory 

s  .  and  (s^^)  —  the  likelihood  coefficient  corresponding  to 

ordinary  optimal  processing  of  the  received  oscillation  x(t)  during  the^ 
kth  scan  cycle  taken  in  isolation.  Processing  of  the  received  oscillation 
during  the  kth  scan  cycle,  as  follows  from  (9.38),  can  add  nothing  to  the 
information  about  the  target  coordinates  ®k_2>  obtained  during 

preceding  scan  cycles.  The  purpose  of  forming  the  function  p(Sj^,...,  s^/x) 

is  to  ensure  that  the  most  reliable  possible  decision  is  made  about  the 
presence  of  the  target  and  about  its  instantaneous  coordinates  s^^  during 

the  kth  scan  cycle.  When  the  decision  is  made  at  the  end  of  the  kth  cycle 
all  of  the  other  coordinates  \_2»***»  ®i  subject  to 

estimation,  and  are  parasitic  random  quantities.  Accordingly,  when  k 
scan  cycles  are  used  in  the  decision  making  process,  the  optimum  output 
should  be  the  function 


P  P  (^tf  •  ’  ’  ^^k—l  — 

=KK  (s*)  1  5  •  •  -  J p^ . s*_,/x)  p  (s^/s,, . . .,  s*_,)  ds, . . 

[i  ^  j; 

or  any  other  function  which  corresponds  one-to-one  to  (9.39). 
region  of  variation  of  the  useful  parameter  s^. 


(9.39) 


S  is  the 


The  factor  in  the  square  brackets  in  the  latter  expression  is  the 
statistically  substantiated  prediction  A(Sj^)  for  the  current  scan 
cycle 

1(S^)=  f ..  .J  /7(s......s*_,/x)p(sjs,....,s*_,)rfs,...ds*^,= 

s  a 


(9.40) 


obtained  on  the  basis  of  observations  over  k-1  preceding  cycles  and  on  the 
basis  of  the  a  priori  conditional  probability  p(s  /s  , s,  -).  The 

K  i.  tC”l 

function  is  the  a  priori  probability  distribution  with  respect  to  the 

kth  scan  cycle. 

If  the  target  coordinates  s^  and  in  adjacent  scan  cycles  are  totally 

connected  random  quantities,  i.e.,  they  are  connected  uniquely  by  a  defined 
functional  relationship. 


Si  fi  (Sj) . Sj^_,  / (Sj), 

✓ 

(9.41) 

where  f ^ . 

•  »  are  certain  uniquely  defined  functions. 

Here 

(S' . s*)  =  />  (sj  8  [s,  -  ft  (s*)] ...  8  [s,_,  -  (s^)l 

(9.42) 

and 

i  (s*)  =  k^p  (S,)  A,  [/,  (s,)) . . .  A*_,  (s^j. 

(9.43) 

Complete  statistical  connectedness  of  the  target  coordinates  in  different 
scan  cycles  can  occur,  for  example,  when  all  that  is  to  be  established 
is  the  presence  of  the  target  on  a  given  trajecjtory,  or  for  a  practically 
stationary  target.  In  the  latter  case 

(9.44) 

(9.45) 

/‘)(s/x)  =  A^/>(s)A,  (s)...A,(s), 

(9.46) 
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which  leads  to  the  usual  operation  of  incoherent  accumulation  of  signal 
energy  over  k  repetition  periods,  where  the  optimum  output  Y(s)  can  be 
the  function 

K(s)  =  cVlnA,(s). 

/=! 


If  the  target  coordinates  in  adjacent  scan  cycles,  on  the  other  hand,  are 
statistically  independent  quantities,  then 


P(Si . S*)  =  />(S,).../?(S^), 

^(s*)  =  P(s*) 


(9.48) 

(9.49) 


and 


^(s/x)=A,p{s^)A^(s*). 


(9.50) 


As  might  be  expected,  if  the  target  coordinates  are  statistically 
independent  in  different  scan  cycles,  no  value  can  be  extracted  from 
the  information  obtained  during  the  preceding  cycles. 


Thus,  depending  upon  the  degree  of  statistical  connection  between  the 
target  coordinates  in  adjacent  scan  cycles,  the  potential  capabilities 
of  a  system  which  considers  all  preceding  scan  cycles  approach  one  of  two 
limits:  the  potential  capabilities  of  a  system  with  a  single  scan  cycle 
and  with  signal  energy  Q^,  or  the  potential  capabilities  of  a  system 
with  incoherent  accumulation  of  approximately  1/a  repetition  cycles  and 
with  total  signal  energy  of  q2.  More  detailed  study  of  the  question  of 

the  potential  capabilities  of  systems  employing  previous  scan  cycles 
requires  the  use  of  statistical  data  concerning  the  possible  tarpt^ 
trajectories  as  the  basis  for  studying  the  a  priori  joint  distribution 
function  p(Sj^,...,  s^)  of  the  sets  of  useful  parameters  s^ .  s.  tor 

i=l  2  ...,k.  The  absence  of  statistical  data  concerning  the  trajectories 
s  ,! .!,s.  of  the  reflecting  objects  makes  it  more  difficult  to  discuss 

the  question  of  resolution.  It  is  clear,  however,  that  the  presence  of 
the  additional  parameters  s^,...,  leads  to  additional  capabilities 

of  resolving  received  signals,  while  the  averaging  operation  (9.39) 
destroys  these  capabilities. 
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Figure  9.3  shows  one  possible  version  of  the  block  diagram  of  a  system 
which  makes  complete  use  of  the  information  over  k  scan  cycles,  which  is 
a  direct  interpretation  of  analytical  expression  (9.39)  for  the  a  posteriori 
probability.  The  notation  in  the  figure  is  as  follows:  C  —  i-system 
which  carries  out  operation  (9.35)  of  converting  the  received  oscillation 
x(t)  to  A^(s^);  X  —  a  multiplication  device  which  multiplies  the  input 

oscillations  (functions);  f  —  integrating  device,  which  carries  out 
operation  (9.39)  of  eliminating  the  parameters  s^^,...,  from  the  output. 

The  outputs  of  the  multiplying  devices  differ  in  several  coefficients  k 

X 

for  the  functions  p(s^,...,  s^/x)  shown  in  the  figure,  which,  of  course, 
does  not  violate  the  optimality  of  the  system  output. 


Fig.  9.3.  Optimum  signal  processing  circuit  which  realizes  target 
trajectory  information. 


This  optimum  processing  circuit  is  extremely  complex.  It  can  be  simplified 
significantly  if  we  avoid  the  attempt  to  reproduce  exactly  the  optimum 
processing  of  the  received  oscillations  over  a  large  number  of  scan  cycles. 
The  gain  achieved  by  considering  previous  scan  periods  will  be  smaller  than 
the  theoretical  maximum,  but  it  can  still  be  sufficiently  large. 
Investigation  of  these  capabilities,  as  well  as  the  a  priori  joint  distri¬ 
butions  for  the  sets  of  parameters  (s. , . . . ,  s,  )  are  beyond  the  scope  of  the 
present  study.  ^ 
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System  (4.68)  can  be  solved  for  a^*.  For  small  errors  and  small 

deviations  of  the  fixing  point,  when  representation  (4.35)  is  acceptable, 
we  can  use  the  approximate  values  of  the  right  part  of  equations  (4.68). 


In  this  case 


(4.69) 


Thus,  in  order  to  obtain  the  best  estimate  a*  it  is  necessary  to  form 
the  two  numbers  Y^(a^.  3q)  and  Y^(a^,  8*),  the  first  of  which 
represents  the  optimum  error  signal  d. 

By  way  of  illustration,  let  us  find  the  optimum  scheme  for  forming  the 
error  signal  and  obtaining  the  estimate  for  the  simplest  signal  e3(t,a) 
which  satisfies  condition  (4.33)  and  contains  the  useful  parameter  ot  and 
parasitic  parameter  e.  According  to  (4.68)  and  (4.23) 
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